Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,308,804	G→A	F222F (TTC→TTT)	alsC ←	D‑allose ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,308,804	0	G	A	100.0%	28.8 / NA	13	F222F (TTC→TTT)	alsC	D‑allose ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/5);  total (8/5)

CCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGCCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGGAATTTCACGTCAATCCCGGAATAGAACGCCGAATTTTTGTTGCCGCCCAGTGCGTAGATGTTGCGCCCGAGCCGCATACGCGTTGTCAGAAACCAAAGGATGAGCGCGACAATTAGTGAGAAGATAACGGGGACAGGTATCCCAATTACGCTGGCGGCAAAGAAGTTCACGAAGTC  >  NC_000913/4308569‑4308980                                                                                                                                                                                                                                            |                                                                                                                                                                                 ccTGCAAAATATTCAGACCGTTGTTTATGGTGCCGATGCTCAACCCGCCAATCAACACAGCGAAAGTGCGACCCTTGCCGCCGAAAAAACTGGTGCCGC‑AATGATAGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGCCCTGCACAAACACAGGAGATGATAACCACCACAAGGATGTGAAATTTCACGTCAATcc                                                                                                                                                                  >  1:307705/1‑251 (MQ=255)                         tgATGGTCCCGATTATCAACCCGCCAATCACCACAGAGAAAAGGCGCCCCTTGCCGCCGAATAAACTGGGGCCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACTACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGAATTTTCACGTCAATCCCACAATAGAACGCCGAAttttt                                                                                                                                            <  1:271772/250‑1 (MQ=255)                                  ngatgatCAACCCGCCAATCACCACAGAGCAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGCCGCCCATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCCTACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGAAGACGCCTGCCAGACCTGCACAAACACCGGAGCTGATAAACACCACCAGGATGTGAAATTTCCCGTCAATCCCGGAATAGAACGCCGAATTTTTGTTGCCGccc                                                                                                                                  >  1:38804/2‑251 (MQ=255)                                                 cAATCACCACAGAGAAAATGCGCCCCGTGCCGCCGAAGAAACTGGTGCCGCCAATGATGGCGCTGGCAAGGGCATAGGTTTCAAAACCCATACCGGCAAGCGGGTCTGCGGCACCTAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGAAATTTCACGGCAATCCCGGAAGAGAACGCCGAATTTTTGTTGCCGCCCAGTGCGTAGATGTTg                                                                                                                   <  1:183988/251‑1 (MQ=255)                                                    tCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGCCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGAAATTTCACGTCAATCCCGGAATAGAACGCCGAATTTTTGTTGCCGCCCAGTGCGTAGATGt                                                                                                                     <  1:149128/246‑1 (MQ=255)                                                            gagaAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGCCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGAAATTTCACGTCAATCCCGGAATAGAACGCCGAATTTTTGTTGCCGCCCAGTGCGTAGATGTTGCGCCCGAGCCg                                                                                                        >  1:22339/1‑251 (MQ=255)                                                                                        aaCTGGTGCCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGAAATTTCACGTCAATCCCGGAATAGAAc                                                                                                                                                       <  1:154696/176‑1 (MQ=255)                                                                                             gTGCCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGAGGATCTGCGGCAACGAGTACTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACAGGAGATGATAAACACCACCAGGATGTGAAATTTCACg                                                                                                                                                                         >  1:215465/1‑153 (MQ=255)                                                                                                                             gttttCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGAAGATGATAAACACCACCAGGATGTGAAATTTCACGTCAATCCCGGAATAGAACGCCGAATTTTTGTTGCCGCCCAGTGCGTAGATGTTGCGCCCGAGCCGCATACGCGTTGTCAGAAACCAAAGGATGAGCGCGACAATTAGTGAGAAGATAACGGGGACAGGt                                        >  1:306390/3‑250 (MQ=255)                                                                                                                                   aaaaCCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGAAAtt                                                                                                                                                                              >  1:123384/1‑110 (MQ=255)                                                                                                                                             ccGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGAAATTTCACGTCACTCCCGGAATAGAACGCCGAATTTTTGTTGCCGCCCAGTGCGTAGATGTTGCGCCCGAGCCGCATACGCGTTGTCAGAAACCAAAGGATGAGCGCGACAATTAGTGAGAAGATAACGGGGACAGGTATCCCAATTACGCTggc                       >  1:53173/1‑251 (MQ=255)                                                                                                                                                                  ccGATTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACCCCACCAGGATGTGAAATTTCACGTCAATCCCGGAATAGAACGCCGAATTTTTGTTGCCGCCCAGTGCGTAGATGTTGCGCCCGAGCCGCATACGCGTTGTCAGAAACCAAAGGATGAGCGCGACAATTAGTGAGAAGATAACGGGGACAGGTATCCCAATTACTCTGGCGGCAAAGAAGTTCACGAAGtn  <  1:262958/251‑2 (MQ=255)                                                                                                                                                                                                                                  ccaGGATGTGAAATTTCACGTCAATCCCGGAATAGAACGCCGAATTTTTGTTGCCGccc                                                                                                                                  >  1:275664/1‑59 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                 CCTGCAAAATATTCAGACCGTTGTTGATGGTGCCGATGATCAACCCGCCAATCACCACAGAGAAAATGCGCCCCTTGCCGCCGAAGAAACTGGTGCCGCCAATGATGGCGCTGGCAATGGCATAGGTTTCAAAACCCATACCGGCAAGCGGTTCTGCGGCACCGAGTCGTGCAGTTGAGACGACGCCTGCCAGACCTGCACAAACACCGGAGATGATAAACACCACCAGGATGTGGAATTTCACGTCAATCCCGGAATAGAACGCCGAATTTTTGTTGCCGCCCAGTGCGTAGATGTTGCGCCCGAGCCGCATACGCGTTGTCAGAAACCAAAGGATGAGCGCGACAATTAGTGAGAAGATAACGGGGACAGGTATCCCAATTACGCTGGCGGCAAAGAAGTTCACGAAGTC  >  NC_000913/4308569‑4308980

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