Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 590,886 C→A G89G (GGG→GGT)  unknown ← Esterase ybfF (EC 3.1.‑.‑)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1590,8860CA98.7% 208.7 / ‑4.5 76G89G (GGG→GGT) unknownEsterase ybfF (EC 3.1.‑.‑)
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (29/46);  total (29/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

TTGTGGCATCAAAAGACTTCAGCATAAACTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑CCCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGCCGAATAAACCGTGAATAAAAACCAAAGTCGGCTTGTTAATTGCTTGTTTTAATTCATTAAAT  >  1/590639‑591129
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
ttGTGGCATCAAAAGACTTCAGCATAAACTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑accc                                                                                                                                                                                                                                                  <  1:141910/251‑1 (MQ=255)
      cATCAAAAGACTTCAGCATAAACTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCaa                                                                                                                                                                                                                                             <  1:842546/250‑1 (MQ=255)
         cAAAAGACTTCAGCATAAACTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATc                                                                                                                                                                                                                                         <  1:865975/251‑1 (MQ=255)
               aCTTCAGCATAAACTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACtt                                                                                                                                                                                                                                     <  1:328122/249‑1 (MQ=255)
                         aaaCTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATa                                                                                                                                                                                                                         >  1:861952/1‑251 (MQ=255)
                          aaCTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATaa                                                                                                                                                                                                                        >  1:816914/1‑251 (MQ=255)
                            cTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      <  1:490077/251‑1 (MQ=255)
                            cTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      <  1:552970/251‑1 (MQ=255)
                            cTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      >  1:184459/1‑251 (MQ=255)
                            cTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      <  1:858632/251‑1 (MQ=255)
                            cTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      <  1:675603/251‑1 (MQ=255)
                            cTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      <  1:329516/251‑1 (MQ=255)
                              gTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      >  1:735440/1‑249 (MQ=255)
                               ttGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGAt                                                                                                                                                                                                                    <  1:682383/250‑1 (MQ=255)
                                      cTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      >  1:746051/1‑241 (MQ=255)
                                      cTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGa                                                                                                                                                                                                             >  1:503411/1‑250 (MQ=255)
                                             cGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAac                                                                                                                                                                                                     <  1:51307/251‑1 (MQ=255)
                                              gTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      >  1:715395/1‑233 (MQ=255)
                                               tAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGa                                                                                                                                                                                                                     >  1:874515/1‑233 (MQ=255)
                                               tAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAaca                                                                                                                                                                                                    >  1:797389/1‑250 (MQ=255)
                                               tAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAaca                                                                                                                                                                                                    >  1:133287/1‑250 (MQ=255)
                                                 aaTATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAAt                                                                                                                                                                                                               >  1:179334/1‑237 (MQ=255)
                                                                   tGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCAt                                                                                                                                                                               <  1:774922/251‑1 (MQ=255)
                                                                      tttGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAg                                                                                                                                                                             <  1:755814/250‑1 (MQ=255)
                                                                       ttGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGtt                                                                                                                                                                           <  1:827991/251‑1 (MQ=255)
                                                                        tGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTc                                                                                                                                                                          <  1:332324/251‑1 (MQ=255)
                                                                                     ggTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAg                                                                                                                                                                                                                       >  1:575661/1‑193 (MQ=255)
                                                                                     ggTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAg                                                                                                                                                                             <  1:458994/235‑1 (MQ=255)
                                                                                     ggTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACtgg                                                                                                                                                             >  1:707954/1‑250 (MQ=255)
                                                                                     ggTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACtg                                                                                                                                                              <  1:543978/250‑1 (MQ=255)
                                                                                      gTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCCTAGAATG‑ACCCACCAAAATCACTTCGCGCAGAGATAAGGAGTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACtgtg                                                                                                                                                            >  1:344862/1‑251 (MQ=255)
                                                                                      gTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACtgtg                                                                                                                                                            >  1:594617/1‑251 (MQ=255)
                                                                                      gTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACtgtg                                                                                                                                                            >  1:560852/1‑251 (MQ=255)
                                                                                      gTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACtgtg                                                                                                                                                            >  1:354065/1‑251 (MQ=255)
                                                                                           aGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAg                                                                                                                                                        <  1:241167/250‑1 (MQ=255)
                                                                                                       tttCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGAt                                                                                                                                           <  1:482012/251‑1 (MQ=255)
                                                                                                       tttCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGAt                                                                                                                                           <  1:306471/251‑1 (MQ=255)
                                                                                                       tttCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGAt                                                                                                                                           <  1:284130/251‑1 (MQ=255)
                                                                                                       tttCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGAt                                                                                                                                           <  1:651143/251‑1 (MQ=255)
                                                                                                              gCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGc                                                                                                                                                                                          <  1:664070/197‑1 (MQ=255)
                                                                                                                 aaTAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATc                                                                                                                                 <  1:234305/251‑1 (MQ=255)
                                                                                                                      cccGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTAccc                                                                                                                            <  1:752315/251‑1 (MQ=255)
                                                                                                                         gTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTa                                                                                                                         >  1:515789/1‑251 (MQ=255)
                                                                                                                           tGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTaaa                                                                                                                       <  1:145046/251‑1 (MQ=255)
                                                                                                                             aaaaCCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAAt                                                                                                                     <  1:848601/251‑1 (MQ=255)
                                                                                                                                gcatccaattttacttGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATgg                                                                                                                   <  1:765762/235‑1 (MQ=255)
                                                                                                                                    catcatGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      <  1:638497/147‑1 (MQ=255)
                                                                                                                                        atGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATGACCCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATttcttc                                                                                                                                                                  <  1:243994/196‑1 (MQ=255)
                                                                                                                                        atGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTg                                                                                                                                             <  1:293074/216‑1 (MQ=255)
                                                                                                                                        atGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAAt                                                                                                                     <  1:739056/240‑1 (MQ=255)
                                                                                                                                             cATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAgg                                                                                                                                                                                                                      <  1:587056/138‑1 (MQ=255)
                                                                                                                                                 ttCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATg                                                                                                 <  1:118067/251‑1 (MQ=255)
                                                                                                                                                              aTGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCaataat                                                                                    <  1:467249/251‑1 (MQ=255)
                                                                                                                                                              aTGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCaataat                                                                                    <  1:709622/251‑1 (MQ=255)
                                                                                                                                                               tGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCG‑CCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATAc                                                                                  <  1:883522/251‑1 (MQ=255)
                                                                                                                                                                 ggCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATAc                                                                                  <  1:221646/250‑1 (MQ=255)
                                                                                                                                                                  gcgcAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATc                                                                                                            >  1:368227/1‑223 (MQ=255)
                                                                                                                                                                       atatCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAatt                                                                           >  1:242380/1‑251 (MQ=255)
                                                                                                                                                                       atatCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAat                                                                            <  1:620209/250‑1 (MQ=255)
                                                                                                                                                                           caatcaCCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTc                                                                       <  1:774886/251‑1 (MQ=255)
                                                                                                                                                                             atcaCCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCat                                                                     >  1:806037/1‑251 (MQ=255)
                                                                                                                                                                                accaccAGTCCTTTTTCCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATc                                                                  <  1:250139/251‑1 (MQ=255)
                                                                                                                                                                                 ccaccaGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATAccc                                                                                <  1:674438/236‑1 (MQ=255)
                                                                                                                                                                                 ccaccaGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCg                                                                 >  1:136125/1‑251 (MQ=255)
                                                                                                                                                                                  caccaGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGc                                                                <  1:517265/251‑1 (MQ=255)
                                                                                                                                                                                       gTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGCCGAAt                                                           >  1:767993/1‑251 (MQ=255)
                                                                                                                                                                                         ccTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGCCGAATaa                                                         >  1:841742/1‑251 (MQ=255)
                                                                                                                                                                                         ccTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGCCGAATaa                                                         >  1:638368/1‑251 (MQ=255)
                                                                                                                                                                                                        ggAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCaataat                                                                                    <  1:586696/209‑1 (MQ=255)
                                                                                                                                                                                                                   cgccAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACGAAAATCTCTTCGCGCAGTGATAAGGATTGAAt                                                                                                                                                                                                               <  1:374363/75‑1 (MQ=255)
                                                                                                                                                                                                                        gTACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATc                                                                  <  1:702995/211‑1 (MQ=255)
                                                                                                                                                                                                                         tACCATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGCCGAATAAACCGTGAATAAAAACCAGAGTCGGCTTGtt                           <  1:40430/249‑1 (MQ=255)
                                                                                                                                                                                                                            cATGGCGGTTTTGCCGCCCATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAac                                                                                                                                                                                                     <  1:260321/76‑1 (MQ=255)
                                                                                                                                                                                                                                            cccATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAACTCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGCCGAATAAACCGTGAATAAAAACCAGAGTCGGCTTGTTAATTGCTTGTTTTAATTCAtt      >  1:194134/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cccATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAg                                                                                                                                                        >  1:320968/1‑105 (MQ=255)
                                                                                                                                                                                                                                            cccATAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGCCGAATAAAc                                                       >  1:602693/1‑202 (MQ=255)
                                                                                                                                                                                                                                                tAGAATG‑ACCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGCCGAATAAACCGTGAATAAAAACCAGAGTCGGCTTGTTAATTGCTTGTTTTAATTCATTAAAt  >  1:545642/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
TTGTGGCATCAAAAGACTTCAGCATAAACTGTTGCACGCTTTCATCGTAAATATGTTGTGCCAAAATTGGTTTGGCTTCCTGACGGGTATGAGGCTGCGTGTTTTTCACCGCAAATAACCCGTTGAAAACCGCATCATGCCCATGTTCGCCGTAGGCAATGGGCGCAATATCAATCACCACCAGTCCTTTTACCAAATCAGGAGAATTCGCCGCCAGTACCATGGCGGTTTTGCCGCCCATAGAATG‑CCCCACCAAAATCACTTCGCGCAGTGATAAGGATTGAATGACATGAAACACATCTTCCGCCATCAAATCATAGTTCATTTCTTCACTGTGAAAGCTCGCGCCGTGATTGCGCAAATCTACCCGTAAAATGGCATAATCTTCACTAAATGCGCGGGCAATAATACCCAAATTATTCATATCGCCGAATAAACCGTGAATAAAAACCAAAGTCGGCTTGTTAATTGCTTGTTTTAATTCATTAAAT  >  1/590639‑591129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑