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%	177.6 / 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 (23/31);  total (23/31)

CAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATG  >  NC_000913/3937033‑3937296                                                                                                                                        |                                                                                                                                 cAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAAttt                                                                                                                               >  2:248065/1‑139 (MQ=255) cAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAAttt                                                                                                                               >  1:1496917/1‑139 (MQ=255)    aCTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTaa                                                                                                                            <  1:1626572/139‑1 (MQ=255)     cTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAAt                                                                                                                           <  1:416192/139‑1 (MQ=255)     cTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAAt                                                                                                                           <  1:248065/139‑1 (MQ=255)        gcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAg                                                                                                                        <  2:1496917/139‑1 (MQ=255)        gcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAg                                                                                                                        <  1:1260673/139‑1 (MQ=255)                  tCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTG‑AAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGAcc                                                                                                             >  2:288045/1‑139 (MQ=255)                          cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                                  <  2:1690405/127‑1 (MQ=255)                          cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                                  >  1:1690405/1‑127 (MQ=255)                           tACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                     <  1:1323917/139‑1 (MQ=255)                            aCCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                    >  2:1571831/1‑139 (MQ=255)                                    aGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATa                                                                                            >  2:1453431/1‑139 (MQ=255)                                       ttGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGt                                                                                         >  2:1536769/1‑139 (MQ=255)                                          gcgcGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                      <  1:528049/139‑1 (MQ=255)                                               aaaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg                                                                                          <  2:1110900/130‑1 (MQ=255)                                               aaaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg                                                                                          >  1:1110900/1‑130 (MQ=255)                                                  ggCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGAt                                                                              <  1:1536769/139‑1 (MQ=255)                                                    cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                            <  1:1571831/139‑1 (MQ=255)                                                    cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAATTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                            <  2:899482/139‑1 (MQ=255)                                                        gAAACCGCAGATAAAGTGCTG‑AAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                       <  1:288045/139‑1 (MQ=255)                                                         aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGCTGATCTGAAACTGGTTGTTAAGCCGTAATTTTAATCAGGCTGTATGACCTCATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                       <  1:1100529/139‑1 (MQ=255)                                                              gCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAata                                                                  <  1:564332/139‑1 (MQ=255)                                                                 gATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTg                                                                                                           >  1:1216056/1‑95 (MQ=255)                                                                 gATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTg                                                                                                           <  2:1216056/95‑1 (MQ=255)                                                                      aGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                      >  2:954717/1‑111 (MQ=255)                                                                      aGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                      <  1:954717/111‑1 (MQ=255)                                                                            gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                      >  2:1137385/1‑105 (MQ=255)                                                                            gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                      <  1:1137385/105‑1 (MQ=255)                                                                                  cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGCCAGATGAACGTGTAATATAAAGAAAAGCATGGCAcgc                                              >  1:1667986/1‑139 (MQ=255)                                                                                  cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgc                                              <  2:1667986/139‑1 (MQ=255)                                                                                   gagaAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgcg                                             <  2:172043/139‑1 (MQ=255)                                                                                    agaAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                                            <  1:1453431/139‑1 (MQ=255)                                                                                    agaAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                                            <  2:230436/139‑1 (MQ=255)                                                                                                    gTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                      >  2:638786/1‑81 (MQ=255)                                                                                                    gTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                      <  1:638786/81‑1 (MQ=255)                                                                                                         cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTAt                                                                                                                 <  2:268492/49‑1 (MQ=255)                                                                                                         cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTAt                                                                                                                 >  1:268492/1‑49 (MQ=255)                                                                                                             tGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTa                   >  2:1711744/1‑139 (MQ=255)                                                                                                              gATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTAt                  >  1:922573/1‑139 (MQ=255)                                                                                                               aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                     <  2:1638059/55‑1 (MQ=255)                                                                                                               aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                     >  1:1638059/1‑55 (MQ=255)                                                                                                               aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATg                 >  1:1311527/1‑139 (MQ=255)                                                                                                                tCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATgg                <  2:922573/139‑1 (MQ=255)                                                                                                                 cTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATg                                                                             <  1:252217/77‑1 (MQ=255)                                                                                                                 cTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATg                                                                             >  2:252217/1‑77 (MQ=255)                                                                                                                      aCTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATcc          <  2:1247292/139‑1 (MQ=255)                                                                                                                          gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAggg                                                   <  1:858179/94‑1 (MQ=255)                                                                                                                          gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAggg                                                   >  2:858179/1‑94 (MQ=255)                                                                                                                              ttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                            <  1:1021934/65‑1 (MQ=255)                                                                                                                              ttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                            >  2:1021934/1‑65 (MQ=255)                                                                                                                              ttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATg  >  1:545570/1‑139 (MQ=255)                                                                                                                                 aGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCg                                                                                    >  2:874753/1‑54 (MQ=255)                                                                                                                                 aGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCg                                                                                    <  1:874753/54‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 CAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATG  >  NC_000913/3937033‑3937296

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: ‑