Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	326,935	0	C	T	59.7%	38.8 / 73.2	67	D105N (GAT→AAT)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (7/20);  new base T (20/20);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 7.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTT  >  NC_000913/326689‑327180                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gcGGACAAAACCCTCCCTGGTAATACCCGTTGATATCGGCCGGGCGGGGTTCGCCCGCCTGCGCGCCCGCTTAAATCATCGCTTCAAACCGCGGATTGCCGC‑GGGTTTGGGGGGTAGCGCCGCTCCCCGGGCCAGCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCgaga                                                                                                                                                                                                                                                    <  2:29358/250‑1 (MQ=255)  cGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTACGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTCGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCg                                                                                                                                                                                                                                                       >  2:268349/1‑247 (MQ=255)         aaaCCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGCGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGCTCCAgcgc                                                                                                                                                                                                                                            >  2:443245/1‑251 (MQ=255)                   gcggATACCCGTGGAAATCGTCGGGGCGCGGGTAGCCCCCCGGCACGCCCGTTTCAATCACCCCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TGGTGACGCTCACCGGGCAATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCAc                                                                                                                                                                                                                                    <  1:139720/246‑1 (MQ=255)                     tgATAACCGTTGAGATCGTACGTGCGCGGGTAGCCCGCCTGCTCGCCCGCTTCAATTATCGCTTCAAACCGGGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCCTCTCCGCCGTAATAGTCGTTTTCACCAATATCGCGAGTCTCGGCCTTGCGGAAGTAGAGCAGGCAGTCGAGGTCGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGTTCCAGCGCATTGCCACgggt                                                                                                                                                                                                                                >  2:380212/1‑249 (MQ=255)                       atacccGTTGAGCTCGTCCGTGCGCGGGTACCCCGCCTGCACGCCCGTTTCAATCATCCCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGCCGCCCACCGGCCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGt                                                                                                                                                                                                                              <  1:94814/247‑1 (MQ=255)                        tAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCACGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGt                                                                                                                                                                                                                              >  2:441611/1‑250 (MQ=255)                          cccGTTGAGACCGTCGGGGCGGGGGTACCCCGCCTGCACGCCCTCTTCAAACATCGCTTCAAACAGCGGATGGCGGCGGGGTTTGGAGT‑TAGTGACGCCCCCGGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGc                                                                                                                                                                                                                           <  2:477927/250‑1 (MQ=255)                           ccGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTCTTGAGATCCAGCGCATTGCCACGGATGTAGcc                                                                                                                                                                                                                          >  2:417730/1‑250 (MQ=255)                                gagctcgcccGTGCGCGGGTACCCCGCCTGCACGCCCGCTTCAACCATCGCTTCAAACAGCGGATGGACGCCGGGTTGGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGc                                                                                                                                                                                                                     <  1:417730/243‑1 (MQ=255)                                gagaccgcccggccgggggTAGCCCGCCTGCACGCCCGTTTCAATAATCGCTTAAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGCCGCTCCCGGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGc                                                                                                                                                                                                                     <  1:396085/235‑1 (MQ=255)                                       tCCGTGCGCGGATAGCCAGCCTGCCCGCCCGCTTCACTCATCGCTTCACACAGCGGCTTGACGCCGGGTTTGGCGG‑TAGTGACGCTCACCGCGCCATCACCGCCGTGATAGTCGTTTTCCCCCATATCGCGAGTCTCGGCCTTGCGGAAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGCCCTGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGCTGTAGCACATGCCGTTGAt                                                                                                                                                                                                              >  2:309483/1‑251 (MQ=255)                                          gTGCGCGGGTAGCCCGCCTGCACACACGCTTCCATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑AAGTGACGCTCACCGGGCCATCACAGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTCGTAGGGCAGGCAGTCGAGGAAGCTCCAGTTATCCAGACCGGGTTCTTGCGCCCAGTTATTGCGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATc                                                                                                                                                                                                             >  2:497807/1‑249 (MQ=255)                                           tGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGc                                                                                                                                                                                                          >  1:226818/1‑251 (MQ=255)                                              gggggtgccccGCCTGCGCGCCCGTTACACTCATCGCTACAAAACGCGGATTGACGCGGGTTTTGGAGG‑TTGTGCCTCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGcgac                                                                                                                                                                                                       <  1:136288/243‑1 (MQ=255)                                                   tAGCCCACCTGCACGCCCGCTTCAATCATCGCTTCACACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGCTAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGCAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCACAGTTATTGAGATCCAGCGCATTGCCACGGCTGTAGCCCATGCCGTTGATAAGCGACGATcc                                                                                                                                                                                                  >  2:376821/1‑251 (MQ=255)                                                              cACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCtt                                                                                                                                                                                       >  1:393056/1‑251 (MQ=255)                                                                       ttCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTc                                                                                                                                                                              >  1:427101/1‑251 (MQ=255)                                                                        tCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTcc                                                                                                                                                                             >  1:553523/1‑251 (MQ=255)                                                                                    tgaaccagcggattgcggggggttttGGAGG‑TAGCGACTCGCCCCGGGCAATCACCGCCGTGATAGTCGGTTTCACCCATAGCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGGCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTGCCGCGTCCGCACTCCATgc                                                                                                                                                                  <  2:122275/227‑1 (MQ=255)                                                                                    taaaaCAGCGGTTGGACGCGGGGTTGGGGGG‑TAGTGCCGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGAGCAGCGACGATCCACCCAGACCTTGACCGCGTCCGCACTCCATg                                                                                                                                                                   <  1:197917/247‑1 (MQ=255)                                                                                    tCAAACGGCGGATTGACGCCGGGTGTGGAGT‑TGGTGCCGCCCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATgc                                                                                                                                                                  <  1:343844/250‑1 (MQ=255)                                                                                    tCAAACAGCGGATTGAGGCGGGGTTGGGAGG‑TAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATgc                                                                                                                                                                  <  1:515178/250‑1 (MQ=255)                                                                                    ccaaacagtgtattgccccggggTTTGGGGG‑TAGTGACGCCCCCCGGGCCACCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACGCCATg                                                                                                                                                                   <  2:427101/229‑1 (MQ=255)                                                                                      aaaCAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATgcggc                                                                                                                                                               >  1:437877/1‑251 (MQ=255)                                                                                       aaCAGCGGATTGCCGCCGGGTTTGGAGG‑TAGTGCCGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCGCTGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGTATTGCCACGGATGTAGCACATGCCGTTGCTCAGCGACGCGCCACCCACACCTTTACCGAGCCCGCAACCCATgcggcg                                                                                                                                                              >  1:489936/1‑251 (MQ=255)                                                                                        acagcggtttgcccgggtgTTTGGAGG‑TGGGGCCTCTCACCGGCCCACCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGcggcgg                                                                                                                                                             <  1:139149/233‑1 (MQ=255)                                                                                         cAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCCTCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGGTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGcggcggg                                                                                                                                                            >  1:532928/1‑250 (MQ=255)                                                                                             ggATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTAtt                                                                                                                                                        >  1:22888/1‑251 (MQ=255)                                                                                                     gCCGTGTTGGGAGT‑TAGTGCCTCTCCCCGCGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACg                                                                                                                                                <  1:139058/251‑1 (MQ=255)                                                                                                           tttGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGGTATTCATAAAAGGTTc                                                                                                                                            >  1:265522/1‑249 (MQ=255)                                                                                                                aGG‑TAGTGACGCTCACCGGGCCATCACCGCCGCGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGAACGGGTTCTTGCGCCCAGTTATTGCGATCCAGAGCATTGCCACGGATGACGCCCATGCCGTTGCTCAGCGAAGCTCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTaatcctcaacagtttcaggtac                                                                                                                                     >  2:288167/1‑230 (MQ=255)                                                                                                                          cgttcccggggCCAACACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGGT‑CCGTTTCATAg                                                                                                                           <  1:443245/243‑1 (MQ=255)                                                                                                                                 cGGGCCATCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGTCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGt                                                                                                                    <  1:380212/251‑1 (MQ=255)                                                                                                                                        tCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACcg                                                                                                                                                                                  >  1:569646/1‑183 (MQ=255)                                                                                                                                          cccgccgTGATAGTCGTTTTCACCCAGATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCCCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGtt                                                                                                           <  2:287687/250‑1 (MQ=255)                                                                                                                                          aCCGCCGTGATCGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGATCTCCCGACCGGGATCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTCTTCCTAAACGGGTCCGGGT‑CCGCTTCATAGGCCCCGTTGTAGcatt                                                                                                           >  2:239715/1‑248 (MQ=255)                                                                                                                                          aCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTAATAGGCCCAGTTGTAGCGtt                                                                                                           >  1:439901/1‑251 (MQ=255)                                                                                                                                                  gATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCATTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACcg                                                                                                                                                                                  <  1:34017/173‑1 (MQ=255)                                                                                                                                                  gATAGTCGTTTTCACCACTATCGAGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCGGACCGGGTTCTTGCGCCCAGTTCTTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGCTCCACCCAGACCTTTACcg                                                                                                                                                                                  >  2:34017/1‑173 (MQ=255)                                                                                                                                                    tAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCCGTTGTCGCGTTTACCCTGt                                                                                                   >  1:343344/1‑249 (MQ=255)                                                                                                                                                      gTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCg                                                                                               <  1:268349/251‑1 (MQ=255)                                                                                                                                                         gTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACcgcg                                                                                                                                                                                <  2:427578/168‑1 (MQ=255)                                                                                                                                                         gTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACcgcg                                                                                                                                                                                >  1:427578/1‑168 (MQ=255)                                                                                                                                                          ttttCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGa                                                                                            >  1:242022/1‑250 (MQ=255)                                                                                                                                                            ttCACCCATATCGCAAGTTTAGGCCGTTCTGGAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTATGCCAGACCGGGTTCTTGGGCCCAGTTACCGAGGTCCAGCGAATTGCGACGGATGGAGCACATGCCGTTGATCAGCGACGATCCACCCAGTCCTTGACCGCGTCCGCACTCCAGGCGGCGGTTATTTATAAACGGGTCAGGTG‑CCGTTTCATAGGCCCGGTTGTAGCGTTTGCCCTGGGGCGGGAAt                                                                                          <  2:297717/250‑1 (MQ=255)                                                                                                                                                                                  gcctgtcgggagtagGGCTGGCAGGCGAGGTAGCTCCAGTTCTCCAGACCGGTTTATTGCGCCCAGGTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTGCCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGTCATCGGGCATCTGGGt                                                                   <  2:497269/241‑1 (MQ=255)                                                                                                                                                                                  gCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTGCCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGGTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGt                                                                   <  2:391106/251‑1 (MQ=255)                                                                                                                                                                                   ccTTGCGGTAGTATGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACTGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTg                                                                  <  1:569712/251‑1 (MQ=255)                                                                                                                                                                                        cGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCTCCCAGTTATTGAGATCCAGCGCATTGCCACGCATGTAGCTCATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCTTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTTCCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGaa                                                             <  1:552870/251‑1 (MQ=255)                                                                                                                                                                                          gtagtaGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGt                                                           <  1:59663/251‑1 (MQ=255)                                                                                                                                                                                            agtagGGCAGGCATTCGAGGTAGCTCCAGTTCTCCAGCCCGGGTTCTTGCGCCCAGTTATTGAGAGCCAGCGCATTGCCACGGATGGAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGGCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGGGCGGAAGta                                                          <  2:158764/250‑2 (MQ=255)                                                                                                                                                                                               agGGCCGGCAGTCGGGGTAGCTCCATTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGACATCCAGCGCATTTCCACGGATGTAGCACATGCCGTTGATCAGCGACGCGCCTCCCTGACCTTTCCCGCTTCCGCCCTCCATGCGGCGGGTATTCATAAACGGTTCAGGTT‑TCGTTTCACAGGCCCAGTTGTAGCGTTTACCCTGCAGCGGGAATGCCAGGTCCGCGGGTATCTggtggcggacgccaaag                                                      <  2:242850/251‑16 (MQ=255)                                                                                                                                                                                                  gcaggtatgcGAGGTAGCTCCATTTCACCAGACCGGTGTCTTGCGCCCCGTTATCGAGTTCCAGCGCACTTCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGCCCTTTTCCGCGTCCGCACTCCATGCGGCGGGTATTCATAACCGGTTCAGGTT‑CCGTGTCACAGGCCCATGTGTAGCGTTTGCCCTGTAGCGGGAATGCCAGGGCAGCGGGCAGCTGTGTTCCGAAGTCAAAGCGa                                                   <  2:489936/242‑1 (MQ=255)                                                                                                                                                                                                    aggcagTCGAGGTCGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGCCCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATa                                                 <  2:199024/251‑1 (MQ=255)                                                                                                                                                                                                         gTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCTTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCg                                            <  1:239715/251‑1 (MQ=255)                                                                                                                                                                                                           gggggTAGCTCCAGTTCTCCGGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGGCCTTTACCGCGTCCGCACTCCATGCGGCGGTTAGTCATAAACGGTTCAGGTT‑CCGTTTCAGAGGCCCAGTTGGAGCGTTGACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCTGAAGTCAAAGCGATAGTCCggg                                          <  2:68604/248‑1 (MQ=255)                                                                                                                                                                                                            gAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGGCCGGgc                                         >  1:196357/1‑251 (MQ=255)                                                                                                                                                                                                              ggTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCTGGTTTCCGTT‑CCTAGGGCCCGTTGTAGCGTTTAACCTGTTGCGGGAATGCCAGGGCAGCGGGCCTCTTGGTGCGGAAGTCCAAGCgttatgcccggccc                                       >  2:145732/1‑238 (MQ=255)                                                                                                                                                                                                                       aGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATTGAGATCCAGCGCAGTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGGTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCa                              <  1:577924/251‑1 (MQ=255)                                                                                                                                                                                                                             ccAGACCTGGTTCTGGCGCCCAGTTATTGAGATCCATCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTGCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTTCCCTGTAGCGGGAATGCCAGGGCAGCTGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCTCTTCAagcagc                        <  2:226818/251‑1 (MQ=255)                                                                                                                                                                                                                              cAGACCTGGTTCTTGCGCCCATTTCTCTTGATCCATCGCAGTGCCTCGGATTTAGCGCATGCCGTTGATCAGCGGCGTTCCACCCAGTCCTTTACCGCGTCCGCACTCCATCCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGTCCCAGTTGTAGCGTTTACCCTGTAGCGTGAATGCCAGGTCAGCTGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCTGGCCGCCCTCTTCAagcagca                       <  2:581775/251‑1 (MQ=255)                                                                                                                                                                                                                                        gcTTGCGCCCAGTTAGTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGt             <  2:22888/250‑1 (MQ=255)                                                                                                                                                                                                                                            gcgcCCAGTTATCGAGACCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGAGCAGCGACGATCCACCCAGACCGGTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATTCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTAtt          <  1:63267/250‑1 (MQ=255)                                                                                                                                                                                                                                              gcCCAGTTATTGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCgg       <  1:182399/251‑1 (MQ=255)                                                                                                                                                                                                                                                ccAGTTAGCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGAt     <  1:309483/251‑1 (MQ=255)                                                                                                                                                                                                                                                   gTTATCGAGACCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACTGAGGTATTCTGATCtt  <  1:450040/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       TCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGG‑TAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT‑CCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTT  >  NC_000913/326689‑327180

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