Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,164,5250CA25.4% 144.9 / 12.5 71A144E (GCA→GAA) nagZbeta N‑acetyl‑glucosaminidase
Reads supporting (aligned to +/- strand):  ref base C (21/32);  new base A (0/18);  total (21/50)
Fisher's exact test for biased strand distribution p-value = 7.46e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGC  >  NC_000913/1164395‑1164654
                                                                                                                                  |                                                                                                                                 
ggttggttGATGGCCAGCGAAATGATCGCTTTGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCAGATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAg                                                                                                                          <  2:48575/140‑1 (MQ=255)
ggttggttGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAg                                                                                                                          <  2:217692/140‑1 (MQ=255)
ggttggttGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAg                                                                                                                          <  2:1711892/140‑1 (MQ=255)
        gATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTa                                                                                                                  <  1:886579/140‑1 (MQ=255)
         aTGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTAt                                                                                                                 >  1:402590/1‑140 (MQ=255)
           ggCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTg                                                                                                               <  2:743765/140‑1 (MQ=255)
           ggCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTg                                                                                                               <  2:743607/140‑1 (MQ=255)
           ggCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTg                                                                                                               <  2:52920/140‑1 (MQ=255)
           ggCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTg                                                                                                               >  2:394669/1‑140 (MQ=255)
            gCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGa                                                                                                              <  2:703102/140‑1 (MQ=255)
            gCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGa                                                                                                              <  1:64053/140‑1 (MQ=255)
              cAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATg                                                                                                            <  1:618406/140‑1 (MQ=255)
               aGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATgg                                                                                                           >  1:729466/1‑140 (MQ=255)
               aGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATgg                                                                                                           >  1:729241/1‑140 (MQ=255)
                gCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGt                                                                                                          <  1:1270592/140‑1 (MQ=255)
                     aTGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCa                                                                                                     <  1:1250197/140‑1 (MQ=255)
                               tGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCTGCCGGTTTATTGATGGTATGCATGAAGCCGGa                                                                                           <  1:881833/140‑1 (MQ=255)
                               tGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGa                                                                                           <  1:881924/140‑1 (MQ=255)
                                       gATATCAGCTTTGCGGCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATGCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGTAGCCGGAATGTAAAc                                                                                   >  2:763421/1‑140 (MQ=255)
                                       gATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCTGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAAc                                                                                   >  1:43842/1‑140 (MQ=255)
                                           tCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGAcc                                                                               >  1:390495/1‑140 (MQ=255)
                                                   gggcCTGTGCTGTTTGTCGGCCATATCCGCGCGGCGTTGGGCGAGGGTTATTATCTTGCCTAACCAAAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGTTATGCTTTAAGCCGGAATCAAAACGACCGGGAAaca                                                                       <  1:791957/138‑1 (MQ=255)
                                                   gcgcGTGGGCTGGTTGTGGGGCATATCAGCGGGGCGTTGGGCGAGGGTTCTTATCTTGCCGATCCACAAAACGCCCTGGCAATGGCCCGCGGGTTTATTGTTGGTATGCATGAAGCCGGATTGAACACGACCGGGAAaca                                                                       <  2:1219170/140‑1 (MQ=255)
                                                   gcgcCTGTGCTGGATGTCGGGAATAACAGCGCGGCGATTGGCGAGCGTTATTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGTTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAaca                                                                       <  2:285899/140‑1 (MQ=255)
                                                         gtgCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTccc                                                                 >  1:1203378/1‑140 (MQ=255)
                                                         gtgCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTccc                                                                 >  1:1307766/1‑140 (MQ=255)
                                                               tagttCGGGCAAAACAGCGCGGCTTTTGGCGTGCTTTCTTATAATGCCGATCCACAAAACGCCCTGGAAATTGCCAGCCGTTTTTTTTAGGTTATGAATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGaca                                                           <  2:1310439/136‑1 (MQ=255)
                                                               gttgcggggCATATCAGCGCGGCGATTGGCGAGGGTTCTTATCATGCCTTTCCACAAATAGCCATGGAATTTGCAAGCCGTTTTATTGTTGTTTTGATTTATGCCGTAATGAAAACGACCGGGAAACACTTCCCAGGaca                                                           <  2:796747/134‑1 (MQ=255)
                                                               gtggTCGGGCATATCCGCGGGGCGTTGGGCGCGGGTGATTATCTGGCCGACCCAAAAAACGCCGTGGCAATGGCCTGCGGTTTTTTTGCTGTCAGGCATGACGCGGCAATGAACACGACCGGGAAACACTTCCCAGGaca                                                           <  1:593208/137‑1 (MQ=255)
                                                               gATGTCGGGCATATCAGCGGGGCGATTGGGGCGCGTTCTTATCATGCCGAACCACAAAAAGCCCTGGCATTTGCCAGCCGGTTTATTGTTGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGaca                                                           <  2:398911/140‑1 (MQ=255)
                                                               gATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAATGCCCTGGCAATATCCAGCCGGTGTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGaca                                                           >  2:520221/1‑140 (MQ=255)
                                                               gATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGaca                                                           >  2:721247/1‑140 (MQ=255)
                                                               gATGTCGGGCATACCAGCGGGGCGTTGGGCGAGCGTTCTTATCATGCCTATCCACAAAACGCCCTGGCATTTGCCAGCGGTTTTTTTTTTGGTTTGAATGATGCCGTAATGAAAACGACCGGGAAACACTTCCCAGGaca                                                           <  2:1292263/140‑1 (MQ=255)
                                                               gATGTCGGGAATATCCGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAACCCCTGGCAATTGCCAGCCGTTTTATTGTTGGTATGCATGAAGCCGTAATGAAAACGACCGGGAAACACTTCCCAGGaca                                                           <  2:510585/140‑1 (MQ=255)
                                                               gATGTCGGGAATACCCGCGCGGCGTTTGGCTTGCGTTCTTATCATGCCGATCCAAAAAAACCCCTGTCATTTGCCCGCGGTTTTATTTATGGTATGAATTACGCCGTAATTAACACGACCGGGAAACACTTCCCAGGaca                                                           <  1:925247/140‑1 (MQ=255)
                                                                  gTCGGGCTTTTCAGCGCGGCGTTTGGGGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCATTTGCCTGCCGGTTTTTTGTTGGTTTGAATGAAGCGGGAATGAAAACGACCGGGAAACACTTCCCAGGACACgg                                                        <  2:502901/140‑1 (MQ=255)
                                                                    cGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAATGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTg                                                      >  1:1477293/1‑140 (MQ=255)
                                                                    cGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTg                                                      >  2:928657/1‑140 (MQ=255)
                                                                    cGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTg                                                      >  1:1479285/1‑140 (MQ=255)
                                                                      ggCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCa                                                    >  1:801724/1‑140 (MQ=255)
                                                                        caattCAGGGCGGCTTTTGGGGAGCTTTATTAACTTGCCGATCCACAAAAATCCCTGGCAATTGCATGCCGGTATATTGATGGATTGCATGAAGCGGTAATCAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGt                                                  <  2:1583453/136‑1 (MQ=255)
                                                                         atatCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTa                                                 >  2:1587412/1‑140 (MQ=255)
                                                                                   ggcttttGGCGAGCTTTCTTAACATGCCGATCAAAAAAACGCCCTGGAAATTGCCAGCGGTTTTATTGATGGTATGAATGAAGCCGGAATTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACt                                       <  1:600477/135‑1 (MQ=255)
                                                                                   ggcgTTTGGGGTGGGTTTTTATCTTGCCTATCAACAAAAAGCCCTGGAATTTGCCTGCGGGTTTATTGTTGGTTTGCTTGAAGCGGTAATTAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACt                                       <  1:1723610/140‑1 (MQ=255)
                                                                                   ggcgTTTGGCTTGCTTTCTTATCTTGCCTAACCAAAAAACGCCCTGGAATTGGCCAGCGGTTTTATTGTTGGTTTGCATGAAGCCGTAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACt                                       <  1:1635721/140‑1 (MQ=255)
                                                                                     gttttGGCGTGGTTTCTTATCATGCCTATCAACAAAATGCCGTGGAATTTGCCAGCGGTTTTTTTGTTGGTTTGCTTTACGCCGTAATTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTca                                     <  1:422606/137‑1 (MQ=255)
                                                                                     cgATTGGCGAGCGTTATTATAATGCCGATCAACAAAACGCCGTGGAATTTGCCAGCCGTTTTTTTGATGTTATGAATGAAGCCGGAATTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTca                                     <  1:1630477/140‑1 (MQ=255)
                                                                                       tttGGGGTGCTTTCTTTTCATGCGTATCCAAAAAAAGCCGTGGAAATTGCCAGCGGGTTTTTTGTTGTTTTGCATGATGCGGTAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  1:692034/139‑1 (MQ=255)
                                                                                       tttGGCGAGCTTTTTTAAAATGCCGATACAAAAAAAGCCCTGGAAATCGCATGCCTGTTTATTGATGGTTTCAATGAAGCCGTAATGAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  1:1683966/139‑1 (MQ=255)
                                                                                       tttGGCGAGCGTTCTTATCTTGCCGATCAAAAAAAAGCCCTGGAATTTGCCAGCGGTTTTTTTGTTGGTTTGATTTAAGCGGTAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  1:767305/139‑1 (MQ=255)
                                                                                       tttGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATGGCCAGCCGGTTTATTGATGGTATGCATTAAGCCGTAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  2:1630950/139‑1 (MQ=255)
                                                                                       tttGGCGAGCGTTATTATATTGCCTATCCACAAAAAGCCCTGGAAATTGCCAGCCGGTTTTTTGATGGTATGCATGAAGCCGTAATGAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  1:1671887/139‑1 (MQ=255)
                                                                                       ttgggCGTGCTTTTTTATCTGGCCGATCAAAAAAAAGCCTTGGCTATTGCCAGCCGGTTTTTTGTTGTTTTGCTTGAAGCGGCAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  1:767397/137‑1 (MQ=255)
                                                                                       ttgggCGTGCGTTCTTATCATGCCTACCCAAAAAATGCCCTGGAATTTGCCCGCGGTTTTTTTGTTGGTATGAATGACGCGGTAATCAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  1:394669/137‑1 (MQ=255)
                                                                                       tgtGGCGAGTTTTTTTATATTGCCTATCCACAAAATGCCGTGGAATTTGCCAGCGGTTTTTTTGTTGTTATGAATGACGCCGTAATTAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  1:71641/138‑1 (MQ=255)
                                                                                       aTTGGCGCGCTTTATAATAATGCCGAACCACAAAAAGCCCTGGAATTTGCCAGCCGGTTTATTTTTGGTATGAATGAAGCCGCAATGAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  1:497893/140‑1 (MQ=255)
                                                                                       aTTGGCGAGCGTTATTATATTGCCTATCAACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGTAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcaca                                   <  2:847576/140‑1 (MQ=255)
                                                                                       aTTGGCGAGCGTTATTAACATGCCGTTCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTTTTGATGGTATGAATGACGCCGTAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCATACTcaca                                   <  2:598435/140‑1 (MQ=255)
                                                                                        ttGGCGTGCGTTCTTATCATGCCGATCCAAAAAAAGCCCTGGCAATTGCCAGCCGGTTTTTTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcacac                                  <  1:888864/140‑1 (MQ=255)
                                                                                         tGGGGAGCGTTCTTATCTTGCCTAACCACAAAAAGCCCTGTAATTTGCATGCGGTTTTATTGATTTTTTCAATGACGCCGTAATCAACACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTcacaca                                 <  1:128412/140‑1 (MQ=255)
                                                                                           cctaGCTTTCTTAACATGCCTAACCAAAAAAAGCCCTGCAAATTGCCACCCGTTTTATTGATGGTATGCATGAACCCGCAATCAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACaaa                               <  2:353533/137‑1 (MQ=255)
                                                                                                            gCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATcc              >  2:692296/1‑140 (MQ=255)
                                                                                                             ccGATCCACAAATAGCCCTGGCAATTGCCAGCGGGTTTTTTGATGGTATGCATGAAGCCGGATTGATAACGACCGGGAAACACTTCCCAGGACACGGTGCTGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCa             <  2:868494/140‑1 (MQ=255)
                                                                                                             ccGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCa             >  2:1437277/1‑140 (MQ=255)
                                                                                                             ccGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCa             >  1:40115/1‑140 (MQ=255)
                                                                                                             ccGATCCAAAAAAACCCATGGCAATTGCCAGCCGTTTTATTGTGGGTATGCATTAAGCCGCAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCa             <  2:871655/140‑1 (MQ=255)
                                                                                                             ccGACCCACAAAAAGCCCTGGCAATTGCCAGCGGTTTTATTGATGGTATGCATTAAGCCGTAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCa             <  2:322187/140‑1 (MQ=255)
                                                                                                                aTCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGt          >  1:401134/1‑140 (MQ=255)
                                                                                                                    aaaaaacGCCCTGGATATTGCCAGCCGGTTTTTTTATGGTATGAATTAAGCCGTATTGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCcac      <  2:266902/138‑1 (MQ=255)
                                                                                                                    aaaaaaaGCCCTGGAAATTGCCAGCCGGTTTATTGATGGTATGAATGAAGCCGTAATCAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCcac      <  1:663953/138‑1 (MQ=255)
                                                                                                                        aaaGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGc  >  2:810947/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                 
GGTTGGTTGATGGCCAGCGAAATGATCGCTATGGATATTGATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGC  >  NC_000913/1164395‑1164654

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