Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,937,168	G→A	*297* (TAG→TAA)	rbsB →	D‑ribose ABC transporter periplasmic binding protein, ribose chemotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,937,168	0	G	A	100.0%	175.4 / NA	54	*297* (TAG→TAA)	rbsB	D‑ribose ABC transporter periplasmic binding protein, ribose chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/28);  total (26/28)

CTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAAC  >  NC_000913/3937037‑3937302                                                                                                                                    |                                                                                                                                       cTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAAt                                                                                                                                 <  1:1755605/139‑1 (MQ=255)   agcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCa                                                                                                                               <  2:1739742/139‑1 (MQ=255)      agcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGt                                                                                                                            >  1:1442397/1‑139 (MQ=255)              tCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAAt                                                                                                                                 <  1:98787/126‑1 (MQ=255)              tCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAAt                                                                                                                                 >  2:98787/1‑126 (MQ=255)                     gCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATgg                                                                                                             <  1:310058/139‑1 (MQ=255)                      cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                            <  1:361205/139‑1 (MQ=255)                      cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                            <  2:1442397/139‑1 (MQ=255)                       tACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                           <  2:601333/139‑1 (MQ=255)                        aCCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                          >  2:192700/1‑139 (MQ=255)                        aCCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                          <  1:1343710/139‑1 (MQ=255)                        aCCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                          >  2:115953/1‑139 (MQ=255)                                  aTTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg                                                                                                <  1:115953/139‑1 (MQ=255)                                  aTTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg                                                                                                >  2:1558209/1‑139 (MQ=255)                                    tGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTc                                                                                              >  1:1757166/1‑139 (MQ=255)                                      gcgcGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg                                                                                                >  2:962341/1‑135 (MQ=255)                                      gcgcGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg                                                                                                <  1:962341/135‑1 (MQ=255)                                       cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAg                                                                                                                              <  2:694016/104‑1 (MQ=255)                                       cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAg                                                                                                                              >  1:694016/1‑104 (MQ=255)                                                cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  <  1:1558209/139‑1 (MQ=255)                                                  tCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  <  2:656440/137‑1 (MQ=255)                                                  tCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  >  1:656440/1‑137 (MQ=255)                                                     aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             <  2:1757166/139‑1 (MQ=255)                                                          gCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAata                                                                        <  1:192700/139‑1 (MQ=255)                                                                   gTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  <  2:1466166/120‑1 (MQ=255)                                                                   gTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  >  1:1466166/1‑120 (MQ=255)                                                                   gTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             <  2:367940/125‑1 (MQ=255)                                                                   gTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             >  1:367940/1‑125 (MQ=255)                                                                       tGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAg                                                           >  2:479434/1‑139 (MQ=255)                                                                        gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                            >  2:281482/1‑89 (MQ=255)                                                                        gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                            <  1:281482/89‑1 (MQ=255)                                                                        gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAgg                                                          >  1:526545/1‑139 (MQ=255)                                                                              cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCt                                                                                                                  <  2:1256215/77‑1 (MQ=255)                                                                              cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCt                                                                                                                  >  1:1256215/1‑77 (MQ=255)                                                                                           gCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  <  2:1396214/96‑1 (MQ=255)                                                                                           gCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  >  1:1396214/1‑96 (MQ=255)                                                                                            cTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACg                                      <  2:86408/139‑1 (MQ=255)                                                                                               aGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                           <  1:1718782/67‑1 (MQ=255)                                                                                               aGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                           >  2:1718782/1‑67 (MQ=255)                                                                                               aGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAata                                                                        >  2:1674101/1‑102 (MQ=255)                                                                                               aGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAata                                                                        <  1:1674101/102‑1 (MQ=255)                                                                                                     cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcgc                               >  2:1674716/1‑137 (MQ=255)                                                                                                     cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcgc                               <  1:1674716/137‑1 (MQ=255)                                                                                                     cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAt                             >  2:1497118/1‑139 (MQ=255)                                                                                                          gATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTAt                        >  2:35389/1‑139 (MQ=255)                                                                                                                  aCTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                          >  1:358277/1‑49 (MQ=255)                                                                                                                  aCTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                          <  2:358277/49‑1 (MQ=255)                                                                                                                      gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                          >  1:1231327/1‑45 (MQ=255)                                                                                                                      gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                          <  2:1231327/45‑1 (MQ=255)                                                                                                                          ttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  >  2:210378/1‑65 (MQ=255)                                                                                                                          ttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  <  1:210378/65‑1 (MQ=255)                                                                                                                          ttAAGCAGGAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATg        >  2:1193069/1‑139 (MQ=255)                                                                                                                             aGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             >  2:781977/1‑67 (MQ=255)                                                                                                                             aGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             <  1:781977/67‑1 (MQ=255)                                                                                                                                agtaatTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcgc                               <  1:961178/110‑1 (MQ=255)                                                                                                                                agtaatTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcgc                               >  2:961178/1‑110 (MQ=255)                                                                                                                                agtaatTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAAc  >  2:1445438/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       CTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAAC  >  NC_000913/3937037‑3937302

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