Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,844,2190TG30.3% 251.6 / 15.4 99F86V (TTT→GTT) adeDcryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (24/45);  new base G (30/0);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 6.99e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAATTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCA  >  NC_000913/3843985‑3844456
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
aaaTTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCTGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTAt                                                                                                                                                                                                                                              >  1:1708735/1‑239 (MQ=255)
       atcaCATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTTCTGGTGTTTGCGCAGAAAAAACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCCGGGGTtttttgtgc                                                                                                                                                                                                                                        >  2:2285827/1‑233 (MQ=255)
            aTTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCt                                                                                                                                                                                                                                  <  1:2285827/239‑1 (MQ=255)
              tAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGc                                                                                                                                                                                                                                <  1:1846127/239‑1 (MQ=255)
                 ccGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGGGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCCGGGGTTTTTGATGCCCACCTGccat                                                                                                                                                                                                                             >  2:1846127/1‑236 (MQ=255)
                     gCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTTTTTTTGCTCACCTGGCTAtttc                                                                                                                                                                                                                         >  2:819276/1‑237 (MQ=255)
                     gCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTGATTGTCACCTGCATATTGa                                                                                                                                                                                                                         >  1:1898518/1‑239 (MQ=255)
                                   aTTGATAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGAGTCCGGGGTTTTTTTTGCTCCCCCTCACATTTAATACCGCAggaaga                                                                                                                                                                                                           >  2:1384499/1‑239 (MQ=255)
                                     tGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATc                                                                                                                                                                                                                      <  1:1211559/226‑1 (MQ=255)
                                               ttCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTGTTGCTCACCTGCACATTTTATCCAGccctgtg                                                                                                                                                                                                            >  2:1211559/1‑220 (MQ=255)
                                               ttCCCGTGGCGACGCTGTTGCAGATTATATTATTGATCATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGTCATGGGTTTTTGTTTGTCACCTCGATATTTGATCCATCatgatgtagcccggc                                                                                                                                                                                                   >  2:1963706/1‑226 (MQ=255)
                                                      ggCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAAc                                                                                                                                                                                        <  1:1268677/239‑1 (MQ=255)
                                                                       ttatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAGTCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCcgcgc                                                                                                                                                                        <  1:1380685/238‑1 (MQ=255)
                                                                       ttatattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCcgcgc                                                                                                                                                                        <  1:1437598/238‑1 (MQ=255)
                                                                          tattatTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAAc                                                                                                                                                                                        <  1:1203837/219‑1 (MQ=255)
                                                                               ttGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCAGGGGTTTTTTGTGGTCACCCCCATATTTGATCCCGCATGGTGAAACCGGGCCCCTTTTTAACCCCCACCCCGCCGCCCgggcctgcc                                                                                                                                                               >  2:1437598/1‑232 (MQ=255)
                                                                                 gATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGGCTgacga                                                                                                                                                              >  1:2305864/1‑238 (MQ=255)
                                                                                 gATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAAATGCGGGGCCCGGGGTTTTTGATTCTCCCCCTCCTACTGAATCCCGGCTGATGAAGCCGGGCACTTTTTAAAACACTACCCTCCCCCCCgggcctgcga                                                                                                                                                              >  2:1825702/1‑230 (MQ=255)
                                                                                  aTAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCcgc                                                                                                                                                                          <  2:816702/225‑1 (MQ=255)
                                                                                  aTAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCcgc                                                                                                                                                                          >  1:816702/1‑225 (MQ=255)
                                                                                    aaTGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACg                                                                                                                                                                                                         <  2:65494/192‑1 (MQ=255)
                                                                                    aaTGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACg                                                                                                                                                                                                         >  1:65494/1‑192 (MQ=255)
                                                                                    aaTGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTTATTCTCACCCGCATATTGAATTCCGCATGGTGACCGCGGGCTCTTTTTAAACCCCTACCCCGCCGCCCGGccctgaga                                                                                                                                                              >  2:2358090/1‑229 (MQ=255)
                                                                                    aaTGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTTATGCTCACCTGCCTATTGAATCCAGCATGAATACGCCGGGCACTTTTGTAAACGCTAacccg                                                                                                                                                                              >  2:1203837/1‑217 (MQ=255)
                                                                                        tctcTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCa                                                                                                                                                                                                                                     <  1:1604463/160‑1 (MQ=255)
                                                                                         ctctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGaa                                                                                                                                                                                                                        <  1:664173/172‑1 (MQ=255)
                                                                                         ctctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGa                                                                                                                                                                                           <  1:573237/201‑1 (MQ=255)
                                                                                         ctctATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGTCAGGGGTTTTTGATGCTCACCCTCATTTTGAAACCAGCCTGATGCCGCCGGGCACTTTTTAAACCGCCTCACGGCCGCGCGGCCGGGAGACCGGTAttt                                                                                                                                                     >  2:5323/1‑239 (MQ=255)
                                                                                              ttCTCGACCTGATCAATGGCGGAGAAACTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTc                                                                                                                                                                                                   <  1:1322341/188‑1 (MQ=255)
                                                                                                  cGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTc                                                                                                                                            <  1:1316144/239‑1 (MQ=255)
                                                                                                  cGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCAGGGGTTATTGATGTTCACCTGCATAtt                                                                                                                                                                                                                           >  2:1604463/1‑160 (MQ=255)
                                                                                                   gACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTTTTTGTGCTCCCCCTCCTATTTAATCCCGCCTGATGTACCCGGTCACTTTTTAAAACGCCAcc                                                                                                                                                                                 >  2:573237/1‑201 (MQ=255)
                                                                                                   gACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCAGGGGTTTTTTATTCTCACCCTCCTTTTTAATCCCGCctga                                                                                                                                                                                                              >  2:664173/1‑172 (MQ=255)
                                                                                                        gATCAATGGCGGAGAAACTTCCGGCCCCATTGTGCTTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTTTTGTTCACCTTCATATTTTATCCAGCATGGTGACGCCGGGCCCttttttaa                                                                                                                                                                                         >  2:1322341/1‑184 (MQ=255)
                                                                                                         aTCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCAGGGGTTTTTGATGCTCACCCGCATATTTGATCCCGCATGATTACCCCCGTCTCTTTTGTAACCACCTCCCCGCCGCCGCGCCTGTCGACCGCTATTTTCGTCCCCCactcaag                                                                                                                                     >  2:2241794/1‑233 (MQ=255)
                                                                                                              tGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTTTTGATGCCCACCTCGCTATTTAATCCAGCAGGATTACCCCGGTCACTTTTTAAAACCCCACCCCGCCGCGCGGCCCTGCGACCCTTATTTTCGGCCCTCATGAAAACgccc                                                                                                                                >  2:911211/1‑236 (MQ=255)
                                                                                                                          ttCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGt                                                                                                                                                          <  1:2304798/201‑1 (MQ=255)
                                                                                                                                    aTTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCAGGGGTTTTTGATGCTCCCCTGCCTATTGAATCCAGGATGATGACCCCGGTCACTTTTGAAAACGCTACCCCTCCGCGCGGCCTGGCGACCGTTATTTGCggc                                                                                                                                                >  2:2304798/1‑199 (MQ=255)
                                                                                                                                               aGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATg                                                                                                                                                                                                                                         <  2:2381029/101‑1 (MQ=255)
                                                                                                                                               aGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATg                                                                                                                                                                                                                                         >  1:2381029/1‑101 (MQ=255)
                                                                                                                                                      tACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGa                                                                                                                                                 <  1:1876553/182‑1 (MQ=255)
                                                                                                                                                             ccGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCCACGGCGGGGGCAGGGGTTATTGTTGCTCACCTGCCAATTGGATCCAGCATGCTGAGTGCGGTCCCTTTTTAAAACGCTTCCCCGCCGCCCGGCCTGACGGCCGTTATTTTCTGCCCTACTGAAATCGTTCAAGGTAGGGGCGGAACCGGTTTTGCCCTGTTTGACCGccggtgcgtc                                                                              >  2:1400927/1‑230 (MQ=255)
                                                                                                                                                             ccGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGt                                                                                                                                                          <  1:1814088/166‑1 (MQ=255)
                                                                                                                                                                gtgtTGGCGCAGAATGCACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATg                                                                                                                                                                                                                                         <  2:1036234/84‑1 (MQ=255)
                                                                                                                                                                gtgtTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATg                                                                                                                                                                                                                                         >  1:1036234/1‑84 (MQ=255)
                                                                                                                                                                gtgtTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatg                                                                                                                                                                                                               >  2:1321888/1‑110 (MQ=255)
                                                                                                                                                                gtgtTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCatg                                                                                                                                                                                                               <  1:1321888/110‑1 (MQ=255)
                                                                                                                                                                gtgtTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTATTTATGCTCACCTGCATATTGAATCCCGCATGATGACGCCGGGCACTTTTTAAACCGCTACCCTGCCGCGCGGCCTGACGACCGGTTTTTGCGCCCCTCATGaa                                                                                                                                       >  2:1876553/1‑182 (MQ=255)
                                                                                                                                                                 tgtTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGc                                                                                                                                                                                                                                <  1:1270483/92‑1 (MQ=255)
                                                                                                                                                                 tgtTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGc                                                                                                               <  1:325445/204‑1 (MQ=255)
                                                                                                                                                                  gtTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTTTTGATGCTCACCTGCATTTTTAATTCATCCTGATTGAGCCGGTCCCTTTTTAAACCGCTACCCTGCCCCGCGGCCTGACGACCGTTATTTTCGACCCTCACGAAATTGTTAATGTGAATGGGGAAGCCGGGTTCGCCTGGGTTTGCCCCCGTTCCGCACaggc                                                                         >  2:527770/1‑239 (MQ=255)
                                                                                                                                                                  gtTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGACACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACagg                                                                          <  1:1963706/238‑1 (MQ=255)
                                                                                                                                                                   ttGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggca                                                                        >  1:494129/1‑239 (MQ=255)
                                                                                                                                                                   ttGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGGGCCCGGGGTTATTGATGACCTCCTGCATATTGAATCCAGCAAGAAGACCACGGTCCCTTTTGAAAACGCGACCCTGCCGCGCGGCCCGAAGACCGTTATTTGCGAACCCCATGAAATCGTCCACGTGATGGGCGAAGCAGGGTTCTCCTGGTTTTCCCCCTGTGCCCAACaggca                                                                        >  2:429674/1‑239 (MQ=255)
                                                                                                                                                                       cgcAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGATGCTCACCCGCATATTTAATCCAGCATGATGACGCCGGTCACTTTTTAAACCCCTACCCTGCCGCGCGGCCCGACGACCGTTATTTGCGAc                                                                                                                                                >  2:1814088/1‑166 (MQ=255)
                                                                                                                                                                       cgcAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTTATGCTCCCCTGCATATTGAATCCAGCATGATGACGCCGGTCCCTTTTGAAACCGCTACCCCGCCGCCCGGGCTGACGACCGTTATTTGCGACCCTCCTGAAATCGTCCACGGGATTGGCGAAGCCGGGTTCGGCCGGTTTGCCCGCCTGGGCCAAAcggcaaggc                                                                    >  2:2305439/1‑239 (MQ=255)
                                                                                                                                                                       cgcAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGCGCCAGGGGTTATTTATGCTCCCCTGCATATTGAAAACAGCATGATTAAGCCGGTTACTTTTGTAACCGCCACCCCTACGCGCGGCCCGACGACCGCTATTTTCGACACCCTTAGAAACCTCAAAGCGATAGGGGAAGCCGGGATCTCCCGGGTTTCGCGCGGTGCCCGACagcgaaggc                                                                    >  2:2014266/1‑239 (MQ=255)
                                                                                                                                                                           gAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGGCGGGGGTTTTTTCTGCTCATCAGTATATTTaaag                                                                                                                                                                                                                      >  2:1270483/1‑90 (MQ=255)
                                                                                                                                                                            aaTACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTc                                                                                                                                                                                                   <  1:931326/110‑1 (MQ=255)
                                                                                                                                                                            aaTACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACaggcaaggc                                                                    <  1:773495/234‑1 (MQ=255)
                                                                                                                                                                              tACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTg                                                                                                                                                    <  1:848879/155‑1 (MQ=255)
                                                                                                                                                                                    gATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAAc                                                               <  1:1736843/231‑1 (MQ=255)
                                                                                                                                                                                    gATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGt                                                           <  1:911211/235‑1 (MQ=255)
                                                                                                                                                                                     aTGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGTGACCCTCATGAAATCGGCCACGTGATGGGCGAAGCCGGATTTGCCCGGGTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTAcct                                                       >  2:100084/1‑236 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTTCCAGGGGTTTTTGATGCTCCCCTGCATATTTTATCCAGCAAGAAGACGCCGGTCACTTTTGaaa                                                                                                                                                                                         >  2:931326/1‑110 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTAt                                                                                                                                                                                                                                              <  2:525729/57‑1 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTAt                                                                                                                                                                                                                                              >  1:525729/1‑57 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCa                                                                                                                                                                                                                                     >  2:591579/1‑66 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCa                                                                                                                                                                                                                                     <  1:591579/66‑1 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAg                                                                                                                                                                                                                   <  2:296995/84‑1 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGaa                                                                                                                                       <  2:539993/160‑1 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGaa                                                                                                                                       >  1:539993/1‑160 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCAAATTGAATTCAg                                                                                                                                                                                                                   >  1:296995/1‑84 (MQ=255)
                                                                                                                                                                                      tGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTGATGCTCACCCGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCCGCCGCGCGGCCTGACGCCCGTTATTTGCGACCCTCATGAAAACGTCAACGGGATGGGCGAAGCCGGGTTCGCCCGGTTTGTCCGCCTGGCCGAACAGGCAAGGCAAAAACAGTAc                                                         >  2:773495/1‑235 (MQ=255)
                                                                                                                                                                                        cTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCCGGGGTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTTTAACCGCTACCCTGCCTCGCGGCCTGCCGGCCGTTAATTGCGACCCTCCTTAAATCGTCCACGTGATGGGCGGAGCCGGCTTCTCCTGGTTTTCCCCCTGTGCCGAACAGGCAAGGCCAAACCAGTACTTAc                                                     >  2:1736843/1‑237 (MQ=255)
                                                                                                                                                                                        cTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATg                                                                                                                                                                                                                                         >  1:538695/1‑60 (MQ=255)
                                                                                                                                                                                        cTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATg                                                                                                                                                                                                                                         <  2:538695/60‑1 (MQ=255)
                                                                                                                                                                                        cTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGa                                                                                                                                                                                           >  1:2024175/1‑106 (MQ=255)
                                                                                                                                                                                        cTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGa                                                                                                                                                                                           <  2:2024175/106‑1 (MQ=255)
                                                                                                                                                                                        cTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGc                                                                                       <  2:43511/203‑1 (MQ=255)
                                                                                                                                                                                        cTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGc                                                                                       >  1:43511/1‑203 (MQ=255)
                                                                                                                                                                                        cTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGGTTTTTGGTGCTCACCTGCAGATTGAATCCAGCATGATGACGCCCGCCCCCTTTGAAACCGCTCCCCTGCCGCGCCGCCTGCCGACCGTTTTTTGGGACCCTCAt                                                                                                                                          >  2:848879/1‑155 (MQ=255)
                                                                                                                                                                                          ccGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCg                                                                                                                                                  >  2:294592/1‑145 (MQ=255)
                                                                                                                                                                                          ccGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCg                                                                                                                                                  <  1:294592/145‑1 (MQ=255)
                                                                                                                                                                                          ccGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGt                                                 <  1:1384499/239‑1 (MQ=255)
                                                                                                                                                                                           cGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTTGTGCTCACCTGCATATTGAATCCAGCATGATGAAGCCGGTCACTTTTTAAACCGCCACCCTGCCGCGCGGCCCGACGACCCGTATTTGCGACCCTCATGGAATCGTCAACGTGATTGGCGAAGCCGGATTTGCCCGGTTTGTCCGCTGTGCCCAAACGGGAAGGCCAAACCAGTACTTACAGGTc                                                >  2:1916620/1‑239 (MQ=255)
                                                                                                                                                                                             gCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTTATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCCGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCa                                               >  1:507342/1‑238 (MQ=255)
                                                                                                                                                                                              cTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAg                                              <  1:2358090/238‑1 (MQ=255)
                                                                                                                                                                                                   cAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTTGTGCCCACCCGCATATTTAATCCAGCAAGATGACCCCGGTCACTTTTTAAACCCCTACCCTGCCGCGCGGCCTGACGCCCGCTTTTTGTGACCCTCATGAAATCGCGCACGTGAGGGGCGGAGCCGGGTTGCCCCGGTTTGCCCGCCGTGCCGCACAGGCCAGGCAAAACCAGAACTTAAAggcaactcctt                                         >  2:910612/1‑229 (MQ=255)
                                                                                                                                                                                                        gATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGcc                                                                                                                                                                                                                                       >  1:1622916/1‑45 (MQ=255)
                                                                                                                                                                                                        gATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCt                                                                                                                                                                                                                                       <  2:1622916/46‑1 (MQ=255)
                                                                                                                                                                                                         aTTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTAc                                   <  1:819276/238‑1 (MQ=255)
                                                                                                                                                                                                           tGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCg                                <  2:494129/239‑1 (MQ=255)
                                                                                                                                                                                                                     ggcgCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCg                                                                                                                                                  <  2:2128259/118‑1 (MQ=255)
                                                                                                                                                                                                                     ggcgCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCg                                                                                                                                                  >  1:2128259/1‑118 (MQ=255)
                                                                                                                                                                                                                        gcAACGGCGGTGCCAGGGTTTATTGATGCTCCCCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCt                    <  1:5323/238‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTaa          <  1:1825702/238‑1 (MQ=255)
                                                                                                                                                                                                                                         gTTTATTGATGCTCCCCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGc    <  1:429674/237‑1 (MQ=255)
                                                                                                                                                                                                                                         gTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGAACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGGGCCa  >  1:1023994/1‑239 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
AAATTTCATCACATTAGCCGGGCTGAATACCAGGAATTGTTAGCCGTTTCCCGTGGCGACGCTGTTGCCGATTATATTATTGATAATGTCTCTATTCTCGACCTGATCAATGGCGGAGAAATTTCCGGCCCAATTGTGATTAAAGGACGTTACATTGCCGGTGTTGGCGCAGAATACACTGATGCTCCGGCTTTGCAGCGGATTGATGCTCGCGGCGCAACGGCGGTGCCAGGGTTTATTGATGCTCACCTGCATATTGAATCCAGCATGATGACGCCGGTCACTTTTGAAACCGCTACCCTGCCGCGCGGCCTGACGACCGTTATTTGCGACCCTCATGAAATCGTCAACGTGATGGGCGAAGCCGGATTCGCCTGGTTTGCCCGCTGTGCCGAACAGGCAAGGCAAAACCAGTACTTACAGGTCAGCTCTTGCGTACCCGCCCTGGAAGGCTGCGATGTTAACGGTGCCA  >  NC_000913/3843985‑3844456

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: