Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,193,996	T→G	D70A (GAT→GCT)	thiE ←	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,193,996	0	T	G	100.0%	187.6 / NA	54	D70A (GAT→GCT)	thiE	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/27);  total (27/27)

GGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTCTGGCACGGAAGTTCTCCGACATATCCGCCATTCCCATTTCAATAG  >  NC_000913/4193726‑4194286                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                   ggATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcgg                                                                                                                                                                                                                                                                       <  1:168430/300‑1 (MQ=255)    tAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcggc                                                                                                                                                                                                                                                                      >  1:701768/1‑298 (MQ=255)                        attccGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAGGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCTCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATGTCGTCATGTGNCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccgcc                                                                                                                                                                                                                                                  <  2:207961/295‑1 (MQ=255)                                 cAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGc                                                                                                                                                                                                                                       >  1:251376/1‑300 (MQ=255)                                      gTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCAAGTAAGCGTTGATAAACAATTGCGCGTTATAGCGGCGGCACAGCGCAATTGCCGCCACGACATCGGCttac                                                                                                                                                                                                                                   >  2:186873/1‑297 (MQ=255)                                      gTTCCAGCCCCTGCGGNGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGNGGCCCAGCGCAATTGCCGNCACGACATCGGCTTCCa                                                                                                                                                                                                                                  >  1:593203/1‑300 (MQ=255)                                       ttCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCTCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCAc                                                                                                                                                                                                                                 <  2:468796/300‑1 (MQ=255)                                        tCCAGCCCCTGCGGTGAAGAAGGCATCTTTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCATGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCAc                                                                                                                                                                                                                                 <  2:58869/299‑1 (MQ=255)                                                     gTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATCCGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGTGGCGGCCCAACGCAATTGCCGCCCCGACATCGGCTTCCACCTCTTCATcgcgac                                                                                                                                                                                                                   >  2:67330/1‑298 (MQ=255)                                                     gTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGcc                                                                                                                                                                                                                   >  1:62265/1‑300 (MQ=255)                                                     gTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTACCCCTCTTCATCgcgc                                                                                                                                                                                                                    >  2:65542/1‑299 (MQ=255)                                                                             gCGTCTGGAACACATGTCCCAGCGCGATATAAGATGTGCGTGCTGCCAGCGCGGCGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCGTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGagag                                                                                                                                                                                           <  2:105881/300‑1 (MQ=255)                                                                                     aaCACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGtacgtacg                                                                                                                                                                                   <  1:540442/300‑1 (MQ=255)                                                                                       cacaTGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGNCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGcc                                                                                                                                                                                 >  2:518137/1‑300 (MQ=255)                                                                                        acaTGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGcc                                                                                                                                                                                 <  2:431986/299‑1 (MQ=255)                                                                                        acaTGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGGTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGcc                                                                                                                                                                                 <  1:415387/299‑1 (MQ=255)                                                                                            gTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCAt                                                                                                                                                                            >  1:183208/1‑300 (MQ=255)                                                                                                     cgATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGAc                                                                                                                                                                   >  1:1808/1‑300 (MQ=255)                                                                                                     cgATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCAAGACATCGGCTTCCACCTCTTCATCGCGCCGATCCTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGAc                                                                                                                                                                   >  2:51930/1‑300 (MQ=255)                                                                                                       atatAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGt                                                                                                                                                                 >  2:683787/1‑300 (MQ=255)                                                                                                                   cGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGAGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAAAAGACGTTCGATCCa                                                                                                                                                         >  2:158894/1‑296 (MQ=255)                                                                                                                       ctgcCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGc                                                                                                                                                 >  1:433812/1‑300 (MQ=255)                                                                                                                        tgcCAGCGCGACGTCGATTTCCATATCGTCATGTTTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGc                                                                                                                                                 <  2:401606/299‑1 (MQ=255)                                                                                                                        tgcCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGc                                                                                                                                                 <  2:401617/299‑1 (MQ=255)                                                                                                                        tgcCAGCGCGAAGGCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGc                                                                                                                                                 <  2:331911/299‑1 (MQ=255)                                                                                                                        tgcCAACGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCATGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGc                                                                                                                                                 <  2:248839/299‑1 (MQ=255)                                                                                                                            aacgcggctttgATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGTCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACATTCGATCCACTGTACGCTGTcc                                                                                                                                            <  2:99119/289‑1 (MQ=255)                                                                                                                                 taCGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGGGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATTCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCcaccac                                                                                                                                       <  2:238839/299‑1 (MQ=255)                                                                                                                                 gacttCGATTTCCAGATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGTTTGCAAATCTTCCTGCCCCAAATGGACACCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCcaccac                                                                                                                                       <  2:238831/296‑1 (MQ=255)                                                                                                                                 gACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCcaccac                                                                                                                                       <  2:363265/300‑1 (MQ=255)                                                                                                                                  aCGTCGATTTCCATATCATCATGTTTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAGATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCATTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACGGACGTTCGATCCACTGTActc                                                                                                                                                 <  2:53209/289‑3 (MQ=255)                                                                                                                                       gATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGt                                                                                                                                  <  1:63012/299‑1 (MQ=255)                                                                                                                                            ccATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTc                                                                                                                            >  1:572450/1‑300 (MQ=255)                                                                                                                                                   gTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCNCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAAc                                                                                                                      <  1:502989/299‑1 (MQ=255)                                                                                                                                                                       gCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAAt                                                                                                 >  1:273894/1‑300 (MQ=255)                                                                                                                                                                       gCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAAt                                                                                                 >  2:487049/1‑300 (MQ=255)                                                                                                                                                                       gCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCAACGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAAt                                                                                                 >  1:683929/1‑300 (MQ=255)                                                                                                                                                                             ggCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCt                                                                                           >  1:453297/1‑300 (MQ=255)                                                                                                                                                                                       gcgGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATNGCGCCGATCTTTGATGCGTAGCTGGAGAGTACNTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCAc                                                                                 <  2:519970/300‑1 (MQ=255)                                                                                                                                                                                                aTTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACNTCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCAAGCTtcctcc                                                                        >  2:507647/1‑300 (MQ=255)                                                                                                                                                                                                      aTCGGTGGCTTGCAAATCTTCCTGCCCCAAATGTATGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACTAAAAGGTACAGGAGGAGAATCAGGCTGATACATCACGCTTCCTCCTTACGc                                                                  <  2:669357/300‑1 (MQ=255)                                                                                                                                                                                                      aTCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGc                                                                  <  1:43822/300‑1 (MQ=255)                                                                                                                                                                                                                      tCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGc                                                  <  1:476971/300‑1 (MQ=255)                                                                                                                                                                                                                        ttCCTGCCCCAAATGGCCGCCATACACCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGACGGCCCAGCGCAATTGCCGCCACGACATAGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCAACAGATACAGTCCTGAACGAAAAGGTACAGGAGGAAAACCAGGc                                                                                            >  1:184872/1‑256 (MQ=255)                                                                                                                                                                                                                        ttCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCct                                                >  1:185177/1‑300 (MQ=255)                                                                                                                                                                                                                                           ccATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTCTGGCACGGAAGTtctac                              >  1:163483/1‑297 (MQ=255)                                                                                                                                                                                                                                           ccATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTCTGGCACGGAAGTTCTcc                              >  1:478978/1‑299 (MQ=255)                                                                                                                                                                                                                                           ccATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCATTACGCAGGTAGATTTCTCCGCCTCTGGCACGGAAGTTCTcc                              >  1:478915/1‑299 (MQ=255)                                                                                                                                                                                                                                           ccATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCNGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTCTGGCACGGAAGTTCTcc                              >  2:504860/1‑299 (MQ=255)                                                                                                                                                                                                                                             aTACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTCTGGCACGGAAGTTCTCCGa                            <  2:296925/299‑1 (MQ=255)                                                                                                                                                                                                                                                        gCTTGATCGCCAGCCGNCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCNGCCTCTGGCACGGAAGTTNTCCGACATATCCGCCAt                >  1:587271/1‑300 (MQ=255)                                                                                                                                                                                                                                                             aTCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTGAGAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCACGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGAttt                                                        <  2:343949/255‑1 (MQ=255)                                                                                                                                                                                                                                                                       gccAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTCTGGCACGGAAGTTCTCCGACATATCCGCCATTCCCATTTCAATAg  <  1:186873/299‑1 (MQ=255)                                                                                                                                                                                                                                                                       gccAGTAAGCGTGGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTCTGGCACGGAAGTTCTCCGACATATCCGCCATTCCCATTTCAATAg  <  1:456430/299‑1 (MQ=255)                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                   GGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTCTGGCACGGAAGTTCTCCGACATATCCGCCATTCCCATTTCAATAG  >  NC_000913/4193726‑4194286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑