Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,070,576	G→A	P394L (CCG→CTG)	yihS ←	sulfoquinovose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,070,576	0	G	A	100.0%	122.1 / NA	34	P394L (CCG→CTG)	yihS	sulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/16);  total (18/16)

TTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTA  >  NC_000913/4070503‑4070653                                                                          |                                                                              ttttCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATAc                                                                         <  1:820912/80‑1 (MQ=255)   ttCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATAcgc                                                                       <  1:820913/80‑1 (MQ=255)    tCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCgg                                                                     >  2:820964/1‑81 (MQ=255)      aGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGa                                                                   >  2:820960/1‑81 (MQ=255)       gATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGAt                                                                  >  2:820954/1‑81 (MQ=255)           aGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCAcc                                                              <  2:820914/81‑1 (MQ=255)           aGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCAc                                                               >  1:820958/1‑80 (MQ=255)           aGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCAc                                                               <  2:820915/80‑1 (MQ=255)            gTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCAc                                                               >  2:820956/1‑79 (MQ=255)                tttCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGc                                                          >  2:820968/1‑80 (MQ=255)                 ttCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCaa                                                        >  1:820957/1‑81 (MQ=255)                     cATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATgca                                                    <  2:820916/81‑1 (MQ=255)                      aTTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATgcag                                                   <  1:820917/81‑1 (MQ=255)                      aTTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATgcag                                                   >  2:820969/1‑81 (MQ=255)                        tAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATgcag                                                   >  1:820962/1‑79 (MQ=255)                               cagcagACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGAt                                           <  1:820919/80‑1 (MQ=255)                               cagcagACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGAt                                           <  1:820918/80‑1 (MQ=255)                                 gcagACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATaaa                                        >  1:820961/1‑81 (MQ=255)                                   agACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAata                                      >  1:820970/1‑81 (MQ=255)                                         gCTGCAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGt                                >  1:820981/1‑81 (MQ=255)                                             cAACCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGttt                              <  2:820920/79‑1 (MQ=255)                                              aaCCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGcc                           <  1:820922/81‑1 (MQ=255)                                              aaCCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGcc                           <  1:820921/81‑1 (MQ=255)                                               aCCGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCg                          <  1:820923/81‑1 (MQ=255)                                                ccGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGt                         >  2:820963/1‑81 (MQ=255)                                                 cGCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTc                        >  2:820967/1‑81 (MQ=255)                                                  gCTGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTcc                       <  1:820924/81‑1 (MQ=255)                                                    tGGAGCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTcc                       >  1:820965/1‑79 (MQ=255)                                                        gCCATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACt                 >  2:820977/1‑81 (MQ=255)                                                         ccATGCCGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACtt                <  1:820925/81‑1 (MQ=255)                                                              ccGGGCGCTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTtggtg           <  2:820926/81‑1 (MQ=255)                                                                   cgcTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGAc      >  2:820974/1‑81 (MQ=255)                                                                     cTAACAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCt    <  1:820927/81‑1 (MQ=255)                                                                       aaCAGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTa  >  1:820972/1‑81 (MQ=255)                                                                          |                                                                              TTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTA  >  NC_000913/4070503‑4070653

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