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

CCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGGCGCG  >  NC_000913/160659‑160915                                                                                                                            |                                                                                                                                      ccGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTGACCAGTTACTTACTTTATTACAATGTAACCGGc                                                                                                                        <  2:822095/139‑1 (MQ=255)                 aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAagag                                                                                                        >  1:436170/1‑139 (MQ=255)                          cGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCtt                                                                                               <  2:349996/139‑1 (MQ=255)                               aGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagca                                                                                          <  1:1049413/139‑1 (MQ=255)                                           aaCGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCt                                                                              >  2:261193/1‑139 (MQ=255)                                                 acacTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTaa                                                                                                           <  2:42607/104‑1 (MQ=255)                                                 acacTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTaa                                                                                                           >  1:42607/1‑104 (MQ=255)                                                   acTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGc                                                                      <  1:770424/139‑1 (MQ=255)                                                      tGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAg                                                                   <  2:912158/139‑1 (MQ=255)                                                        aCAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGAt                                                                                                                  <  2:1394150/90‑1 (MQ=255)                                                        aCAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGAt                                                                                                                  >  1:1394150/1‑90 (MQ=255)                                                               tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAc                                                          >  1:164088/1‑139 (MQ=255)                                                                gtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAcc                                                         <  1:1495233/139‑1 (MQ=255)                                                                   aGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAccccc                                                      <  1:261193/139‑1 (MQ=255)                                                                     gAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCct                                                    <  2:36447/139‑1 (MQ=255)                                                                      aaaaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCctc                                                   <  2:530085/139‑1 (MQ=255)                                                                             tATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                                            >  2:1249694/1‑75 (MQ=255)                                                                             tATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                                            <  1:1249694/75‑1 (MQ=255)                                                                                    gcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGa                                     >  1:395652/1‑139 (MQ=255)                                                                                     cACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCtt                                         >  2:212757/1‑134 (MQ=255)                                                                                     cACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCtt                                         <  1:212757/134‑1 (MQ=255)                                                                                     cACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGAc                                    >  2:519150/1‑139 (MQ=255)                                                                                     cACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGAc                                    <  2:395652/139‑1 (MQ=255)                                                                                                aaTGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTg                            >  2:883267/1‑136 (MQ=255)                                                                                                aaTGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTg                            <  1:883267/136‑1 (MQ=255)                                                                                                      aaTTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTAttt                   >  2:409508/1‑139 (MQ=255)                                                                                                        ctaactcactcgttactttattACCATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCATAATTTCCACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTtct                 <  1:133917/127‑1 (MQ=255)                                                                                                          aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                                            >  1:712525/1‑46 (MQ=255)                                                                                                          aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                                            <  2:712525/46‑1 (MQ=255)                                                                                                          aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTc                                                                                              <  1:663680/60‑1 (MQ=255)                                                                                                          aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTc                                                                                              >  2:663680/1‑60 (MQ=255)                                                                                                          aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTa                                                           >  1:1446384/1‑95 (MQ=255)                                                                                                          aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTa                                                           <  2:1446384/95‑1 (MQ=255)                                                                                                          aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCa               >  1:1418820/1‑139 (MQ=255)                                                                                                          aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCa               >  2:308198/1‑139 (MQ=255)                                                                                                          aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCa               >  2:521899/1‑139 (MQ=255)                                                                                                           agtaactggttactttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCAt              >  2:133917/10‑139 (MQ=255)                                                                                                            ccAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAc                                                          >  1:545635/1‑94 (MQ=255)                                                                                                            ccAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAc                                                          <  2:545635/94‑1 (MQ=255)                                                                                                                ttacttacTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAccccc                                                      <  2:1212351/94‑1 (MQ=255)                                                                                                                ttacttacTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAccccc                                                      >  1:1212351/1‑94 (MQ=255)                                                                                                                ttacttacTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGAtg         >  1:1259573/1‑139 (MQ=255)                                                                                                                     tactttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGgcg    <  1:521899/139‑1 (MQ=255)                                                                                                                       ctttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGgcgcg  >  2:1541994/1‑139 (MQ=255)                                                                                                                       ctttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGgcgcg  <  2:436170/139‑1 (MQ=255)                                                                                                                       ctttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGgcgcg  <  1:308198/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                      CCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGGCGCG  >  NC_000913/160659‑160915

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