Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	71,566	0	C	T	34.8%	78.3 / 19.0	46	L134L (CTG→CTA)	thiB	thiamine ABC transporter substrate binding subunit
Reads supporting (aligned to +/- strand):  ref base C (4/26);  new base T (16/0);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 4.89e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGTT  >  NZ_CP009273/71318‑71810                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         gtggtgTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCa                                                                                                                                                                                                                                                         <  1:429564/250‑1 (MQ=255) ggggtgTAACTCAGTACAAGTCTGCTTTCCCCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCTTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCa                                                                                                                                                                                                                                                         <  2:384753/248‑1 (MQ=255)    gtgtAACTAAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGtt                                                                                                                                                                                                                                                      <  1:125066/250‑1 (MQ=255)                          tttCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGCTTTTTAGTTTGTTTTTTTTTTTAACG‑AAggg                                                                                                                                                                                                                               >  1:242441/1‑249 (MQ=255)                           ttCCCCTTTTAAAAACAGGCCGTAGCCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCCTCAGCGCCATAGC‑CTTTTGGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGc                                                                                                                                                                                                                               <  2:680959/249‑1 (MQ=255)                            tCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCg                                                                                                                                                                                                                              <  1:366115/249‑1 (MQ=255)                               ccTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCATCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTTTTTTTGTTATAAATG‑AAAGCGGAGt                                                                                                                                                                                                                          >  1:616522/1‑250 (MQ=255)                                   ttAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTTTTTGTCATAAACG‑AAggaggacaata                                                                                                                                                                                                                       >  1:611915/1‑240 (MQ=255)                                     aaaaaCAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTATAAAAG‑AAGGGGGAGAAGCCa                                                                                                                                                                                                                     >  1:567512/1‑249 (MQ=255)                                             gccgagggCTTCGCCCCAGCTTTTGGTGACGGTGCCCGTTTTCTTCCCCAGTTTCTGCCAGGC‑TTGGGGGGCGTCTCCGCCATAGC‑CTTTTTGATTCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGACCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAc                                                                                                                                                                                                            <  2:207077/244‑1 (MQ=255)                                                  aGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTGTTCTTAACG‑AAAGGGAAGAAGGCATAAACaaaacggac                                                                                                                                                                                                       >  1:423559/1‑244 (MQ=255)                                                     cTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCTGTTTTTTTTTTTTCTTAAAG‑GAGAAGAAAGAGCCATAACAAAAC‑GGtgagg                                                                                                                                                                                                    >  1:589332/1‑245 (MQ=255)                                                      ttCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTTTTTTTTTCTTTAAAA‑AAAGGGGAGGAGCCCTAAACAAAC‑AGTACGa                                                                                                                                                                                                    >  1:210043/1‑248 (MQ=255)                                                            ccAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTGTTTTAAACG‑AAGGCGAAGTAGGCATAATCAAAC‑CGTACGaaaagggct                                                                                                                                                                                            >  1:490175/1‑244 (MQ=255)                                                               gCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTCGTTTGTTTTTTTCATAAAAG‑AAGGCGAAGTAGGCCTAATTAAAA‑GGGACGAAAGTATCatta                                                                                                                                                                                         >  1:215087/1‑250 (MQ=255)                                                                 cTTTGGTGACCGTGCCCGTTTTTTTCGCCGGTTACTGCCAGGC‑TTGTGGGGCGTCTTCGCCATAGA‑CTGTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCattatt                                                                                                                                                                                       <  1:561790/250‑1 (MQ=255)                                                                 cTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTGTTTTTTTTCTAAAAG‑AAGGGGGAGTAAGCATAATCAAAA‑GGTACGAAAGGGTTattttt                                                                                                                                                                                       >  1:328759/1‑250 (MQ=255)                                                                        gaccgtgaccgtTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTCTTTTTATAAACG‑AAGGGGAAGTAGCCCTAATTAAAA‑CGGACGCAAGAGGCATTATTTCAaccg                                                                                                                                                                                >  1:202934/1‑250 (MQ=255)                                                                        accgttgccctttttCTTCCCCATTTTCTCCAAGGC‑TTGGGGGGCTTCATCGCCTTGGC‑CTTTTTGCACCCATACCAACAGCCCCAGCCCCGTTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAgccg                                                                                                                                                                                <  2:190513/239‑1 (MQ=255)                                                                         acccgtgtctgtttCCTCCCCAGTTTCCGGCAAGCATTGGGGGGCG‑CAACCCCCTAGATCTTTTTGCATCAAAGCCAACAGACCCACCCCCGGGGCAGGGGTCCGCGGTCCCTAAAAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAgccg                                                                                                                                                                                <  2:522986/239‑1 (MQ=255)                                                                                ccgtTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGGTTGTTCTTTTCATTAACA‑AAGGCGAAGTAGTCCGCATCAACCAGGTACAAAACTGTC‑TTATTCCTGCCGCCGGGGAc                                                                                                                                                                        >  2:377044/1‑250 (MQ=255)                                                                                          cgccgttttCCGCCAGGC‑TTGGGGGGGTTCATCGCCATAGA‑CTTTTTGCACCCATAGCAACAGCCCCACCCCCGGTGTCCTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACg                                                                                                                                                                 <  2:648089/241‑1 (MQ=255)                                                                                          cGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTc                                                                                                                                                                                                                                                          <  1:249453/160‑1 (MQ=255)                                                                                           cccAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGC‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCAt                                                                                                                                                             <  1:626258/249‑1 (MQ=255)                                                                                             cAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTGTCATAAAAA‑AAAGGGGAGGAGGCATAAACAAAA‑GGGAAGAAAGAGGCATTATTTCAGCCGCCCGGAACGCTTACGGcctac                                                                                                                                                           >  1:487869/1‑246 (MQ=255)                                                                                              gtttacTGCCAGGC‑TGGGGGGGCTCCATCGCATTAGC‑CTTTTTGCACCCATAGCAACAGACCCAGCCCGGGGGTAGTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCg                                                                                                                                                          <  2:83468/245‑1 (MQ=255)                                                                                                               gtgGGGCGTCATCGCAATGGA‑CTTTTTGCACCCATACCAACAGCCCCACCCCCGGTGTACTGGTGCGCGGACCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCttt                                                                                                                                          <  1:538727/249‑1 (MQ=255)                                                                                                                     ggTCTTCGCCATAGC‑TTTTTGGCTTCCAAAGCACCAGACCCAGCCCCGGTGTCCGGGTGCGCGGATCCTGATAAATCTCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCaaa                                                                                                                                   <  1:345543/249‑1 (MQ=255)                                                                                                                       tCATCGCCATAGT‑CTTTTTGCACCCATAGCACCAGCCCCACCCCCGTTGTACTGGTGCGCGGACCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACa                                                                                                                                 <  2:48121/250‑1 (MQ=255)                                                                                                                                  gA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTGTTTTTGTTATAAAAGAAAGGGGAAG‑AGCCATAAACAAAC‑GGTACGGAAATGTCAATATTTCCGCCCCCGGGGACGGTAACGGCATCCGCTG‑CCCCACCCCCTTTTTCAAACAGTCCGGttttt                                                                                                                      >  1:674465/1‑249 (MQ=255)                                                                                                                                               ccATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCt                                                                                                         <  1:457421/250‑1 (MQ=255)                                                                                                                                                    gCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTAATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGGTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAACC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCTAACAGTCCGGTTTTACTGGCGG‑CGTCTAACa                                                                                                     <  1:523296/249‑1 (MQ=255)                                                                                                                                                        cAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTTTTTTTTATAAAAG‑AAGGGGGAGGAGCCATAAATAAAC‑GGGACGAAAAGGTCATTATTTCCGCCGCCCGGGACGTTTACGGCATCCGCCG‑CCCCACCCCCTTTTGCAAACAAACCCGGTTTTCTTTCGG‑CGTCTCACAGggtg                                                                                                >  1:370712/1‑250 (MQ=255)                                                                                                                                                              cggccgtggggtACGGGGGCGCGGTCCCTTAAAATTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTGCTTGCAGGCGTTGTGGCGGGTTTTTCAGTTTGTTCTGGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATcc                                                                                          <  2:642037/240‑1 (MQ=255)                                                                                                                                                                  cccGGTGTACTGGTGCGCGGATCCTGATAATTCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGcc                                                                                      <  1:367471/250‑1 (MQ=255)                                                                                                                                                                   ccGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGccc                                                                                     <  1:252402/250‑1 (MQ=255)                                                                                                                                                                       tgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTGTTTTTTTCATAAACG‑AAGGCGAAGGAGCCATAATCAAAC‑GGTACGAAAGGGTCATTATTCCAGCCGCCGGGGACGCTTACGGCATCCGCCG‑CCACACCGCCTTTTGCAAACACACCGGTGTTACTTGCGG‑CGGCTAACAGGGTGTTATCCAGCCCCacc                                                                                 >  1:254198/1‑248 (MQ=255)                                                                                                                                                                          aCTGGTGCGCTGATCCTGATAAAGCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGcacc                                                                              <  2:246132/250‑1 (MQ=255)                                                                                                                                                                                                   aCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCt                                                      <  1:447516/249‑1 (MQ=255)                                                                                                                                                                                                              ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTGGTTTTTTCTTTTCATTAAAG‑AAAGCGAAGAAGGCATAATCAAAA‑GGGTCTAAAGTGGCATTATTACAGCCGCCGGGAAACGTAACAGCATCCCCTGCCCCCACCCCTTTT‑GCAAAAAGACCGGGTTTTCTTGAGGCCGGCTACAAGGTTGTTTTACA‑CCCCCGCACACCACCAGGCCTACTGTTTTTTTCTTcctttccctgc                                          >  2:673667/1‑241 (MQ=255)                                                                                                                                                                                                                 aTCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAATGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACg                                        <  2:650645/249‑1 (MQ=255)                                                                                                                                                                                                                              aGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGaaa                          <  2:47811/250‑1 (MQ=255)                                                                                                                                                                                                                              aGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGaaa                          <  1:529864/250‑1 (MQ=255)                                                                                                                                                                                                                                          cTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCa              <  2:567512/250‑1 (MQ=255)                                                                                                                                                                                                                                                   gggTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGccac       >  1:525792/1‑248 (MQ=255)                                                                                                                                                                                                                                                      tttttCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGtt  >  2:37275/1‑250 (MQ=255)                                                                                                                                                                                                                                                      tttttCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGtt  >  1:296432/1‑250 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                         GTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGC‑TTGTGGGGCGTCATCGCCATAGA‑CTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACG‑AAGGCGAAGTAGCCATAATCAAAC‑GGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTG‑CCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGG‑CGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTGAGAAGCGAAACGCCATCTTCCAGCGCCACCAGTT  >  NZ_CP009273/71318‑71810

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