Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,490,095	G→A	intergenic (+424/‑807)	aldA → / → cybB	aldehyde dehydrogenase A/superoxide oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,490,095	0	G	A	100.0%	54.2 / NA	20	intergenic (+424/‑807)	aldA/cybB	aldehyde dehydrogenase A/superoxide oxidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/4);  total (16/4)

TGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCG  >  NC_000913/1489862‑1490309                                                                                                                                                                                                                                          |                                                                                                                                                                                                                       tGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTg                                                                                                                                                                                                       >  1:123775/1‑251 (MQ=255)   ccAATGATATAGGACGAGACTATTTCTTCATCGGTATAACCCAATGACTCGTTCTTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGTATGCCCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGCCAAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGcc                                                                                                                                                                                                     >  2:256765/1‑251 (MQ=255)      aTGATCTCGGAAGCGACTATTTCTTCATCGGTCTAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGCGAATCGATACCAGGTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGTAATGATATTTTCTGCCGCt                                                                                                                                                                                                  >  2:131109/1‑251 (MQ=255)        gATATCGGAAGAGCCTATTTCTTCATCGGTATAACCAAATGACTCGCTATTGGTGGTTGCTTGTTTAAGAGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGCATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCCGACACATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTgc                                                                                                                                                                                                >  1:182151/1‑251 (MQ=255)                   agaCTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCAc                                                                                                                                                                                     >  1:76985/1‑251 (MQ=255)                              tcATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCttt                                                                                                                                                                           >  1:11698/1‑250 (MQ=255)                               cATCGGTATAACCAAATGACTCGTTATTGGAGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTAc                                                                                                                                                                         >  1:214966/1‑251 (MQ=255)                                                                ttGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGcc                                                                                                                                          >  1:38455/1‑249 (MQ=255)                                                                  gCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGt                                                                                                                                       >  1:217234/1‑250 (MQ=255)                                                                                           ctTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCAAGTATAGGCATGAATGGTCGTCATCGTGc                                                                                                             >  2:107323/1‑251 (MQ=255)                                                                                             tCGGCAGTCTCTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTTTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTg                                                                                                              <  1:83222/248‑1 (MQ=255)                                                                                                  aGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGATGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGGGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTc                                                                                                      >  2:199804/1‑251 (MQ=255)                                                                                                     cACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGTATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGTGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTAt                                                                                                   <  1:256765/251‑1 (MQ=255)                                                                                                                  gaATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCACCAGTAGTGTAGGGAAGGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTAGTCATCGTGCCGACTTCTATCCCGAAACTGTc                                                                                       >  1:146936/1‑250 (MQ=255)                                                                                                                      cGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGCATGGTCGTCATCGTGCCGACTTCTATCCCGACACTGTCATGc                                                                                   >  2:93629/1‑250 (MQ=255)                                                                                                                          aCCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGc                                                                              <  2:123775/251‑1 (MQ=255)                                                                                                                                         gAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGgcgc                                                               >  1:116542/1‑251 (MQ=255)                                                                                                                                                                    gcATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGAACAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGGGGTGCATGACGCCAc                                    >  2:47316/1‑251 (MQ=255)                                                                                                                                                                                     gCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCg                    <  1:91109/250‑1 (MQ=255)                                                                                                                                                                                                      tCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTgtcg  >  1:164309/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                       TGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTAGTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCG  >  NC_000913/1489862‑1490309

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