Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	39.5%	63.1 / 12.8	38	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (1/22);  new base G (15/0);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.03e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCG  >  NC_000913/941478‑941954                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      ggcACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGa                                                                                                                                                                                                                                       <  2:296592/251‑1 (MQ=255)            ccGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGcc                                                                                                                                                                                                                            <  1:183546/251‑1 (MQ=255)                       cgccgGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGTCAGAAATCTAAt                                                                                                                                                                                                                 <  1:400362/251‑1 (MQ=255)                                  aCCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGat                                                                                                                                                                                                     <  2:70615/251‑1 (MQ=255)                                  aCCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGa                                                                                                                                                                                                      <  2:219845/250‑1 (MQ=255)                                                      gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCtat                                                                                                                                                                                <  1:48891/250‑1 (MQ=255)                                                                                            actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                                                                                                                         <  2:225991/251‑1 (MQ=255)                                                                                               cctccGCGCAAATTGCGGAAGGTGTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGGTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                                                                                                                      <  2:207163/251‑1 (MQ=255)                                                                                                        aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGG‑GCTTACTTTTTTTAAAAAGCACGCCAGAAAACTAAAA‑CCCGCCTGATC‑AACATACATCCCC‑GCTATACCAACACCGGGGCCCGGCGCGAAGAagaggggggcaccctttgtccgggg                                                                                                                             >  2:410850/1‑230 (MQ=255)                                                                                                         aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGG‑GCTTACTATCTTTAACAAGCACGCCCGAAAACTAAAG‑CCCCCACGTTC‑ATCATCGATCCGC‑GCCCTACCCACACCCGGGCCGGGGGCGCAGAAGAGGGGAACCCCATTCGTcccggg                                                                                                                             >  1:219193/1‑246 (MQ=255)                                                                                                             gCGGAAGGTTTGAACTATACCGACGGCGGCTGGGCAGATGGCAACAGCTCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                                                                                                                         <  2:341683/250‑1 (MQ=255)                                                                                                                   ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGG‑ACTTACTATCTTTAACAGGCACGCCCGAAAACTAATG‑CCCGCAAGATC‑ATCATCGATCCCC‑GCTATAACCACACCGGTGCCGGGCGCGGAGAAGAGGGGGTCCCCATTCGTCCGGGGACAGaatgcgc                                                                                                                  >  1:538512/1‑245 (MQ=255)                                                                                                                             ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa                                                                                                          <  2:412829/249‑1 (MQ=255)                                                                                                                             ataCCTACGGCGGCTGGGCAGATGGCAACAGCACGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa                                                                                                          <  2:22710/249‑1 (MQ=255)                                                                                                                                     ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑GCTTAATATCTTGAACAAGCACCCCAAAAAAATAAAT‑CCCGCATGATCCACCA‑CGATTCGCACCTTTCCCACCACCG‑CGCCGGGCGCGAAGATGAGTGGATTCCGATTTATCCCGGAACAGATGCCGCACTTGTTAAAGGTcctgg                                                                                                >  1:173255/1‑247 (MQ=255)                                                                                                                                      gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGTGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg                                                                                               <  2:413729/251‑1 (MQ=255)                                                                                                                                        ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑ACTTACTATCTTTAAAAGGCACGCCAGAAAATTCATG‑CCCCCATGATC‑CTCATCGATCCGC‑GCCATACCGACACCCGGGCCGGGCGCGGAGAAGGGTGGGTCCCTCTTCGGCCGGGCACAGATGCCGCACAGGTTTACgggccgtcgg                                                                                              >  2:203005/1‑242 (MQ=255)                                                                                                                                                aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc                                                                                                <  1:51323/240‑1 (MQ=255)                                                                                                                                                aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGG‑ACTTTCTATTTTGAACAGGCACCCCAGAAAACTAATG‑CCCGCCTGATT‑ATCATCGATCCGC‑GCTCTACCGACACCCGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCTTCCGGGAACAGATGCCGCCCCGGGTAACGGTCTGGc                                                                                                >  2:51323/1‑240 (MQ=255)                                                                                                                                                    ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑ACTTACTATCTTGAACAGGCACGCCAGAAAACTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGCCGCACCGGGTAAAGGGCTGGCGGAAGGCATGATCAc                                                                                 >  1:69504/1‑251 (MQ=255)                                                                                                                                                      cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCCTTA‑TATATTTAACAGGGCCGCCAGAAATCTAATGGCCCCCA‑GATC‑TTCATCGGTCCGC‑GCCATACCGACACCCGGGCCGGGGGCGGAGATGAGGGGATCCCTATTTGTCCCGGGACAGATGCCCCACCGGGTAAAGGGCTGGCGTACggactgattactg                                                                               >  2:288635/1‑239 (MQ=255)                                                                                                                                                             ccgacGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGTATGAGTGGCGGTGGTGTG‑ACTTTCTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGTATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCt                                                                        <  2:478169/247‑1 (MQ=255)                                                                                                                                                                  ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑ACTTACTCTCTTGAACAGGCACGCCAGAAATCTAAAG‑CCCCCCTGATC‑ATCATCCATCCGC‑GCTATAACGTCACCGGGGCCGGGGGCGAAGATGAGTGGGACCCTATTCGTCCGGGAAACGAAGCCGCACTGGGTAACGGGCTGGCGTACGTCATGATACCTTAAAAACtggtgg                                                                   >  2:371671/1‑251 (MQ=255)                                                                                                                                                                    agatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACGCCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGAt                                                                 <  2:219193/249‑1 (MQ=255)                                                                                                                                                                          aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCa                                                           <  1:176959/251‑1 (MQ=255)                                                                                                                                                                            aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑ACTTACTCTCTTTAACAGGCACGCCCGGAAACTAAAG‑GCCCCCCGATT‑ATCATCCCACCCA‑GCTATAACGACCCCGGGGCCGGGCGGGAAAAAGAGGGGATCCCTTTTCTGCCGGGAACAGATACCCCACTGGGTAAAGGTCCGGCGTAAGGCATTATCACTGAAAACCTGGTGGGTCAGGccat                                                         >  1:528938/1‑248 (MQ=255)                                                                                                                                                                              cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑ACTTACTATCTTGAAAAGGCACGGCAGAAAACTAATG‑CCCGCATTATC‑ATCATCGATCCGC‑GCCATACCGACCCCGGGGCCGGGCGCGAAGAAGAGTGGAACCCTATTCGTCCGGGGACCGATGCCGCACTGGTTAAAGGTCTGGCGGACGTCATGATCACTGAAAACCTGGTGGGTAAGGCATTcc                                                       >  1:72760/1‑251 (MQ=255)                                                                                                                                                                                gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑ACTTACTATCTTGAACAAGAACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGGAGATGAGGGGGTCCCCATTCGTCCGGGGACAGAAGCCGCACTGGGTTACGGGCCGGCGTACGGCATGCTCACCGAAAAACCGGTGGGACAGGCAATCCTc                                                     >  2:188586/1‑251 (MQ=255)                                                                                                                                                                                  aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGa                                                   <  1:311563/251‑1 (MQ=255)                                                                                                                                                                                        gTAGTGCTGATTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAata                                             <  2:66351/251‑1 (MQ=255)                                                                                                                                                                                              cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc                                                                                                <  1:403334/194‑1 (MQ=255)                                                                                                                                                                                                gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAATTATTGCGtt                                      <  2:163065/250‑1 (MQ=255)                                                                                                                                                                                                       aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑ACTTACTAACTTTAACAGGCACGCCCGAAATCTAATA‑CCCGCAAGATC‑ATCATCGATCCGC‑GCTATAACGACACCGGTGCCGGGCGCGAAGAAGAGGGGATCCCCATTCGTCCCGGGACAGATGCCGCACTGGGTAACGGTCTGGCGGACGTCATGATCACTGAAAACCTGGGGGGTCAGGCAATTCTCGATAAATATTGCGTTGGCTACg                               >  1:263223/1‑250 (MQ=255)                                                                                                                                                                                                          aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑AACTACTATCCTGAAAAAGCAAGCCCCAAAAATAATG‑CCCGCATGATT‑ATCATCGATCCGC‑GCCATAACCCCAACCGTGCCGGGCGCGGAGATGAGAGGAGTCCTCTTTGTCCGGGGACAAAAACCCCCCTAGTTAACGGTCTGGCGGCCGTCCTGATCAATGAAAAACCGGGGGGTCAGGCAATACTCGATAAAAATTGCGTTGGATAccaca                            >  2:332187/1‑246 (MQ=255)                                                                                                                                                                                                           aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGG‑GCTTACTATCTTGAACAGGCCCGCCCGAAATCTAAAG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTAATCGTCCGGGGACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCCTGATCACTGAAAAACCGGTGGGTCAGGCATTCCTCGATAAATATTGTGTTGGGTACCATg                            >  2:393423/1‑249 (MQ=255)                                                                                                                                                                                                                    aaaCGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTg                 <  2:27925/251‑1 (MQ=255)                                                                                                                                                                                                                          gAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTgccagcc           <  2:69504/251‑1 (MQ=255)                                                                                                                                                                                                                                   gCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCg  >  2:63597/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      GGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTG‑ACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATC‑ATCATCGATCCGC‑GCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCG  >  NC_000913/941478‑941954

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