Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,456,758 A→T I162I (ATA→ATT)  unknown → ParA‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,456,7580AT98.1% 176.1 / ‑4.1 53I162I (ATA→ATT) unknownParA‑like protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (30/22);  minor base C (1/0);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

TGGTGCGCCACGTATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATTCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCAATTGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCCGGGTTAGTCATCAGTCAAACGAATAAGAAGACAAATTTAGCTAAAGAAGTTATTTCGTTTATCAACGAAAACTTTGAAATACCGTTATTAAAGGGTAGTACGGCGGTACGCGTCAGCTATGCAGAAGTCTTAAGTTCCGGTAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGCACTAGGAGTTCAATATGAGTAAATTAACAAAAGGTACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACGGCTATGCCACAGCTCAAAGA  >  1/1456521‑1456992
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
tggagctccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaag                                                                                                                                                                                                                               <  1:454616/238‑5 (MQ=255)
   agctccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggc                                                                                                                                                                                                                             <  1:756345/240‑7 (MQ=255)
       ccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggcagtac                                                                                                                                                                                                                        <  1:423631/245‑12 (MQ=255)
       ccccgcaTTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggcagtac                                                                                                                                                                                                                        <  1:806558/245‑12 (MQ=255)
                      gCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAg                                                                                                                                                                                                         >  1:227459/1‑251 (MQ=255)
                         aTGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTa                                                                                                                                                                                                      <  1:542422/251‑1 (MQ=255)
                         aTGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTa                                                                                                                                                                                                      <  1:480745/251‑1 (MQ=255)
                          tGACCGATTCGATTAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggc                                                                                                                                                                                                                             >  1:675372/1‑221 (MQ=255)
                                    gattAAACTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCt                                                                                                                                                                                                       >  1:397924/1‑239 (MQ=255)
                                         aaCTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggcagtacaacgg                                                                                                                                                                                                                   <  1:735347/222‑17 (MQ=255)
                                           cTTGCCGATTACATTCTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAgctc                                                                                                                                                                                    <  1:85244/251‑4 (MQ=255)
                                                          cTGATCCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGtttt                                                                                                                                                                     >  1:177639/1‑251 (MQ=255)
                                                                       aaCCATCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaggc                                                                                                                                                                                                                             >  1:291841/1‑176 (MQ=255)
                                                                            tCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   <  1:590249/251‑1 (MQ=255)
                                                                            tCACAATTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   <  1:683301/251‑1 (MQ=255)
                                                                                 aTTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGttattagaaagc                                                                                                                                                                                                                             >  1:437352/1‑166 (MQ=255)
                                                                                 aTTTGATATTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAAGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGa                                                                                                                                              >  1:425982/1‑251 (MQ=255)
                                                                                      atatTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGc                                                                                                                                          <  1:546383/250‑1 (MQ=255)
                                                                                      atatTTGGGCATGCAAAGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCa                                                                                                                                         >  1:404339/1‑251 (MQ=255)
                                                                                                      aGATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   >  1:49132/1‑225 (MQ=255)
                                                                                                       gattcgatagcgCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTaatcaaatca                                                                                                                         >  1:740660/12‑250 (MQ=255)
                                                                                                       gATTCGATAGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAgc                                                                                                                                                                                      >  1:33978/1‑188 (MQ=255)
                                                                                                           cgatagAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   <  1:242778/215‑1 (MQ=255)
                                                                                                                gAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAAt                                                                                                               >  1:237924/1‑251 (MQ=255)
                                                                                                                     ggTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCg                                                                                                                                               <  1:662222/214‑1 (MQ=255)
                                                                                                                       tGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAAc                                                                                                         <  1:8772/250‑1 (MQ=255)
                                                                                                                            gCACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAgc                                                                                                                                                                                      <  1:620220/168‑2 (MQ=255)
                                                                                                                             cACGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTa                                                                                                  <  1:525285/251‑1 (MQ=255)
                                                                                                                               cGTATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAgg                                                                                                >  1:200031/1‑251 (MQ=255)
                                                                                                                                 tATGCAAATTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGa                                                                                               <  1:96834/250‑1 (MQ=255)
                                                                                                                                      aaaTTAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCaa                                                                                         >  1:47557/1‑251 (MQ=255)
                                                                                                                                         ttAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAAGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAatat                                                                                      <  1:299934/251‑1 (MQ=255)
                                                                                                                                         ttAATGATAAGCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTAGTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAatat                                                                                      <  1:430768/251‑1 (MQ=255)
                                                                                                                                                   gCTCAAAGCGGGTTTAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGa                                                                                               >  1:530502/1‑232 (MQ=255)
                                                                                                                                                             ggtttAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   >  1:351420/4‑170 (MQ=255)
                                                                                                                                                             ggtttAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAggca                                                                  >  1:703895/4‑249 (MQ=255)
                                                                                                                                                             ggtttAGTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAggca                                                                  >  1:50307/4‑249 (MQ=255)
                                                                                                                                                                   gTCATCAGTCAAACAAATAAACAAACCAATTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTaa                                                            >  1:674597/1‑251 (MQ=255)
                                                                                                                                                                           tcaaacaaataaacaaaccaaTTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATcc                                                    <  1:335871/232‑1 (MQ=255)
                                                                                                                                                                                caaataaacaaaccaaTTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAgg                                                                                                >  1:471212/15‑202 (MQ=255)
                                                                                                                                                                                caaataaacaaaccaaTTTAGCAAAAGATGTTATTGAGTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCggg                                                 >  1:334441/15‑249 (MQ=255)
                                                                                                                                                                                                         agatgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   >  1:381207/13‑126 (MQ=255)
                                                                                                                                                                                                           atgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAg                                                                                              >  1:587947/11‑177 (MQ=255)
                                                                                                                                                                                                           atgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAAcag                     >  1:508724/11‑248 (MQ=255)
                                                                                                                                                                                                           atgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAAcag                     >  1:666918/11‑248 (MQ=255)
                                                                                                                                                                                                           atgttattgagTTTATTAAAGATAACTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGGATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAAcagc                    <  1:438021/241‑4 (MQ=255)
                                                                                                                                                                                                             gttattgagTTTATTAAAGATAACTTTGAAATCCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACAGCTa                  >  1:551538/9‑251 (MQ=255)
                                                                                                                                                                                                                         attaaagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTaaa                                                                               >  1:517246/10‑178 (MQ=255)
                                                                                                                                                                                                                         attaaagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACAGCTATGCCACAGCTCa      >  1:585959/10‑251 (MQ=255)
                                                                                                                                                                                                                            aaagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAaaga                                                                                                                                                   >  1:58835/7‑107 (MQ=255)
                                                                                                                                                                                                                            aaagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATggg                                      >  1:653384/7‑216 (MQ=255)
                                                                                                                                                                                                                             aagataaCTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACAGCTATGCCACAGCTCAAaga  <  1:334682/246‑1 (MQ=255)
                                                                                                                                                                                                                                    cTTTGAAATTCCGTTATTAGAAGGCAGTACAACGGTTCGCGTCAGCTATGCCGAAGTATTAAGCTCCGGCAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGTACTAGGAGTTCAATATGAGTAAATTAACAAAAGGCACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAAc                       >  1:835765/1‑223 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
TGGTGCGCCACGTATTGAAAAAGCCATGACCGATTCGATTAAACTTGCCGATTACATTCTGATTCCACTGAAACCATCACAATTTGATATTTGGGCATGCAAAGATTCAATTGAGCTGGTGCAAGCACGTATGCAAATTAATGATAAGCTCAAAGCCGGGTTAGTCATCAGTCAAACGAATAAGAAGACAAATTTAGCTAAAGAAGTTATTTCGTTTATCAACGAAAACTTTGAAATACCGTTATTAAAGGGTAGTACGGCGGTACGCGTCAGCTATGCAGAAGTCTTAAGTTCCGGTAATACGGTTTTTGAAAGTCAGTCTAAAGAAGCGAAAGCAGAGATTAATCAAATCACAAATGAAATTTTAACTGCACTAGGAGTTCAATATGAGTAAATTAACAAAAGGTACACTAAGCAAATCCGGGCTTTCTATGGGTAGCATTGAAAAAACGGCTATGCCACAGCTCAAAGA  >  1/1456521‑1456992

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