Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,817,685	0	C	G	33.3%	113.2 / 19.1	60	intergenic (+143/‑87)	nadE/cho	ammonia‑dependent NAD(+) synthetase/excinuclease Cho
Reads supporting (aligned to +/- strand):  ref base C (32/6);  new base G (0/19);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 2.78e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGACACCATGCCGCTCTATATCGGCAA  >  NZ_CP009273/1817440‑1817929                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      cGTCCCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATTATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGgtg                                                                                                                                                                                                                                                   <  2:768781/250‑1 (MQ=255)    aCCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTTTaa                                                                                                                                                                                                                                                >  2:997966/1‑250 (MQ=255)         aaCAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑ttt                                                                                                                                                                                                                                          >  2:17338/1‑250 (MQ=255)           cAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑ttttt                                                                                                                                                                                                                                        >  2:238258/1‑248 (MQ=255)               tCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATc                                                                                                                                                                                                                                    <  1:905115/250‑1 (MQ=255)                             gagaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCtttt                                                                                                                                                                                                                      >  1:809171/1‑250 (MQ=255)                                aaCTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCtttttt                                                                                                                                                                                                                    >  1:630499/1‑249 (MQ=255)                                                                 ccGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAAc                                                                                                                                                                                  >  1:124560/1‑250 (MQ=255)                                                                           ccGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCg                                                                                                                                                                         >  1:734697/1‑249 (MQ=255)                                                                                 tCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCa                                                                                                                                                                   >  1:217622/1‑249 (MQ=255)                                                                                    ttgtgttatGGAAAAATTAATATTTTCCACATTTTGG‑TTGTGGAAAAAAAAAGTAATTTGTCCCCCCCGATTGCGGTGGGAACGCTTAATTCGCCCT‑AAAAAACTA‑GAAAAATCAAATACTTGGCGGGAG‑T‑GA‑AATAGG‑CCTGACGGGCTTAGCCCGTCGGAGGTTTTACCT‑TTTGCAATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAg                                                                                                                                                               <  2:253426/242‑1 (MQ=255)                                                                                      gatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTg                                                                                                                                                             >  1:603777/1‑250 (MQ=255)                                                                                       atTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTgg                                                                                                                                                            >  2:33177/1‑250 (MQ=255)                                                                                       atTTCTGAAAAAAGAAATAATTTGCACATAGTGGTTTGT‑GCAAAAAAGAGTAATTTTTCCACGCCGGAGGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TAAAAATACAATACATGGCAGGAG‑T‑CA‑AATAGG‑CCTGACAGGCGTAGCCCGTCAGAGGGTGAAACT‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTgg                                                                                                                                                            <  1:686448/250‑1 (MQ=255)                                                                                          tctaaaaaaaGAAATATTTTGCACATTTTGGATTGT‑GAAAAAAAGATTATTTTTTTCACGCCGGATGGGGCGTAAACCCCTTATTCGACCT‑AAAAAACTA‑TGAAAATTAAATACATTGCAGGAG‑T‑CA‑AATAGG‑CCTGAAGGGCGTACCACGTCAGAGGTTTTACCC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTa                                                                                                                                                          <  1:784978/244‑1 (MQ=255)                                                                                             gaaaaaaTTAATATTTTGCACATATGGGTTTTG‑GAAAAAAAGAGTAGTTTGTCCCCGCGGGAGGGGGCGTAACCCCTTTATCCGACCT‑ATAAAACTA‑TGAAAATCAAATACATTGAAGAGG‑T‑CG‑ATTGGG‑CGGGAAGGGCGTAGCCCGTCAGGGGTTGTACCT‑TTTGCCACCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGc                                                                                                                                                      <  2:882679/248‑1 (MQ=255)                                                                                                  aaGTAATAATTTGCAAATTTTGGATGTG‑GCAAAAAAGAGTATTTTTTTCCCGCCGGAGGCGGGGTGAACCCCTTATCCGACCT‑AAAAAACTA‑GGCAAATTAAATACATTGCAGGAG‑T‑CA‑AATAGG‑CGTGACAGGCGTAGCACGTCAGAGGGTGTACCT‑TTTGTCATCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTa                                                                                                                                                 <  1:50652/250‑1 (MQ=255)                                                                                                   aGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTaa                                                                                                                                                >  1:164197/1‑250 (MQ=255)                                                                                                          atTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTcc                                                                                                                                         >  1:865788/1‑250 (MQ=255)                                                                                                             tcccaatatTGGATTGG‑GCAAAAAGAATTAATTTTTCCACGCGGAATGCGGCTTGACCGCTTTATCCACCCA‑AAAAAACAA‑GCAAAATACAATACATGGCAGGAG‑T‑CA‑AATAGG‑CTGAACGGGTGTACCACGTCAGTGGGGTTACCC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCcgc                                                                                                                                       <  1:714161/244‑1 (MQ=255)                                                                                                               cacaTATTGGCTTGT‑GCAAAAAAGATTATTTTGTTCACGCGGGATGCGGCTTGAACCCTTTATTCGACAT‑AAAAAACTA‑TACAAATTAAATACATTGCAGGAG‑T‑CG‑AAGAGG‑CCTGCCAGGCGTAGCCCGTCAGAGGGTGTACCT‑TTTGTCATCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCgg                                                                                                                                     <  1:987637/249‑1 (MQ=255)                                                                                                                 caTATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTc                                                                                                                                  >  2:819363/1‑250 (MQ=255)                                                                                                                     ttGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATACAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGCGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGact                                                                                                                              >  1:1029071/1‑248 (MQ=255)                                                                                                                            tgtataaaaaaGTGTATTTTGTACCCGCCGGAGGCGGCGTACACTCTTTATCCGCACT‑AAAAAACTA‑TAAAAATAAAATACTTTGCGGGTG‑C‑CT‑AATAGG‑CCTGACGGGTGAACCGCTCCGGGTGGTCTATCT‑TTTTTCTCCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGaa                                                                                                                         <  2:832146/243‑1 (MQ=255)                                                                                                                            taagaaaaaaaGAATTATTTGCTCGCGCCGGATGCGGCGGAATCCTTTTATTAACCCT‑AAAAAACTA‑GACAAATACAACATATGGGGGGAG‑T‑AA‑AATGGC‑CATGAGGGGCGTAGCCCGTCAGAGGGTGTACCT‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGaa                                                                                                                         <  2:331144/243‑1 (MQ=255)                                                                                                                            gt‑gCAAAAAAGATTAATTTTTCCCCGCCGAGGGGGGGGTGACCGCCTTATTC‑ACCTAAAAAAACTA‑TAAAAATTAAATACATGGCGGGAG‑T‑CG‑AATAGG‑CTTGACGGGTGTAGCCGGTCAGCGGGTGTACCT‑TTTGTCACCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAg                                                                                                                        <  1:987442/249‑1 (MQ=255)                                                                                                                                  aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCgg                                                                                                                  >  2:874858/1‑250 (MQ=255)                                                                                                                                        agTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCCAATTTGAAGCTGCGGCAAttt                                                                                                            >  2:594274/1‑250 (MQ=255)                                                                                                                                         atattttttacACGCGGAATGGGGCATAAACGCTTTATTCACCCA‑AAAAAAAGA‑GCAAAATTAAATATATTGAAGGA‑‑T‑CGAAATGGG‑CCGGACGGGGGTAGCCCGTAAGGCGGTTTACCT‑TTTGTCATCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTa                                                                                                           <  2:138474/240‑1 (MQ=255)                                                                                                                                            aTTTTTTCCCGCGGGATGCGGCTTGACCGCCTTATTCGACCT‑AAAAAACAA‑TGAAAATTCAATACATGGCAGGAG‑T‑CG‑AATAGG‑CTTGACAGGCGTAGCACGTCAGAGGTTGTACCT‑TTTGTCACCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGa                                                                                                        <  1:40514/250‑1 (MQ=255)                                                                                                                                              ttGTACCCGCTGGATGTGCGGTAAACGCTTAATTCGACAT‑AAAAAACAA‑TGCAAATACAATACTTTGATGGGG‑A‑AG‑AATAGG‑CCGGCAAGGCGTAGCACGGCAGGGGGTGTACCT‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGaa                                                                                                       <  2:611391/249‑1 (MQ=255)                                                                                                                                                 ttCCCGCCGGATGCGGCGTGCACGCCTTATTCGCCCT‑AAAAAACAA‑TGCAAATTCAATACATGGCAGGGG‑T‑CA‑AATGGG‑CCTGACGGGCGTACCCCGTCAGAGGTTTTACCT‑TTTTTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATc                                                                                                   <  1:889138/250‑1 (MQ=255)                                                                                                                                                 ttCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGACTATc                                                                                                   >  1:600735/1‑250 (MQ=255)                                                                                                                                                         gggtcccgcgAGACCGCCTTATTCCCCCT‑ATAAAAATA‑ACCAATTTCAATACATTG‑ACGAGTT‑CA‑TATGGG‑CATGACGGGCGAAGCCCTTCGGAGGGTGTACCT‑TTTGCCACCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACAt                                                                                           <  2:419104/243‑1 (MQ=255)                                                                                                                                                         ggATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACAt                                                                                           >  2:422945/1‑250 (MQ=255)                                                                                                                                                           aTGCGGAGTAAACGCTTTATGCAACTA‑AAAAAACAA‑GGAAAAATAAATACA‑TGGAGGGG‑TGCG‑AATGGG‑CCGGACGGGCGTAGCCCTGCAGAGGGT‑TAACCTTTTGCCACCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACAttt                                                                                         <  2:409358/250‑1 (MQ=255)                                                                                                                                                                gcggaaaCGCCTAATTCCACCA‑AAAAACCTA‑AAAAAATTAAATCAAGTGCAGGAT‑T‑GA‑AATAGG‑CCTGACAGGCGTAGCACGTAAGGTGGTGTATCT‑TTTGTCACCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGtt                                                                                    <  2:539932/245‑1 (MQ=255)                                                                                                                                                                 cttaaccGCTTTATGCGACCT‑AAAAAACTA‑TGAAAATTAAATAAATTGCAGGAG‑C‑CA‑AATAGG‑CTTGACGGGGGTAGCCCGTCAGAGGGTGTACCT‑TTTTTCACCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTc                                                                                   <  1:589371/244‑1 (MQ=255)                                                                                                                                                                  gTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCa                                                                                  >  2:459790/1‑250 (MQ=255)                                                                                                                                                                          attATTC‑ACCA‑AAAAAAATATTGCAAATTCAAGACATTGCAGGGT‑C‑CG‑AATAGG‑CCTGAGGGGCTTAGCCCGTCAGACTGTAAACCT‑TTTGTCTTCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTaa                                                                          <  2:221782/249‑1 (MQ=255)                                                                                                                                                                            tATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAAt                                                                         >  2:331917/1‑249 (MQ=255)                                                                                                                                                                             atttccacaaaaaaaaCTA‑TGAAAATTAAACAGATTACGGGAG‑T‑C‑‑AAGAGGCCCTGACGGGCGTACCCCGTCAGATGTTCTCCCT‑TTTGCCATCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAAt                                                                         <  2:630499/237‑1 (MQ=255)                                                                                                                                                                              ttCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCATCTTTTTTAATAGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGACAACAGTGGTACGGCGTTTATCTTCTCCGCGGCTCGAATTTGAAGTTGAGTCAATTTTTGAATATCCCGAACATTTACGTTCATTCCTTatggac                                                                      >  2:1042779/1‑245 (MQ=255)                                                                                                                                                                                cGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACt                                                                     >  2:151588/1‑249 (MQ=255)                                                                                                                                                                                  aCCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTac                                                                  >  2:351343/1‑250 (MQ=255)                                                                                                                                                                                  aCCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGCCATCGACCCGCCTCTTTTTTACTCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTac                                                                  >  1:3124/1‑250 (MQ=255)                                                                                                                                                                                        taAAACTA‑TCCAAATTCAATACATGGCGGGAG‑T‑CG‑AATAGG‑CTTGACAGGCGTAGCCCGTAGGGGGGTGTACCT‑TTTGTCACCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTacccac                                                              <  1:615614/249‑1 (MQ=255)                                                                                                                                                                                         aaaaCTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTacccaccc                                                            >  1:675399/1‑250 (MQ=255)                                                                                                                                                                                              tA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGAc                                                         >  1:532153/1‑248 (MQ=255)                                                                                                                                                                                                        tCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTAt                                              >  2:625588/1‑250 (MQ=255)                                                                                                                                                                                                                    aGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGAGGGTGTACCT‑TTTGTCACCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGt                                  <  1:21690/250‑1 (MQ=255)                                                                                                                                                                                                                            cG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGt                            >  1:682802/1‑250 (MQ=255)                                                                                                                                                                                                                                aTAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGa                          >  2:74973/1‑249 (MQ=255)                                                                                                                                                                                                                                       cTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGACACCATg                   >  2:651511/1‑250 (MQ=255)                                                                                                                                                                                                                                       cTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGACACCATg                   >  2:109868/1‑250 (MQ=255)                                                                                                                                                                                                                                         gACGGGCGTAGCACGCAAGAGGTTGTACCT‑TTTTTCATCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGACACCATGcc                 <  1:238258/250‑1 (MQ=255)                                                                                                                                                                                                                                             ggCGTAGCACGTCAGACGTTGTAACT‑TTTGTCACCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGACACCATGCCGctc             <  1:908675/250‑1 (MQ=255)                                                                                                                                                                                                                                                   gCACGGCAGACGGTGTACCC‑TTTGTCATCGCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGACACCATGCCGCTCTATATc       <  1:4219/250‑1 (MQ=255)                                                                                                                                                                                                                                                    cccGTGAGAGGGTTTACCT‑TTTGCCATGCCCCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGACACCATGCCGCTCTATATCg      <  1:206695/248‑1 (MQ=255)                                                                                                                                                                                                                                                          aGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGACACCATGCCGCTCTACATCGGcca  >  2:920343/1‑246 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      CGTACCTCAACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGT‑GCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCT‑ATAAAACTA‑TGCAAATTCAATACATTGCAGGAG‑T‑CG‑AATAGG‑CCTGACAGGCGTAGCACGTCAGACGGTGTAACC‑TTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACTGGATAGATAACCAGCATTCGGAGTCAACAGTGGTACGGCGTTTAACTTCTCCGCGGCTCGAATTTGAAGCTGCGGCAATTTATGAATATCCCGAACATTTACGTTCATTCCTTAATGACTTACCCACCCGACCAGGGGTGTATCTGTTTCATGGTGAAAGTGACACCATGCCGCTCTATATCGGCAA  >  NZ_CP009273/1817440‑1817929

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