Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,856,028	0	A	G	29.6%	73.1 / 12.1	54	pseudogene (936/1617 nt)	BW25113_RS19145	PTS sugar transporter subunit IIB
Reads supporting (aligned to +/- strand):  ref base A (15/23);  new base G (16/0);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.50e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTGC‑GGG‑AA‑TCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCG  >  NZ_CP009273/3855783‑3856277                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                   aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTGC‑GGG‑AA‑TCa                                                                                                                                                                                                                                                               <  2:748476/250‑1 (MQ=255)   cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTGC‑GGG‑AA‑TAAG‑c                                                                                                                                                                                                                                                            <  1:1149692/250‑1 (MQ=255)      aaTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GGATT‑AG‑GGTGGG‑GGG‑AA‑tagacaa                                                                                                                                                                                                                                                           >  1:1003707/1‑244 (MQ=255)       aTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTGC‑GGG‑AA‑TCAG‑CAAGc                                                                                                                                                                                                                                                        <  2:757836/250‑1 (MQ=255)       aTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGCGAAATGAACAGGAAGG‑AAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GC‑AA‑GGGTGC‑GGG‑GA‑CCAG‑Cagtc                                                                                                                                                                                                                                                        >  2:328197/1‑247 (MQ=255)           cGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTGC‑GGG‑AA‑TCAG‑CAAGCCCG‑c                                                                                                                                                                                                                                                   <  1:1369249/250‑1 (MQ=255)            gATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TTAAA‑GATTGC‑‑GG‑AA‑ACAG‑CAAAACCG‑cc                                                                                                                                                                                                                                                  >  1:1542029/1‑250 (MQ=255)            gATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GT‑AA‑GGTTGC‑GGG‑GA‑ACAG‑CCAACCCG‑cc                                                                                                                                                                                                                                                  >  1:304378/1‑250 (MQ=255)                   gTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTC‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑AGGTGG‑GGG‑AA‑ATAG‑CCAGCCCC‑CCCATTTTa                                                                                                                                                                                                                                           >  1:560014/1‑250 (MQ=255)                            atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GT‑AA‑AGTTGCGGGG‑AG‑TCAC‑AAACCCCC‑CCCCTTTTccccctttt                                                                                                                                                                                                                                   >  2:298647/1‑245 (MQ=255)                                ccGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGCGGG‑TC‑AA‑GGTTGG‑GGG‑AT‑ACCG‑CAAGCCCC‑CCCCTTTTccccagcttcccg                                                                                                                                                                                                                               >  2:363218/1‑241 (MQ=255)                                  gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGCT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GT‑AA‑GGGTGC‑GGG‑AA‑ACAA‑CAAACACC‑CCACCCTTTCCCC‑A‑TTTTTCgggg                                                                                                                                                                                                                          >  1:902427/1‑248 (MQ=255)                                        ttGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑TGGTGG‑GGG‑AA‑ACAG‑AAAGCCCG‑CCACTTTTACCCGCA‑TTTTCAgggcggg                                                                                                                                                                                                                       >  1:575094/1‑243 (MQ=255)                                        ttGCTGAACATCGGGATCCAGTTTCGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGCGGG‑TT‑AA‑GGTGGG‑GGG‑AA‑ACGG‑AAAACCCC‑CCCCTTTTCCCCC‑A‑TTTT‑CGGTGCCGcgg                                                                                                                                                                                                                    >  1:1008887/1‑248 (MQ=255)                                           ctaaaCATCGGGCTCCAGTTTGGCGGTAAAC‑CCTAGTCAATCACCCCTCCCCCCGTGTCGCCCA‑CCACACAAAAGAGATACATTCCGGTGAACATG‑GGGGCCGCCAGCAGGGGGTGTACCGCAAACAGCAAGGGTAAAATGAACAGAAAGG‑TAAATTCCAGCGGTT‑CGGTAATCCCCC‑CCAGCAGGGC‑GG‑TT‑AG‑GGTTGG‑GGG‑AA‑AAAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGa                                                                                                                                                                                                                 <  2:730958/247‑1 (MQ=255)                                           cTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GT‑AA‑AGTTGG‑GGG‑AACACAC‑ACAGCCCG‑CC‑CTTTTACCCA‑ATTTTtccggtccggagg                                                                                                                                                                                                                    >  1:1121382/1‑245 (MQ=255)                                              taaactgggtaccAGTTTTGCGGAAAAACCCTAGTAAAACGGCCCCCCCCCTATTTCCCCCC‑CCACCCAAAAGGAAAACATTGCGG‑CGAAATTGGGGGCCGCCAGCGGGGTGTGCACCGAAAACAGCAGGGGGAAAAGAAACAGGAAGG‑T‑AATTCCCGGGGTTGCGTTAATCCCCA‑CCAGCAGGGG‑GG‑TT‑AG‑GGTGGC‑GGA‑AA‑TCAG‑CAGGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGa                                                                                                                                                                                                                 <  2:719718/237‑1 (MQ=255)                                              aaCATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑AGATGT‑GGG‑AA‑TCAA‑CCAACCCC‑CCCCCTTTATCCC‑A‑TTTTCTCGTGCCGCCGTGCA‑g                                                                                                                                                                                                              >  1:773775/1‑249 (MQ=255)                                                cATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GT‑AA‑GGGTGC‑GGG‑AA‑TAAG‑CAAACCCC‑CCCCCTTTACCCG‑A‑ATTTTCGGTGCGGCAGGGAA‑GTa                                                                                                                                                                                                            >  1:420394/1‑249 (MQ=255)                                                 aTCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑AGGTTC‑GGG‑GA‑AAAA‑CAAACCCA‑CCACTTTTTCCCC‑A‑TTTTCCGGGGCGGCCGGGaaagaca                                                                                                                                                                                                           >  1:835698/1‑245 (MQ=255)                                                      gaccctgTTTTGCGGAAAAA‑CCAAGCCACTCACCCCCCCCCCCATGTGCCCCC‑CCACACAAAAGAGATACATTACGGCCGAAGTT‑GAGGCCGCCGGCACGGGGTACACCACAAACCGCAGCGGGAAAATGAACAGAAAGG‑TAAATCCCAGGGTTG‑CGATAATCCCCC‑CCAGCATGGC‑GT‑TT‑AA‑GTTTGG‑GGA‑AA‑TGAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCt                                                                                                                                                                                                     <  2:566137/244‑1 (MQ=255)                                                       aTCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑TA‑GGGTGC‑GGG‑GA‑GCAG‑CAAGCCCG‑CCCCTTTTACCCC‑A‑TTTTTCAGTTCGGGCGGGAA‑GGACATCGcc                                                                                                                                                                                                     >  1:317965/1‑248 (MQ=255)                                                          cAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGGTGG‑GGG‑GA‑TCAG‑CAAACCCG‑CCCCCTTTTCCCC‑A‑TTTTCCGGTGCGGCGGGGGA‑GGACATCGCTcaa                                                                                                                                                                                                  >  1:1133040/1‑249 (MQ=255)                                                               ttGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑TA‑GGGTGC‑GGG‑GA‑TCAG‑CAAACCCG‑CCACTTTTTCCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACAACGCTAaaaaaaat                                                                                                                                                                                            >  1:783202/1‑248 (MQ=255)                                                                     tAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GT‑AA‑GGTTGC‑GGG‑GA‑TCAG‑CAAACCCG‑CCCCCTTTTCCCC‑A‑TTTTTCCGTGCGGCCGTGAA‑GGACATCGCTAAAAAAATG‑CCCACg                                                                                                                                                                                     >  1:1468708/1‑250 (MQ=255)                                                                             tAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GT‑AATGGTTG‑‑GGG‑GA‑TTAA‑CAAACCCG‑CCACTTTTATCCG‑A‑GTTTTCCGTGGGGGCGTGAA‑ATAAACCGCTCAAAAAAAA‑CCCACGGCGccc                                                                                                                                                                               >  1:527971/1‑248 (MQ=255)                                                                             tAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑GAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑TA‑AGTTTC‑GGG‑AA‑TTAG‑CAAACCCG‑CCACCTTTACCCG‑A‑TTTTTCCGTGGGGCAGGAGA‑GGACATTACGTAAGAAATT‑CCCACCGCAccc                                                                                                                                                                               >  1:1475788/1‑248 (MQ=255)                                                                                          ccccccccATGGTGCCCACCCACACAAAAGAGAACAATTGCGGGCGACAGG‑GAGGCCGCCAACAGGGGGGCTCCCAAAAACAGCAA‑GGGGAAAAGAAAGGGGAGG‑AAAATTTCCGGGGTTGCGGTAATCCCCA‑CCAGCGGGGC‑GG‑TA‑AG‑GGTGGG‑GGG‑AA‑TCAG‑AAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGt                                                                                                                                                                  <  2:363806/245‑1 (MQ=255)                                                                                           cccGCACATGTCGCCCA‑CCACACAAAAGAGATACATTACGGTCACCATT‑GGGGCCGCCAGCACGGCGTGTCCCGCAAACAGCAACGGTAAAATGAACAGGAAGG‑AAAATTCCAGCGTTT‑CGTTAATCCCCC‑CCAGCGTGGC‑GT‑TT‑AA‑GGTGGC‑GGG‑AA‑TCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTa                                                                                                                                                                <  1:813920/249‑1 (MQ=255)                                                                                              gCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑TA‑AGGTGC‑GGG‑GA‑ACAG‑CAAACCCC‑CCCCCTTTTCCCC‑G‑TTTTCCGGGGGGGCAGGGAA‑GGACATCCCCAAAAAAAAGCCCCCCGGCCCCAAAAATTTTTGGTTTAcg                                                                                                                                                              >  1:438956/1‑249 (MQ=255)                                                                                                ccATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑TA‑GGTTTC‑GGG‑AA‑ACAG‑CAAACCCG‑CCCCTTTTTCCCC‑A‑TTTTTCCGGGCGGGCGGGAA‑GTACATCGCCAAAAAAAAG‑ACCACGGGCCCCAACATTTTTTAGTTAACggg                                                                                                                                                           >  1:936270/1‑248 (MQ=255)                                                                                                              cacaccaAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GT‑AA‑AGGTGG‑GGG‑GA‑ACAG‑CAAAGCCC‑CCACCTTTTCCCG‑ATTTTTCCGGGCGGGCGGGGAG‑GAAAATGCCTAAAAAAATA‑CCACCGGCCCCAAAACTTTTTTGTTTACCG‑GCAGGGGaaaaacc                                                                                                                                              >  2:1055002/1‑247 (MQ=255)                                                                                                               accacaAAGATACACATGACGGTCACCATG‑GAGCCCCCCGCCACGCCGTGAACCGAAAACAGAAGCGGAAAAATGAACGGGAGGA‑TAATTCCCGGCGGTT‑CGGTAATCCCCA‑CCAGTATGGC‑GG‑TT‑AG‑GGTTGC‑GGG‑AA‑CCAG‑AAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc                                                                                                                                             <  2:449469/246‑1 (MQ=255)                                                                                                               acaccaAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑GT‑AA‑AGGTGG‑GGG‑GA‑TTAG‑CAAAACCC‑CCCCTTTTACCCC‑A‑TTTTTCGGGGGGGCGGGGGA‑GTACATCGCCTAAGAAAAA‑CCCACGGCCCCAAAAATTTTTTAGTTACCGTGCCGGGGCAAACAGc                                                                                                                                             >  1:844344/1‑250 (MQ=255)                                                                                                                       gagaTACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGGTGC‑GGG‑GA‑TCAG‑CAAACCCC‑CCACCTTTTCCCG‑C‑TTTTCCGGTGGGGCAGTGGA‑GGAAATCGCTAAAAAAATA‑CCCACGGCCCCAAAGATTTTTTAGTTACCGTGCAGGGCAAAACCGCCCTCC‑gg                                                                                                                                     >  1:720549/1‑249 (MQ=255)                                                                                                                        agaTACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑TA‑AGTTGC‑GGG‑GA‑TAAC‑CCAACCCGCCCCCTTTTCCCCA‑C‑TTTTCGGGGTGGGCGGGGAG‑GGACACCCCCAAAAAAACT‑CCCCCGCCCCCAAAAATTTTTTGTTTACCGTTCCGGGCCAACCCCCCTTCC‑ggg                                                                                                                                    >  2:1005639/1‑250 (MQ=255)                                                                                                                                   ggTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTGC‑GGG‑GA‑TCAG‑CAAGCCCG‑CCCCTTTTACCCG‑A‑TTTTTCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCCCGGCACCCAAGATTTTTTAGGTACCGTGCAGGGCAAAAACGCCCTTC‑GGGAAACACGAttt                                                                                                                         >  1:432450/1‑249 (MQ=255)                                                                                                                                      cGACATG‑GAGGCCGCCGCCGCGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGT‑TAAATTCCAGCGGTT‑CGGTAATCCCCC‑CCAGCATGGC‑GT‑TA‑AG‑GGTTGC‑GGA‑AA‑TCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGc                                                                                                                     <  1:1458904/250‑1 (MQ=255)                                                                                                                                        aaatagggggCCGCCGGCAGGGCGTGCGCCGAAAAAAGAAAGTGTGAAAGAAAGAGGAAGT‑AAAATCCCGGGGGCT‑GGGTAATCCCCA‑CCGGCTGGTG‑GG‑TG‑AG‑GGTGGG‑GGG‑AA‑ACAA‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCg                                                                                                                    <  2:839532/242‑1 (MQ=255)                                                                                                                                                   cgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑AGGTGC‑GGG‑GA‑ACAA‑CAAGGCCC‑CCCCTTTTTCCCG‑A‑ATTTTCGGTGCGGCAGGGAAGGAAAACCGCTAAAAAAAGC‑CCCACGCCCCCAACATTTTTTGAGTTACCGTGCGGGCAAAAACCCCCTCCC‑GGGAAAAAAGATTTCA‑CGGCCC‑gggcgcc                                                                                                        >  2:1155243/1‑245 (MQ=255)                                                                                                                                                       aGCGCGGTGTATCCCGAAAACAGAAACGGAAAAATAAAAAGGAAGG‑AAAGTCCCAGCGTTT‑CGGTAATCCCCC‑CCAGC‑TGGG‑GG‑TT‑AA‑GGGTGG‑GGGAAA‑TCAG‑AAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAAc                                                                                                    <  2:97617/250‑1 (MQ=255)                                                                                                                                                        gccgggCGTGAGCAGCAAACAG‑AGAGGTAAAAAGAACAAGAAGGTTAAACTCC‑GGGGTT‑TCGAAATTCCCACCCAGCA‑GGC‑CG‑TT‑GG‑GGGTGG‑GGG‑AA‑ACAGCCAGGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACt                                                                                                   <  2:1511234/246‑1 (MQ=255)                                                                                                                                                         cgcggcggGTCCCACAAACAGCACCGGTGAAATGAAAGGGAAGG‑TAAATTCCAGCGGTT‑GGGTAATCCCCA‑CCAGCAGGGC‑GG‑TT‑AA‑GGTGGC‑GGG‑AA‑TCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTg                                                                                                  <  1:1505840/248‑1 (MQ=255)                                                                                                                                                                    ccGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTGC‑GGG‑GA‑TCAG‑CAAACCCG‑CCCCTTTTACCCG‑G‑TTTTTCGGTGCGGCAGGGAA‑GGAAATCGCCTAAGAAAAG‑CCCACGGCACCCAAGATTTTTGAGTTACCGTGCAGGGCCAAACCGCCTTCCGGGGAACAA‑GGTTTCAGCGGCGCGGGCGC‑CAAAATGAACCCttctc                                                                                        >  2:1234164/1‑246 (MQ=255)                                                                                                                                                                                    taaaaTAAACAGGAAGT‑AAAATCCCGGCGTTT‑CGGTAATCCCCC‑CCAGCATGGC‑GG‑TT‑AA‑GGTGGC‑GGG‑AA‑TCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTAc                                                                       <  1:386461/248‑1 (MQ=255)                                                                                                                                                                                        aTGAAAAGGAAGT‑AAAGTCCCAGCGGTG‑CGGTAATCCCCA‑CCGGAAGGGG‑GG‑TT‑AG‑GGTGGG‑GGA‑AA‑CCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCt                                                                   <  1:783860/250‑1 (MQ=255)                                                                                                                                                                                                 aaGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTTC‑GGG‑AA‑TCAG‑CACGCCCC‑CCCCTTTTACCCG‑A‑TTTTCCGGATCGGCAGTGGA‑GAACATCGCTAAAAAAAAG‑CCCACCGCCCCCAAGATTTTTTTGGTACCGGGCCGGGCAAAACCGCCCTCC‑CGGGACAAACATTTCCGAGGCTC‑CGGGGCCCCACTGAGATCATGCCGACGCTGTGCCCCGGACCTCTGCATggccc                                                          >  2:1162906/1‑246 (MQ=255)                                                                                                                                                                                                        aactCCAGGGGTT‑CGGTAACCCCCA‑CCAGCATGGC‑GG‑TT‑AG‑GGTTGC‑GGA‑AA‑TCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCa                                                     <  1:1465432/246‑1 (MQ=255)                                                                                                                                                                                                           tCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTGC‑GGG‑AA‑TCAG‑CAAGCCCG‑CCCCTTTTTCCCG‑A‑TTTTTCGGTGCGGCCGTGAA‑GTACATCGCTAAAGAAAAG‑CCCACCGCACCCAAAATTTTTTAGTTACCGCGCCAGGCAAAAACGCCCTCC‑GGGGACCAAGAATTTAGCCGCCC‑CGCGGCCACACCGAACACCTTCACAAGCTGCCCCCCATACAACTGAATgcgccccgccaccc                                                 >  2:663799/1‑238 (MQ=255)                                                                                                                                                                                                            ccAGCGGTT‑CGTAAACCCCCC‑CCAGCATGGC‑GT‑TT‑AG‑GTTTGC‑GGA‑AA‑GCAT‑CACGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAg                                                 <  1:751388/249‑1 (MQ=255)                                                                                                                                                                                                                       ggTATTCCCCC‑CCAGGATGGC‑GG‑TG‑AG‑GGTTGC‑GGG‑AA‑CCCG‑CAGCCCCG‑CCCCTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                                      <  1:328197/250‑1 (MQ=255)                                                                                                                                                                                                                       ggTAATCCCCC‑CCAGAAGGGC‑GT‑TT‑AA‑GGTGGC‑GGA‑AA‑TCAG‑CAGCCCCG‑CCACTTTTACCCG‑A‑TTTTCCAGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAcc                                       <  1:655140/249‑1 (MQ=255)                                                                                                                                                                                                                       gacaaTCCCCA‑ACAGCGGGGC‑GG‑TA‑AG‑GGTGGG‑GGA‑AA‑TCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAcc                                       <  1:19829/246‑1 (MQ=255)                                                                                                                                                                                                                                            gggggttaaa‑GGTGGG‑GGA‑AA‑CCAC‑CAGGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTa                     <  1:322785/242‑1 (MQ=255)                                                                                                                                                                                                                                                                  aa‑TCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAATGGTGTAAACCg  >  2:809358/1‑250 (MQ=255)                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                   AACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTGCCCA‑CCACACCAAAGAGATACATTACGGTCGACATT‑GAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGG‑TAAATTCCAGCGGTT‑CGGTAATTCCCA‑CCAGCATGGC‑GG‑TT‑AA‑GGTTGC‑GGG‑AA‑TCAG‑CAAGCCCG‑CCACTTTTACCCG‑A‑TTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATG‑CCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCC‑GGGAACAACGATTTCAGCGGCTC‑GGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGGTGTAAACCG  >  NZ_CP009273/3855783‑3856277

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