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%	158.1 / NA	42	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 (21/21);  total (21/21)

ATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTC  >  NC_000913/268332‑268486                                                                              |                                                                              atcCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCa                                                                            >  1:52810/1‑81 (MQ=255)    cGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATg                                                                          <  1:52775/80‑1 (MQ=255)     gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATgca                                                                        <  2:52776/81‑1 (MQ=255)     gCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATg                                                                          >  2:52813/1‑79 (MQ=255)      cTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATgcag                                                                       <  2:52778/81‑1 (MQ=255)      cTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATgcag                                                                       <  2:52777/81‑1 (MQ=255)       tCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATgcag                                                                       >  1:52814/1‑80 (MQ=255)       tCGGCCATCACCTTAAG‑CCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATgcag                                                                       >  2:52819/1‑79 (MQ=25)        cGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATgcag                                                                       >  2:52826/1‑79 (MQ=255)         ggCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGcagcag                                                                    <  1:52779/81‑1 (MQ=255)         ggCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGcagcag                                                                    >  1:52820/1‑81 (MQ=255)           ccATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAgcaggc                                                                  <  1:52781/81‑1 (MQ=255)           ccATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAgcagg                                                                   >  1:52816/1‑80 (MQ=255)           ccATCACCTTAAGCCCCAATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAgcaggc                                                                  <  2:52780/81‑1 (MQ=255)             aTCACCTTAAACCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAgcaggcag                                                                <  2:52782/81‑1 (MQ=255)               cACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGAt                                                              >  1:52823/1‑81 (MQ=255)                aCCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGAtt                                                             <  2:52783/81‑1 (MQ=255)                 ccTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGAttt                                                            <  1:52784/81‑1 (MQ=255)                     aaGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGc                                                         <  1:52785/80‑1 (MQ=255)                      aGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCgg                                                       >  2:52822/1‑81 (MQ=255)                      aGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCgg                                                       <  2:52786/81‑1 (MQ=255)                            cATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGata                                                 <  1:52787/81‑1 (MQ=255)                             aTTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGatat                                                >  1:52815/1‑81 (MQ=255)                             aTTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGatat                                                <  2:52788/81‑1 (MQ=255)                                gTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCGGCAGGCAGATTTTGCGGCTGATAttct                                             >  1:52818/1‑81 (MQ=255)                                gTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttc                                              <  2:52789/80‑1 (MQ=255)                                    aCGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgctt                                         >  1:52828/1‑81 (MQ=255)                                    aCGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgctt                                         >  2:52833/1‑81 (MQ=255)                                       gACCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgcttctg                                      <  2:52790/81‑1 (MQ=255)                                       gACCGAAACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgcttctg                                      <  2:52791/81‑1 (MQ=255)                                        aCCGACACAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATAttctgcttctgc                                     >  1:52821/1‑81 (MQ=255)                                            acacAGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCg                                 <  1:52792/81‑1 (MQ=255)                                              acaGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTg                               >  1:52842/1‑81 (MQ=255)                                                aGCACACACCGGATGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGa                             >  2:52825/1‑81 (MQ=255)                                                             tGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCa                >  1:52835/1‑81 (MQ=255)                                                             tGGACGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCa                <  1:52793/81‑1 (MQ=255)                                                                 cGGGCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGc            >  2:52824/1‑81 (MQ=255)                                                                   ggCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTg          >  1:52844/1‑81 (MQ=255)                                                                   ggCTGGAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTg          <  2:52794/81‑1 (MQ=255)                                                                       ggAACAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAg      >  2:52843/1‑81 (MQ=255)                                                                           cAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTc  >  1:52831/1‑81 (MQ=255)                                                                           cAGTTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTc  <  1:52795/81‑1 (MQ=255)                                                                              |                                                                              ATCCGCTCGGCCATCACCTTAAGCCCCCATTGTAGACGGACCGACACAGCACACACCGGATGGACGGGCTGGAACAGCTCATGCAGCAGGCAGATTTTGCGGCTGATATTCTGCTTCTGCGCCGTGGACAACTGTCCCCCACCGCTGGTAGGTTC  >  NC_000913/268332‑268486

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