Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,120	0	T	C	25.6%	82.9 / 12.8	45	intergenic (+103/‑95)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (30/2);  new base C (0/11);  total (30/15)
Fisher's exact test for biased strand distribution p-value = 1.36e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACT  >  NC_000913/3780877‑3781365                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                            aTCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGTCGCTGTATGTACT‑CCCTCACTCTT‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGc                                                                                                                                                                                                                                                     >  2:679061/1‑251 (MQ=255) aTCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCCCCCCCAGGGAGCGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AACCG‑CCACCGCACTCTTTTTTTCTCCCT‑CGc                                                                                                                                                                                                                                                     >  1:4098/1‑251 (MQ=255)   cGAAAAAGAGATGAAAGTGGCAATGACGCTGAGTGGCGCGGAATCGATCAGCGAAATTACGCAAGATGCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGc‑cc                                                                                                                                                                                                                                                  <  1:679020/251‑1 (MQ=255)     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)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                            ATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACT  >  NC_000913/3780877‑3781365

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