Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	160,782	C→T	*235* (TAG→TAA)	sfsA ←	sugar fermentation stimulation protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	160,782	0	C	T	100.0%	104.8 / NA	33	*235* (TAG→TAA)	sfsA	sugar fermentation stimulation protein A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/17);  total (16/17)

GCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGGCGCG  >  NC_000913/160661‑160915                                                                                                                          |                                                                                                                                      gCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAg                                                                                                                      >  1:384343/1‑139 (MQ=255)     taAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGAt                                                                                                                  <  1:824783/139‑1 (MQ=255)             aGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAag                                                                                                          >  2:1015453/1‑139 (MQ=255)             aGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAag                                                                                                          >  1:1198018/1‑139 (MQ=255)                 gCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTaa                                                                                                           <  2:951403/134‑1 (MQ=255)                 gCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTaa                                                                                                           >  1:951403/1‑134 (MQ=255)                                       tAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCg                                                                                <  2:384343/139‑1 (MQ=255)                                         aaCGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCt                                                                              >  1:687601/1‑139 (MQ=255)                                              aaCACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTa                                                                         >  2:207011/1‑139 (MQ=255)                                               acacTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTttatta                                                                                                                                     <  1:1029128/78‑1 (MQ=255)                                               acacTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTaa                                                                        >  1:101167/1‑139 (MQ=255)                                                    tGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAg                                                                   <  2:868250/139‑1 (MQ=255)                                                     gACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTAACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGa                                                                  <  1:872048/139‑1 (MQ=255)                                                      aCAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTa                                                                         <  2:1179045/131‑1 (MQ=255)                                                      aCAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTa                                                                         >  1:1179045/1‑131 (MQ=255)                                                             tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAaga                                                                                                         <  1:688946/92‑1 (MQ=255)                                                             tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAaga                                                                                                         >  2:688946/1‑92 (MQ=255)                                                             tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                                                      <  1:1143558/95‑1 (MQ=255)                                                             tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                                                      >  2:1143558/1‑95 (MQ=255)                                                             tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAc                                                          >  2:1049438/1‑139 (MQ=255)                                                                 aGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAccccc                                                      <  1:207011/139‑1 (MQ=255)                                                                  ggAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACa                                                                                                                                   <  2:217241/61‑1 (MQ=255)                                                                  ggAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACa                                                                                                                                   >  1:217241/1‑61 (MQ=255)                                                                       aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGt                                                                                                                     <  1:658001/70‑1 (MQ=255)                                                                       aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGt                                                                                                                     >  2:658001/1‑70 (MQ=255)                                                                       aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                                            >  2:1119621/1‑79 (MQ=255)                                                                       aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                                            <  1:1119621/79‑1 (MQ=255)                                                                          gTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGa                                             <  1:849168/139‑1 (MQ=255)                                                                           tATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAg                                            <  2:1198018/139‑1 (MQ=255)                                                                                   cACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGAc                                    <  2:687601/139‑1 (MQ=255)                                                                                                  tAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTAt                     <  1:994479/139‑1 (MQ=255)                                                                                                        aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCt                                       >  2:234724/1‑115 (MQ=255)                                                                                                        aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCt                                       <  1:234724/115‑1 (MQ=255)                                                                                                         aCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCAt              >  2:303358/1‑139 (MQ=255)                                                                                                                     ctttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGgcgcg  <  2:875080/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                      GCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGGCGCG  >  NC_000913/160661‑160915

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