Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	268,409	C→T	E184E (GAG→GAA)	ykfA ←	CP4‑6 prophage; putative GTP‑binding protein YkfA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	268,409	0	C	T	100.0%	147.4 / NA	41	E184E (GAG→GAA)	ykfA	CP4‑6 prophage; putative GTP‑binding protein YkfA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/29);  total (12/29)

TCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGT  >  NC_000913/268333‑268480                                                                             |                                                                        tcCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCAt                                                                     <  1:62326/81‑1 (MQ=255) tcCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCAt                                                                     <  1:62325/81‑1 (MQ=255)    gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATg                                                                    >  2:62379/1‑79 (MQ=255)       cGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATgcag                                                                 >  1:62382/1‑79 (MQ=255)        ggCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATgcagca                                                               <  1:62327/80‑1 (MQ=255)         gCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAgcagg                                                             >  1:62384/1‑81 (MQ=255)                ccTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGAttt                                                      <  1:62328/81‑1 (MQ=255)                 cTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGAtttt                                                     <  1:62329/81‑1 (MQ=255)                    aaGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTg                                                    <  1:62331/79‑1 (MQ=255)                    aaGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTg                                                    <  2:62330/79‑1 (MQ=255)                     aGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCgg                                                 >  2:62378/1‑81 (MQ=255)                       cccccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCt                                               <  1:62333/81‑1 (MQ=255)                       cccccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCt                                               <  1:62332/81‑1 (MQ=255)                        ccccATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTg                                              <  1:62334/81‑1 (MQ=255)                             ttGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAtt                                         >  2:62373/1‑81 (MQ=255)                                tAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttct                                       <  2:62335/80‑1 (MQ=255)                                  gACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgc                                     <  2:62336/80‑1 (MQ=255)                                   aCGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgctt                                   <  1:62337/81‑1 (MQ=255)                                   aCGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgctt                                   >  1:62402/1‑81 (MQ=255)                                    cGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgcttc                                  >  2:62387/1‑81 (MQ=255)                                     ggACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgcttct                                 <  1:62338/81‑1 (MQ=255)                                         cGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTgcgc                             <  1:62339/81‑1 (MQ=255)                                          gACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGcc                            >  2:62398/1‑81 (MQ=255)                                             acaGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTg                         <  1:62340/81‑1 (MQ=255)                                              caGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTgg                        <  1:62341/81‑1 (MQ=255)                                              caGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTgg                        <  1:62342/81‑1 (MQ=255)                                              caGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTg                         >  1:62385/1‑80 (MQ=255)                                                gCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGAc                      <  2:62344/81‑1 (MQ=255)                                                gCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGAc                      <  1:62343/81‑1 (MQ=255)                                                 cacacaCCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACa                     <  2:62345/81‑1 (MQ=255)                                                    acacCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACt                  <  2:62347/81‑1 (MQ=255)                                                    acacCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACt                  <  1:62346/81‑1 (MQ=255)                                                       tcGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTc               <  1:62348/80‑1 (MQ=255)                                                       ccGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTc               <  1:62349/81‑1 (MQ=255)                                                        cGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTcc              >  2:62376/1‑81 (MQ=255)                                                         ggATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTccc             <  2:62350/81‑1 (MQ=255)                                                            tGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCa          <  1:62351/81‑1 (MQ=255)                                                              gACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCAcc        <  1:62352/81‑1 (MQ=255)                                                                 gggCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCt     >  2:62400/1‑81 (MQ=255)                                                                 gggCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCt     >  2:756287/1‑81 (MQ=255)                                                                    cTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGt  <  1:62353/81‑1 (MQ=255)                                                                             |                                                                        TCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGT  >  NC_000913/268333‑268480

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