Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,552,362	T→C	intergenic (+371/+36)	fdnI → / ← yddM	formate dehydrogenase N subunit gamma/DNA‑binding transcriptional regulator YddM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,552,362	0	T	C	100.0%	161.1 / NA	48	intergenic (+371/+36)	fdnI/yddM	formate dehydrogenase N subunit gamma/DNA‑binding transcriptional regulator YddM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/18);  total (30/18)

TCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTC  >  NC_000913/1552117‑1552597                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            accagtaaaTTGCAGTGAACGTATAGGAATGAGAAGCGGAGGAAATCAGACAATGTTGTGTTAGCCATCAGTTGCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAgat                                                                                                                                                                                                                                        <  2:307109/243‑3 (MQ=11)    aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTACCGCATCCGGCACTGTTGCGTTTGTCATCAGATTCAAATGGCGATGAAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCATTATCCGGCATGAACACACATCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGAGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGtcct                                                                                                                                                                                                                                      >  1:249242/1‑247 (MQ=18)     aTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTa                                                                                                                                                                                                                                     >  1:70809/1‑250 (MQ=255)            gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAATGCCTTATTCGGCATGTAAACCCTTTAAATCCATTAGATTGCAGTGAACGTGTAGGCCTGATAACCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTTAATGgcgc                                                                                                                                                                                                                             >  2:353359/1‑251 (MQ=18)                 gAACGTTTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTtat                                                                                                                                                                                                                        <  1:176531/251‑1 (MQ=255)                 gAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTtat                                                                                                                                                                                                                        >  1:182274/1‑251 (MQ=255)                               gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATACAATAGGTTGCAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCCCCCGTTCTAAATGGCGCTTTATAAAGGGCCAttttt                                                                                                                                                                                                          >  2:227138/1‑251 (MQ=18)                               gATAAGCGGAGCGCATCAGGCAATTTTGCGTGTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCTTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCTGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTATCGCAGCATGCAATGTTGCGTTTTTCACCAGTTGTAAATGGCGGTTTAGAAAGTGCCAtttt                                                                                                                                                                                                           <  2:441487/250‑1 (MQ=17)                                  aaGCGTAGCGAATCAGGTAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa                                                                                                                                                                                                       <  2:10644/251‑1 (MQ=255)                                      gTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTgt                                                                                                                                                                                                   >  1:74336/1‑250 (MQ=255)                                       tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGACCGTGTAAGCCTGATAAGCGTAGCGCCTCAGGCAATGTTGCGTTTGTCACCAGTTCCAAATGGCGCTTTTTTAAGTGCCATTTTTTTAtttcg                                                                                                                                                                                                  >  2:145192/1‑248 (MQ=17)                                       tAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg                                                                                                                                                                                                  >  1:17536/1‑251 (MQ=255)                                                   cAATGTTGCGTTTGTCATCAGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGc                                                                                                                                                                                      <  2:102876/251‑1 (MQ=255)                                                                            aaaTGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGtttt                                                                                                                                                             >  1:38662/1‑251 (MQ=255)                                                                                      gTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTtctgcctct                                                                                                                                                   >  1:358446/1‑251 (MQ=255)                                                                                        aaaaGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCCGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTCTTGCGTAACCAGACGGCGTAATCGCGACAAATCCACTGTTTTTtctgcctctgc                                                                                                                                                 >  1:25625/1‑251 (MQ=255)                                                                                              cGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACt                                                                                                                                           >  1:239054/1‑251 (MQ=255)                                                                                                       cATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtt                                                                                                                                  <  2:118241/251‑1 (MQ=255)                                                                                                         tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGtttt                                                                                                                                >  1:71796/1‑251 (MQ=255)                                                                                                                     cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGCTTCAGcc                                                                                                                    >  2:125619/1‑251 (MQ=255)                                                                                                                        aGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCcac                                                                                                                  <  2:468187/250‑1 (MQ=255)                                                                                                                                  cATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTgcggcg                                                                                                        >  1:144688/1‑250 (MQ=255)                                                                                                                                             cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATc                                                                                             >  1:366071/1‑250 (MQ=255)                                                                                                                                                     aCGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGCTTGCAGTGAACGTGTAGGCCTGATACGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGAGGCGAACTGCCGATCACCACGGa                                                                                     >  1:406901/1‑250 (MQ=255)                                                                                                                                                              ccGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTGGTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATc                                                                           <  2:331552/251‑1 (MQ=255)                                                                                                                                                                     gAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCAttt                                                                    >  1:103691/1‑251 (MQ=255)                                                                                                                                                                       aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTAGAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACAGCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATTACCACGGAAAGCTTAATCGCCATTTCt                                                                  <  1:354889/251‑1 (MQ=255)                                                                                                                                                                                aaaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCTTCaaa                                                         <  1:145192/251‑1 (MQ=255)                                                                                                                                                                                 aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATCAGCGTAGCGCATCAGGCACTGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAATGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGCCTCCAGGCGTTTTCCAGCTTCAGCCACATTTGCGGCGAACTGCCGCTCACCACGCAAAGTTTAATCGCCCTTCCTGTCGTCaa                                                          >  2:105670/1‑249 (MQ=37)                                                                                                                                                                                  aTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAgc                                                       >  2:349765/1‑251 (MQ=255)                                                                                                                                                                                    ccAATAGATTGCAGTGAACGTGTCGGCCTGATAAGCGTAGCGCATCAGGCACTGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAgctg                                                     >  2:51984/1‑251 (MQ=255)                                                                                                                                                                                      aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAgctgct                                                   <  2:76994/251‑1 (MQ=255)                                                                                                                                                                                       aTAGATTGCAGTGAACGTGTAGGCCTGCTAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTCATCGCCATTTCTGGCGTCAAAGCTGCtt                                                  >  2:421578/1‑251 (MQ=255)                                                                                                                                                                                           aTTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTcc                                              <  1:421578/251‑1 (MQ=255)                                                                                                                                                                                            ttGCAGTGAACGGGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCAAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTTTGGCGTCAAAGCTGCTTTTCCt                                             >  1:43078/1‑251 (MQ=255)                                                                                                                                                                                                 gTGAACGTGTAGGCCTGATAAGCGTAGCGCATCCGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCCGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTc                                          >  2:27339/1‑249 (MQ=255)                                                                                                                                                                                                       gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCg                                  >  1:159980/1‑251 (MQ=255)                                                                                                                                                                                                       gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGcc                                                                                                                    >  1:300950/1‑169 (MQ=255)                                                                                                                                                                                                       gtgtAGGCCTGAGAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTTCCAGACTCCAGGCGTTTTGCAGATTCAGcc                                                                                                                    <  2:300950/169‑1 (MQ=255)                                                                                                                                                                                                                    aaGCGTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGa                      <  2:366071/250‑1 (MQ=255)                                                                                                                                                                                                                        gTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATc                                   <  1:156627/233‑1 (MQ=255)                                                                                                                                                                                                                        gTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATc                                   >  2:156627/1‑233 (MQ=255)                                                                                                                                                                                                                        gTAGCGCATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACGCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAg                 <  2:71796/251‑1 (MQ=255)                                                                                                                                                                                                                           gcgcATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGcg               >  1:403490/1‑249 (MQ=255)                                                                                                                                                                                                                              cATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAttt                                                                                                              <  2:266889/152‑1 (MQ=255)                                                                                                                                                                                                                              cATCAGGCAATGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACAttt                                                                                                              >  1:266889/1‑152 (MQ=255)                                                                                                                                                                                                                              cATCAGGCAAGGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCGAttt           <  2:403490/251‑1 (MQ=255)                                                                                                                                                                                                                                       aTGTTGCGTTTGTCACCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCCTTGctc  >  2:265128/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            TCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTAATCGCGACACATCCACTGTTTTTTCTGCCTCTGCCAGACTCCAGGCGTTTTGCAGATTCAGCCACATTTGCGGCGAACTGCCGATCACCACGGAAAGTTTAATCGCCATTTCTGGCGTCAAAGCTGCTTTTCCTGTCAGCAATCGGCTTGCCGTTGAGGGAGCAATTTCCATTGCTC  >  NC_000913/1552117‑1552597

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