Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,152	0	A	G	33.3%	70.3 / 13.3	43	intergenic (+135/‑63)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (26/2);  new base G (0/14);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 2.27e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACT  >  NC_000913/3780927‑3781401                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            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)                                                                                                                                                                                                                   cccgg‑gggggTGGGCCGGGGGGGGGG‑GAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTAcc                        <  2:653249/247‑1 (MQ=255)                                                                                                                                                                                                                         gggAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCCCGAGTTTACCGCCGtt                  >  1:54224/1‑251 (MQ=255)                                                                                                                                                                                                                                         gAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCCTCCTGACt  >  1:124662/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            CAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGC‑GGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACT‑CCCTCACTCTG‑AAACG‑ACACCGCACTCTTTTTTTCTCCCT‑CGC‑CCCTC‑CG‑GGGAGAGGGCCGGGGTGAGGG‑G‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACT  >  NC_000913/3780927‑3781401

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