Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,197,739	A→G	S167P (TCC→CCC)	ymfE ←	e14 prophage; putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,197,739	0	A	G	100.0%	262.1 / NA	73	S167P (TCC→CCC)	ymfE	e14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (41/32);  total (41/32)

CTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCT  >  NC_000913/1197643‑1197832                                                                                                 |                                                                                              ctGATAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACAGGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAt                                                                                            <  1:457126/100‑1 (MQ=255) ctGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAACAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAt                                                                                            >  2:860054/1‑100 (MQ=255)  tGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAtt                                                                                           <  2:1203071/100‑1 (MQ=255)   gggaAGGATAAAATGACAAAAAAGAATAAAGACAGGAAAGAACCTTCAACATGAAAAATATCAAGTTGTTTGAAAAAAAAGATTATTAGGAAGGGAATTa                                                                                          <  2:123525/98‑1 (MQ=255)      aaTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATg                                                                                       <  1:754871/100‑1 (MQ=255)      aaTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATg                                                                                       <  1:755041/100‑1 (MQ=255)       aTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                      <  1:389751/100‑1 (MQ=255)        tGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                     >  1:1144280/1‑100 (MQ=255)        tGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                     >  1:745187/1‑100 (MQ=255)         gATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCaa                                                                                    >  1:528637/1‑100 (MQ=255)          aTAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                   <  2:1165179/100‑1 (MQ=255)           tAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                   >  2:757236/1‑99 (MQ=255)            aaaaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCacagc                                                                                 >  1:218712/1‑96 (MQ=255)             aaaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTAt                                                                                >  2:275553/1‑100 (MQ=255)             aaaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTAt                                                                                >  1:318028/1‑100 (MQ=255)               aTTACAAAAAAGAGTAAAGACAGCAAAGAACGTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                              >  1:1057277/1‑100 (MQ=255)               aTTACAAAAAAGAATAAAGACAGCAAAAAACCTTCAACATGAAAAATATCCATATGTATGAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                              <  2:207171/100‑1 (MQ=255)                    aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAt                                                                         >  2:768015/1‑100 (MQ=255)                       aaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCa                                                                      <  2:798318/100‑1 (MQ=255)                        aaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaa                                                                     >  1:60422/1‑100 (MQ=255)                        aaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaa                                                                     <  2:983776/100‑1 (MQ=255)                           aaTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaa                                                                  >  1:580107/1‑100 (MQ=255)                            aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaaa                                                                 >  2:673355/1‑100 (MQ=255)                            aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCAATTTGTTAGCAAAAAAAGATTAATAGGCAGGGAATTAATGCAATTATCTAAAATTCaaaaaa                                                                 >  2:1172329/1‑100 (MQ=255)                              aaaGACAGCAAAGAACCTTCAACATGAAAAGTAACCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAGTCAAAAAata                                                               >  2:1029352/1‑100 (MQ=255)                              aaaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGAGAATTAATGCAATTATCGAAAATTCAAAAAata                                                               >  2:279792/1‑100 (MQ=255)                                aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATc                                                             <  2:630122/100‑1 (MQ=255)                                 gACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATcc                                                            >  1:535975/1‑100 (MQ=255)                                  aaagaaaaGAACCTTCAAAATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATAAAAAAAAAAAAAAATATCCa                                                           <  2:661443/95‑1 (MQ=255)                                     gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                        >  1:901307/1‑100 (MQ=255)                                       aaaGAACCTTCAACATTAAAAATATTCATTTGTTTGAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAAt                                                      <  2:119146/100‑1 (MQ=255)                                       aaaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaaa                                                       >  2:1282747/1‑99 (MQ=255)                                       aaaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAAt                                                      >  2:809879/1‑100 (MQ=255)                                        aaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATa                                                     >  1:1188900/1‑100 (MQ=255)                                          gAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGt                                                   >  1:494466/1‑100 (MQ=255)                                            aCCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtat                                                 <  1:360336/100‑1 (MQ=255)                                                taaaCATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATAAAAAAAAAAAAAAATATCCAAAAATAGTATACtt                                             <  2:1469019/98‑1 (MQ=255)                                                   aCATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTAt                                          <  2:528634/100‑1 (MQ=255)                                                      tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                       <  1:817856/100‑1 (MQ=255)                                                      tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                       >  1:993736/1‑100 (MQ=255)                                                      tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                       >  1:605407/1‑100 (MQ=255)                                                       gAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                      <  1:429231/100‑1 (MQ=255)                                                         aaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCaga                                    <  1:1353289/100‑1 (MQ=255)                                                          aaaTATCCGTTTGTTTGAAAAAAAAAAGTAGTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAATAATATCCCAAAATAGTATACTTTATTCCataa                                   <  2:1258905/100‑1 (MQ=255)                                                           aaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAGATAGTATACTTTATTCCagaag                                  <  2:1387186/100‑1 (MQ=255)                                                             tatCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAagag                                <  2:318026/100‑1 (MQ=255)                                                              atCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGt                               >  2:96273/1‑100 (MQ=255)                                                              atCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGt                               <  1:1130502/100‑1 (MQ=255)                                                               tCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGtt                              >  1:851412/1‑100 (MQ=255)                                                                 cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                            >  1:736781/1‑100 (MQ=255)                                                                 cATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAACATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCa                            >  1:1104224/1‑100 (MQ=255)                                                                     tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                        >  1:63569/1‑100 (MQ=255)                                                                     tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTAAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                        <  2:36441/100‑1 (MQ=255)                                                                          gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                   >  1:1160400/1‑100 (MQ=255)                                                                          gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                   <  1:1335415/100‑1 (MQ=255)                                                                          gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                   >  1:1477591/1‑100 (MQ=255)                                                                          gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                   <  1:30218/100‑1 (MQ=255)                                                                           cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                  >  1:1530984/1‑100 (MQ=255)                                                                            aaaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGttt                 >  2:293600/1‑100 (MQ=255)                                                                             aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                <  2:1144273/100‑1 (MQ=255)                                                                             aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                <  1:709389/100‑1 (MQ=255)                                                                             aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                <  1:642828/100‑1 (MQ=255)                                                                               aaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGATCAATATAATGTTTGTc              >  1:450465/1‑100 (MQ=255)                                                                                aaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCt             <  1:150086/100‑1 (MQ=255)                                                                                   gATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCa          <  1:447903/100‑1 (MQ=255)                                                                                    attattAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCaa         >  1:257095/1‑100 (MQ=255)                                                                                    attattAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCaa         >  2:1264053/1‑100 (MQ=255)                                                                                     ttattaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAt        >  1:932725/1‑100 (MQ=255)                                                                                     ttattaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAt        >  2:461998/1‑100 (MQ=255)                                                                                      tattaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtt       >  1:219626/1‑100 (MQ=255)                                                                                       attaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttt      >  2:611535/1‑100 (MQ=255)                                                                                        ttaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAtttt     <  2:221487/100‑1 (MQ=255)                                                                                         taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt    >  1:586382/1‑100 (MQ=255)                                                                                          aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc   <  2:237459/100‑1 (MQ=255)                                                                                          aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc   <  1:255186/100‑1 (MQ=255)                                                                                          aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc   >  1:170992/1‑100 (MQ=255)                                                                                          aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc   <  2:1188893/100‑1 (MQ=255)                                                                                          aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc   <  1:757240/100‑1 (MQ=255)                                                                                          aggaagggAATTAATGCAATTATCGAAAATACAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc   <  2:1345558/100‑1 (MQ=255)                                                                                           ggaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCt  >  2:546943/1‑100 (MQ=255)                                                                                           ggaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCt  >  2:727179/1‑100 (MQ=255)                                                                                                 |                                                                                              CTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCT  >  NC_000913/1197643‑1197832

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