Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,133	0	C	G	23.8%	101.5 / 10.6	42	intergenic (+116/‑82)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (27/5);  new base G (0/10);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 2.04e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACT  >  NC_000913/3780890‑3781365                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            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)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            TGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACT  >  NC_000913/3780890‑3781365

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