Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	68,807	A→G	I414I (ATT→ATC)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	68,807	0	A	G	100.0%	45.6 / NA	15	I414I (ATT→ATC)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/8);  total (7/8)

CATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGC  >  NC_000913/68681‑68927                                                                                                                               |                                                                                                                         cATCACGTAATTAACCGCCATCACCTGATCTGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTGATCAcc                                                                                                                    <  1:234615/133‑1 (MQ=255)                   gATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTGATCACCCCTTTCAGGAGTTGGTTa                                                                                                  >  1:224783/1‑133 (MQ=255)                              cGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTGATCACCCCTTTCAGGCGTTGGTTAGCGTTCGgtgt                                                                                       >  2:461548/1‑133 (MQ=255)                                    aGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTGATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTgcggcg                                                                                 <  2:326871/133‑1 (MQ=255)                                           cATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTGATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTa                                                                          >  1:614809/1‑133 (MQ=255)                                                      gTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTGATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAg                                                               <  2:495343/133‑1 (MQ=255)                                                       tGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTGATCACCCCTTTCAGGCGTAGGTTAGCGTTCGATGTGCGGCGGCCCTTAAACCAGTCGAGc                                                              >  2:253706/1‑133 (MQ=255)                                                              aaGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTGATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCg                                                       <  1:431469/133‑1 (MQ=255)                                                                             aTCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTGATCACCCCTTTCAGGCGTTAGTTAGCGTTCGATGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCagag                                        >  1:607518/1‑133 (MQ=255)                                                                                                      cGTCGGTAGCTAGGTTAAGATCGGTGATTACCTCTTTCAGGGGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGt               <  2:258509/133‑1 (MQ=255)                                                                                                      cGTCGGTAGCGAGGTTAAGATCGGTGATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGt               >  2:347334/1‑133 (MQ=255)                                                                                                                gAGGTTAAGATCGGTGATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGa     <  1:603774/133‑1 (MQ=255)                                                                                                                 aGGTTAAGGTCGGTGGTCCCCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGaa    <  1:253706/133‑1 (MQ=255)                                                                                                                  ggTTAAGATCGGTGATCACCCCTTTCAGGCGTTAGTTAGCGTTCGATGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTAGCCCATGCTTCGGTCAGCGCCGGAAg   <  2:607518/133‑1 (MQ=255)                                                                                                                   gTTAAGATCGGTGATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTTGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGc  >  2:371836/1‑133 (MQ=255)                                                                                                                               |                                                                                                                         CATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGC  >  NC_000913/68681‑68927

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