Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0252682,233,2390TA34.6% 102.7 / 33.3 84intergenic (+139/+234)yohP/dusCsmall membrane protein YohP/tRNA dihydrouridine(16) synthase DusC
Reads supporting (aligned to +/- strand):  ref base T (14/37);  new base A (27/0);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 1.05e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTT  >  NZ_CP025268/2233111‑2233378
                                                                                                                                 |                                                                                                                                            
tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCGCCCCCAGCACCCCCggg                                                                                                                                   >  2:2469822/1‑140 (MQ=255)
    atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGGCCAAAGCCCCCCCCCCCAccccccccccgcgcg                                                                                                                               >  1:338125/1‑133 (MQ=255)
    atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCGCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCGCCGCCCCCCGCACCCCCCGGcgct                                                                                                                               >  2:177772/1‑140 (MQ=37)
         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGCGGCCAAAACCCCCCCCCCCCAGCACCCCCCGGCACAGGCGa                                                                                                                          >  1:4518164/1‑140 (MQ=25)
         gTCCGAAAAAAACGATAAAATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCTCCGCCGCCCCCTGCACCCCCGGGCTCAGGCGa                                                                                                                          >  2:1840921/1‑140 (MQ=21)
            cGAAAATAACGATAATATTTCAGTCTGTCATCATGTTGCGCTCCTGTCATGTGCATTTCTTCATATAATCACTGTCGCAAGGAGCGC‑GCGGGGGGGGGCCCACAGCCCCCGCCCCCCCCACCCCCGCGCCCCGGCGacac                                                                                                                       >  1:3250180/1‑137 (MQ=255)
            cGAAAATAACGAGAATATTTCCGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGGCCATCGCCGCCGCCCCCAGCACCCCCCGGCTCTGGCGAACa                                                                                                                       >  2:1137758/1‑140 (MQ=25)
            cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGCGCGC‑GCAGGGGGCGGCCAAACGCCGCCCCCCCCAGCACCCCCCGGCTATGGCGACCa                                                                                                                       >  2:3429528/1‑140 (MQ=25)
            cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAA‑AGCCCCCCCCCCCACCCCCCCCCCGCACAGGCGAACa                                                                                                                       >  1:5086741/1‑140 (MQ=255)
            cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAACGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACa                                                                                                                       >  2:5478160/1‑140 (MQ=255)
                    aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATACTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATCGCCGCCGCCCCCAGCACCCCCCGGCGCTGGCGAACAAAAAcgcc                                                                                                               >  2:533402/1‑140 (MQ=21)
                    aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGGGGGGCCCAAGCCCGCGCCCCCAGCACCCCCGCGCCCTGGGGGACaaaaaagcc                                                                                                               >  1:6233075/1‑135 (MQ=255)
                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGcgcg                                                                                                       <  2:6074919/140‑1 (MQ=255)
                              ttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGt                                                                                                     <  2:2540438/140‑1 (MQ=255)
                               tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTACATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                                    <  1:5944908/140‑1 (MQ=255)
                                  gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGccc                                                                                                 <  1:4596370/140‑1 (MQ=255)
                                      ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATCTACTCACTGGCGCTAGGAGCGC‑GCCGGGGGCGGCGAATCGACGCCGCCCCCAGCACCCCCGGGCTCGTGCCTACAAAATAGCCGCTGCGCGGGGGCCGCgg                                                                                             >  2:1446302/1‑140 (MQ=255)
                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATAGACTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCCGGCTCTGGCGAACAAAAACGCCGCAGCGCGGTGCCCTCGGc                                                                                            >  2:6191067/1‑140 (MQ=18)
                                          tCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTa                                                                                        <  2:482024/140‑1 (MQ=255)
                                          tCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTa                                                                                        <  2:2947400/140‑1 (MQ=255)
                                          tCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGA‑TTa                                                                                        <  2:5607223/140‑1 (MQ=255)
                                           ccTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTAt                                                                                       <  2:3033912/140‑1 (MQ=255)
                                               ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGGGCAAGGAGCGC‑GCCGGGGGCGGCCAATCGCCGCCGCCCCGTGCAGCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTTTCCCt                                                                                   <  1:2894524/140‑1 (MQ=18)
                                               ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAACGGCGCCGCCCCCAGCACCCCCCGGCTCTGGCGGACAAAAAAGCCGCTGCGCGGTGGCCTAGGC‑ttttacag                                                                                   >  2:3119261/1‑134 (MQ=18)
                                               ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCCCCCCCCCCAACACCCCCCGGCTCAGGCCGACAAAATCCCCGCTTCGCGGTGCCCACGGC‑TTATacca                                                                                   >  2:1788382/1‑139 (MQ=18)
                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCGCCGCCCCCTGCACCCCCGGGCTCAGGCGAACAAAATCGCCGCTGCGCGGTGCCCGCGGC‑GTATCCCATa                                                                                 >  2:2790123/1‑140 (MQ=18)
                                                    cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGGCAATCCCCGCCGCCCCCAGCACCCCCCGGCTCGGGCGAACAAAAAAGCCGCCGCGCGGTGCCCACGGC‑GTATaccagaccg                                                                              >  2:217282/1‑134 (MQ=18)
                                                    cTCCTGTCATGTGCATTGCGACATATAATCACTGGCGCTAGGAGCGC‑GCAGGGGGCGGCCAATCGACGCCGGCCCCATCAACCACGGGCACAGGCGCACAAAATCCCCGCAGCGCGGTGCCCTCGGC‑TTATCCCTTACgg                                                                              >  2:6029523/1‑140 (MQ=17)
                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCACGGC‑GTATACCATACGGcaa                                                                           >  2:5514631/1‑138 (MQ=17)
                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCGCCGCCCCCCGCCCCCCCCGGCAAAGGCGACCAAAAACGCCCCAGCGCGGGGCCCGCGGC‑GTATCCCATACGGcaa                                                                           >  2:3531872/1‑138 (MQ=17)
                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACGCCCCCCCCCCCAGCACCCCCCGGCTCAGGCGAACAAAATCGCCGCTGCGCGGGGCCCACGGC‑ATATCCCGTAAGGcga                                                                           >  2:484390/1‑138 (MQ=17)
                                                         gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTAcc                                                                         <  2:1223201/140‑1 (MQ=255)
                                                               tgCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCCGGCGCTGGCGCACAAAATCCCCGCTGCGCGGTGCCCACGGC‑GTATCCCGTAAGGCAAACGGGGCg                                                                   >  2:4391359/1‑140 (MQ=14)
                                                               tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCAAACCCCCCCCCCCCAGCACCCCCCGGCACAGGCGAACAAAATACCCCCTGCGCGGTGCCCACGGC‑AAATACCATACGGCGACCCCGTCg                                                                   >  2:598304/1‑140 (MQ=9)
                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGCCCAACCCCGCCCCCCCCAGCACCCCCCCGCACTTGCGGACAAAAAAGCCGCGGCGCGGGGGCCGCGGC‑GTATACCATAAGGCGAACCGGTCggg                                                                 >  2:3918679/1‑140 (MQ=255)
                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGGGGCCAAAAGCCGCCGCCCCCCGCACCCCCCGGCGCTGGCGAACAAAAAACCCCCGGCGCGGTGGCCACGGC‑GGATACCATAACGCTACCcgggaggg                                                                 >  2:2313847/1‑135 (MQ=9)
                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCggg                                                                 <  2:752402/140‑1 (MQ=255)
                                                                 cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCggg                                                                 <  1:3755834/140‑1 (MQ=255)
                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCgggag                                                               <  2:2637648/140‑3 (MQ=255)
                                                                    tGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGgcgc                                                              <  2:2649059/140‑1 (MQ=255)
                                                                    tGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGAACAAAAAAGCCGCTTCGCCGTGCCCAAGGC‑ATATACCATACCGCTACCCGGTAGGgcgc                                                              >  2:3546706/1‑140 (MQ=14)
                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCCATAGCCGCCCCCCCCAGCCCCCCCCGGCTCTGGCGAACAAAAAAGCCCCGGCGCGGGGGCCACGGC‑ATATACCATAACGCAACCCGGTCGGGCGCGa                                                            >  2:5236769/1‑140 (MQ=9)
                                                                           tataATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                                       <  2:6136056/140‑1 (MQ=37)
                                                                           tataATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                                       <  2:5703274/140‑1 (MQ=37)
                                                                           tataATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                                       <  2:1023299/140‑1 (MQ=37)
                                                                             taATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACa                                                     <  2:6321261/140‑1 (MQ=25)
                                                                              aaTCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAACATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCGGACGGCTACCGGGTCGGGCGCGAGGTAACAt                                                    <  2:1215651/140‑1 (MQ=255)
                                                                                tCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCCAACCCCCCCGCCCCCCGCACCCCCCGGCGCTGGCGAACAAAAAAGCCCCGGCGCCGGGGCCACCGC‑GAATACCATACGGCAACCCGGGCGGGCGCGCGGTAAccaac                                                  >  2:4986354/1‑136 (MQ=2)
                                                                                tCACTGGCGCAAGGAGCGC‑GCCGGGGGCGGCCCAAACCCGCCCCCCCCCGCACCCCCCGGCACAGGCGGACAAAATACCCCGGGCGCCGTGCCCTCGGCGTAAGCCCTTA‑GGCAACCCGGTCGGGCGCGAGGTATCATcc                                                  >  2:2635728/1‑140 (MQ=2)
                                                                                tCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAACCCCGCCCCCCCCAGCACCCCCCGGCTCAGGCGAACAAAAAAGCCCCAGCCCCGTGCCCCCCGC‑ATAACCCAAAAGGCAAACCGGTCCGGCGCCAGGTAACAAcc                                                  >  2:4465241/1‑140 (MQ=2)
                                                                                   cTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATAGCCGCCCCCCCCAGCCCCCCCCGGCTCTGGCGAACAAAAACGCCGCTGCGCGGTGCCCACGGC‑GTATACCATACGGCTACCCGGTCGGGCGCGAGGTAACaaaccag                                               >  2:1674609/1‑134 (MQ=11)
                                                                                   cTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAAAAGCCGCCGCCCCCAGCACCCCCGGGCTCTTGCGCACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCATAACGCTACCGGGTCGGGCGCGAGGTAACATcccag                                               >  2:1675074/1‑138 (MQ=11)
                                                                                      gcgcAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaa                                            <  2:5318038/140‑1 (MQ=32)
                                                                                        gcAAGGAGCGC‑GCAGGGGGCCGCCAAACGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGAACAAAAAAGCCCCGGCGCGGTGCCCACGGC‑TTATCCCTTACGGCTTCCGGGGGGGGCGCGAGGGAACATCCCAGTaaaa                                          >  2:769329/1‑140 (MQ=11)
                                                                                          aaGGAGCGC‑GCTGGGGGCGGCCCAACGCCGCCGCCCCCTGCACCCCCGGGCTAATGCGAACAAAATACCCGCAGGGCGGTGGCCTCGGC‑GTTTACCATACGGCAAACGGGTCGGGGGCGAGGGAAACACCCAGTaaaaag                                        >  2:705589/1‑138 (MQ=2)
                                                                                          aaGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                                        <  1:3405821/140‑1 (MQ=32)
                                                                                           aGGAGCGC‑GGAGGGGGGGGCCCATAGCCGCCGCCCCCTGCCCCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑GTATCCCTAACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                                       >  2:6009230/1‑140 (MQ=11)
                                                                                             gAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                                     <  2:5104303/140‑1 (MQ=32)
                                                                                                cgc‑gcTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                                  <  1:967010/140‑1 (MQ=17)
                                                                                                cgc‑gcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                                  <  2:1440438/140‑1 (MQ=32)
                                                                                                cgc‑gcAGGGGGCGGCCAATAGCCGCCGCCCCCAGCACCCCCGGGCACTTGCGAAAAAAATCGCCGGTGCGCGGTGCCCTCGGC‑TTATCCCATACGGCTACCGGGTAGGGCGCGAGGGTACAACCCTGTAAAAAGCGcccg                                  >  2:5539769/1‑139 (MQ=11)
                                                                                                cgc‑gcAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                                  <  1:40257/140‑1 (MQ=17)
                                                                                                 gc‑gcAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCATACGGCTACCGGGTCGGGCGCGAGGTACCATCCCTGTGAAACGCGCCCTc                                 <  1:733461/140‑1 (MQ=255)
                                                                                                     cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCGGTAAAACGCGCCCTCAGc                              <  2:3055755/140‑1 (MQ=17)
                                                                                                       gggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc                            <  2:1809798/140‑1 (MQ=32)
                                                                                                          ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca                         <  2:3309420/140‑1 (MQ=32)
                                                                                                          ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca                         <  2:1283432/140‑1 (MQ=32)
                                                                                                          ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca                         <  2:6393353/140‑1 (MQ=32)
                                                                                                          ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcaca                         <  2:735868/140‑1 (MQ=32)
                                                                                                             ggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATACCAGTAAAACGCGCCCTCAGCCCACAAcc                      >  2:5023845/1‑140 (MQ=12)
                                                                                                             ggcCAAACGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGAAAAAAATAGCCGCTGCGCGGTGCCCACGGC‑TTAAACCTTAAGGCTACCGGGTCGGGCGCGAGGTAACATACCAGTAAAAAGCGCCCACAGCCCACAAcc                      >  2:1883005/1‑140 (MQ=11)
                                                                                                              gcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa                     >  2:3059787/1‑140 (MQ=32)
                                                                                                              gcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑GTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa                     >  2:6439394/1‑140 (MQ=17)
                                                                                                              gcCAATCGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGAACAAAATAGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAACCa                     >  1:5005321/1‑140 (MQ=12)
                                                                                                                cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAACCAtg                   >  1:241834/1‑140 (MQ=17)
                                                                                                                cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCGTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCCGCCCACATCCCtg                   >  1:3496030/1‑140 (MQ=12)
                                                                                                                    cgccgccgcCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg               <  1:84771/140‑1 (MQ=32)
                                                                                                                    cgccgccgcCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTggg               >  2:1279/1‑140 (MQ=32)
                                                                                                                         ccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAACCATGTTGGCTGcc          >  1:1046095/1‑140 (MQ=17)
                                                                                                                         ccgccCCCTGCACCCCCGGGCTCTGGCGAAAAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcc          <  2:2040521/140‑1 (MQ=17)
                                                                                                                           gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc        <  2:2277808/140‑1 (MQ=32)
                                                                                                                             ccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCgg      >  1:5724534/1‑140 (MQ=32)
                                                                                                                              cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGc     <  2:1852861/140‑1 (MQ=32)
                                                                                                                              cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGc     <  1:1788382/140‑1 (MQ=32)
                                                                                                                                 tGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCtt  <  1:3389029/140‑1 (MQ=32)
                                                                                                                                 tGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCtt  <  2:2508419/140‑1 (MQ=32)
                                                                                                                                 |                                                                                                                                            
TGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGC‑TTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTT  >  NZ_CP025268/2233111‑2233378

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: ‑