Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,131	0	T	C	34.2%	62.4 / 10.9	42	intergenic (+114/‑84)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (24/3);  new base C (3/11);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 3.11e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACT  >  NC_000913/3780881‑3781365                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCt                                                                                                                                                                                                                                                >  2:495578/1‑251 (MQ=255)          tGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGgaga                                                                                                                                                                                                                                     >  1:248141/1‑251 (MQ=255)           gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGagag                                                                                                                                                                                                                                    >  2:227035/1‑251 (MQ=255)                                               cAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCa                                                                                                                                                                                              >  1:800847/1‑251 (MQ=255)                                                     aaTTACGCAAGCTTCGCTGGTGCAGGGGCTGGGTCAAGAGTTGCATGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTCTGTACT‑CCCCCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCCC‑CG‑GGGCGAGGGCCGGGGAGAGGG‑G‑ACAAGCCCGCTCACCTCACCAATTTTc                                                                                                                                                                                        >  2:303955/1‑251 (MQ=255)                                                      nttcgtaaagtttggtgtgTGGAGGGGGTGGG‑AAAGGGTTCGCTGCGGGCCGTGGCCCCCGTGGCGAAGGGGACGGCGGCATGGCGGTTGGCCCCCCCTCCCCATCCCCCCCCCCCTCGGGGGGGGGGACCCTTTGGGCGCTG‑GTGTTCC‑CCCCCCCCCGGAAAACG‑CCCCCCCACTTTTTTTTCCTCCCC‑CCC‑CCCTC‑GG‑GGGGGGGGGCCGGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                                                                                                                                                       <  1:495550/233‑1 (MQ=255)                                                                atcgggtGGTGAGGGGGCTGGGAAAAGGGGTCGCTGG‑GGACCTGCCCCCCAGGGAGAAAGGGAGGGCGGCAAGTTCGTTTGCCCCCCCCCCCCCTACCCCTCTCCCGAGGGGGGGGGGGCCTTTTGGGCGCTTTTTTCCCC‑CCCTCACTCTG‑AAACG‑CCCCCGC‑CTTTTTTTTTCCCCCC‑CCC‑CCCCCGGG‑GGGGGGGGGCGGGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                                                                                                                                             <  1:452515/245‑1 (MQ=255)                                                                                      tAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                                                                                                                       >  2:654307/1‑251 (MQ=255)                                                                                                            ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa                                                                                                                                  >  2:744090/1‑251 (MQ=255)                                                                                                              cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCCC‑CG‑GGGAGAGGGCCGGGGTGCGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc                                                                                                                                >  1:2840/1‑251 (MQ=255)                                                                                                                          aaGGGAATGCGGCAAGTTCGTTGGCGCCCCCCCCCCCAACCCCTCCCCGACGGGGAGAGGGGACCGTGCGGCGTGTTTTGTCCC‑CCCGCACTCGG‑AAACC‑ACCCCGCCCTTTTTTTTCCCCCCT‑CCC‑CCCCC‑GG‑GGGAG‑GGGCCGGGGGGGGGG‑GAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGc                                                                                                                    <  1:744043/251‑1 (MQ=255)                                                                                                                                             gTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAACCATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcataat                                                                                                 >  1:378512/1‑251 (MQ=255)                                                                                                                                                                  cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCCCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                                                                                                                                       >  2:33258/1‑176 (MQ=255)                                                                                                                                                                                    ggggACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATatcat                                                          >  1:304031/1‑251 (MQ=255)                                                                                                                                                                                    ggggACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCCCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATCCCGCAAAATACTGGCAACATCATCCGTCTCTGCGCTAATACCGGCTTTCGTATGAATatcct                                                          >  2:172717/1‑249 (MQ=255)                                                                                                                                                                                    ggggACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTGTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCTGGGTGATGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTCTGCTACACACCGTAATTTACGAACCAGAAATTTCGCCACATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCGGCATatcag                                                          >  2:737263/1‑250 (MQ=255)                                                                                                                                                                                     ttgtcacaagagcagtgtatgttCTCCCCCCACCTTG‑ACACG‑AC‑CCCCACTCTTTTTCTCCCCCC‑CCC‑CCGCG‑GG‑GGGCGGGGGGCGGGGTGGGGG‑GAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATatcat                                                          <  2:303480/236‑1 (MQ=255)                                                                                                                                                                                       ttgtgaaatgcgagagttatccccccctccctCCT‑AAACGCACCCCCCACT‑TTTTTTTCCCCCC‑CCC‑CCCCC‑CG‑GGGGAGGGGCCGGGGGGGGGGAA‑AAAAGGCCGAACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATatcatc                                                         <  2:74785/227‑1 (MQ=255)                                                                                                                                                                                            ttCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACGGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGa                                                  >  2:551397/1‑251 (MQ=255)                                                                                                                                                                                                     gtatgtaCT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATatcat                                                          >  1:58271/1‑234 (MQ=255)                                                                                                                                                                                                          taCT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCCC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTgg                                    >  1:372917/1‑251 (MQ=255)                                                                                                                                                                                                                 ccccccacctacagccccacCCCTCTTTTTTTTGCCCCCC‑CCC‑GCCCC‑GG‑GGGAGGGGCGCGGGGGGGAGG‑G‑CAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACg                                <  2:688080/234‑1 (MQ=255)                                                                                                                                                                                                                    ccctcg‑AAACG‑CCCCCCCCCTCTTTTTCCCCCCCC‑CCC‑CCCCC‑GG‑GGGGGGGGGCGGGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAg                           <  1:419390/246‑1 (MQ=255)                                                                                                                                                                                                                    aCTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTgggcagaaaag                           >  2:419410/1‑243 (MQ=255)                                                                                                                                                                                                                      ccccggaaccgcacccGCCCTTTTTTTGTCGCCCC‑CCC‑CCCCC‑CG‑GGGAGGGTGCCGGGGGGGGGG‑A‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAg                           <  1:202657/237‑1 (MQ=255)                                                                                                                                                                                                                      acttggacgagcccccGCTCTTTTTTTTTCCCCCC‑CGC‑CCCCC‑GG‑GGGAGAGGGGGGGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAg                           <  2:372937/237‑1 (MQ=255)                                                                                                                                                                                                                      accctacacgcaccccgcatcttttttctCCCCCC‑CGC‑CCCCC‑GG‑GGGGGGGGGCGGGGGTGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAg                           <  1:654267/229‑1 (MQ=255)                                                                                                                                                                                                                       cctcgcaacagcccccccacttgtttttCCCCCTGTGC‑CCCCC‑GG‑GGGAGGGCGCGGGGGGGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAg                           <  1:551363/228‑1 (MQ=255)                                                                                                                                                                                                                       ccccaaaaacaccaccgtctttttttttCCCCCT‑GCC‑CCCAT‑CG‑GGGAGGAGGCCGGGTTGGGGG‑GAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAg                           <  2:608968/230‑1 (MQ=255)                                                                                                                                                                                                                        ccactgaagagccaccgctctTTTTTTTACCCC‑CGCGCCCTC‑CGTGGGGGAGGGCGGGGTAGAGGT‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAg                           <  1:172707/234‑1 (MQ=255)                                                                                                                                                                                                                               g‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GTGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCCCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCCCTAATACCGGCTTTCGTCTGCATATCATCGCACCCATGGGCTTTGCCTGGGACGATAAGCCCCTGcccc                 >  1:620821/1‑251 (MQ=255)                                                                                                                                                                                                                               g‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGCACCAGAAATTCCGCCAAATACTGGCAACATCATcc                                                                                               >  2:780157/1‑173 (MQ=255)                                                                                                                                                                                                                               g‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGa                                  >  1:303463/1‑234 (MQ=255)                                                                                                                                                                                                                               g‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGAc                                 >  1:74780/1‑235 (MQ=255)                                                                                                                                                                                                                               g‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGcgcc                 >  1:562064/1‑251 (MQ=255)                                                                                                                                                                                                                               g‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGcgccg                >  1:688038/1‑251 (MQ=255)                                                                                                                                                                                                                                      gacctctTTTTTCCCCCCC‑CCC‑CCCCC‑GG‑GGGGGGGGGCGGGGTTGGGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCgg           <  1:562263/248‑1 (MQ=255)                                                                                                                                                                                                                                         ctctTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATatcat                                                          >  1:327542/1‑201 (MQ=255)                                                                                                                                                                                                                                         ctctTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGc        >  2:202669/1‑251 (MQ=255)                                                                                                                                                                                                                                               ttttCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGGTGGGTGTTGCCTGGGACGATAATCGCCTGCGCCGCGCGGGCTTGGACt  >  1:117247/1‑251 (MQ=255)                                                                                                                                                                                                                                               ttttCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGAtt                                           >  1:307148/1‑210 (MQ=255)                                                                                                                                                                                                                                               ttttCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACAGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACt  >  2:562298/1‑251 (MQ=255)                                                                                                                                                                                                                                                               ccTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGcgcc                 <  2:5382/222‑1 (MQ=255)                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               AAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACT  >  NC_000913/3780881‑3781365

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