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%	120.0 / NA	38	*235* (TAG→TAA)	sfsA	sugar fermentation stimulation protein A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/18);  total (20/18)

CTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGT  >  NC_000913/160662‑160909                                                                                                                         |                                                                                                                                cTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGt                                                                                                               >  1:891143/1‑139 (MQ=255)           cAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGAt                                                                                                            >  2:627755/1‑132 (MQ=255)           cAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGAt                                                                                                            <  1:627755/132‑1 (MQ=255)              aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGTGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTaa                                                                                                                       <  2:414179/118‑1 (MQ=255)              aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGTGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTaa                                                                                                                       >  1:414179/1‑118 (MQ=255)              aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGt                                                                                                               >  1:984907/1‑126 (MQ=255)              aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGt                                                                                                               <  2:984907/126‑1 (MQ=255)              aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGAt                                                                                                            <  2:289892/129‑1 (MQ=255)              aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGAt                                                                                                            >  1:289892/1‑129 (MQ=255)                gcctttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTGCTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                                                <  2:891143/136‑1 (MQ=255)                gCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                                                >  1:757843/1‑139 (MQ=255)                              cAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGa                                                                                  >  2:1192327/1‑139 (MQ=255)                                   aTCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAttt                                                                             <  1:1192327/139‑1 (MQ=255)                                        aaCGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTaa                                                                                                     >  1:654796/1‑110 (MQ=255)                                        aaCGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTaa                                                                                                     <  2:654796/110‑1 (MQ=255)                                              acacTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTttatta                                                                                                                               >  2:269609/1‑78 (MQ=255)                                                  ttGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCa                                                              >  1:517655/1‑139 (MQ=255)                                                         cTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTc                                                       >  2:699888/1‑139 (MQ=255)                                                            tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAc                                                    >  1:714724/1‑139 (MQ=255)                                                             gtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAcc                                                   <  2:757843/139‑1 (MQ=255)                                                                   aaaaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                                                <  2:987787/88‑1 (MQ=255)                                                                   aaaaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                                                >  1:987787/1‑88 (MQ=255)                                                                     aaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTg                                           >  1:1129969/1‑139 (MQ=255)                                                                      aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGAtttttt                                                                                                       >  2:1289950/1‑78 (MQ=255)                                                                      aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGAtttttt                                                                                                       <  1:1289950/78‑1 (MQ=255)                                                                            tttGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTg                                                                     <  2:1255429/106‑1 (MQ=255)                                                                            tttGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTg                                                                     >  1:1255429/1‑106 (MQ=255)                                                                            tttGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCt                                    <  2:1129969/139‑1 (MQ=255)                                                                             ttGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCtt                                   <  2:74550/139‑1 (MQ=255)                                                                             ttGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCtt                                   >  1:1132579/1‑139 (MQ=255)                                                                                 gcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGa                               <  2:517655/139‑1 (MQ=255)                                                                                         ccAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGtt                       <  2:1132579/139‑1 (MQ=255)                                                                                                           aGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAcccc                                                 >  1:797333/1‑95 (MQ=255)                                                                                                           aGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAcccc                                                 <  2:797333/95‑1 (MQ=255)                                                                                                             ttacttacTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAAttt                                                        >  2:585273/1‑86 (MQ=255)                                                                                                             ttacttacTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAAttt                                                        <  1:585273/86‑1 (MQ=255)                                                                                                             ttacttacTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAg                                      >  2:710137/1‑104 (MQ=255)                                                                                                             ttacttacTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAg                                      <  1:710137/104‑1 (MQ=255)                                                                                                              tacttactTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGAtgt  >  1:594838/1‑139 (MQ=255)                                                                                                                         |                                                                                                                                CTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGT  >  NC_000913/160662‑160909

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