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%	176.1 / 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 (18/36);  total (18/36)

CTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAAC  >  NC_000913/3937037‑3937302                                                                                                                                    |                                                                                                                                       cTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAAt                                                                                                                                 <  2:1016380/139‑1 (MQ=255)    gcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAg                                                                                                                              >  1:580502/1‑139 (MQ=255)          aCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                                        >  1:722701/1‑139 (MQ=255)          aCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                                        <  2:613593/139‑1 (MQ=255)          aCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                                        <  2:237896/139‑1 (MQ=255)              tCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGAc                                                                                                                    >  1:829752/1‑139 (MQ=255)                   cAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGAt                                                                                                               <  2:690078/139‑1 (MQ=255)                    aGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATg                                                                                                              <  2:722701/139‑1 (MQ=255)                      cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                            <  2:392848/139‑1 (MQ=255)                      cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                            <  2:132847/139‑1 (MQ=255)                         cccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGAc                                                                                                         <  1:734482/139‑1 (MQ=255)                                  aTTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg                                                                                                <  1:736137/139‑1 (MQ=255)                                    tGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGt                                                                                                                         >  2:746096/1‑112 (MQ=255)                                    tGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGt                                                                                                                         <  1:746096/112‑1 (MQ=255)                                       cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAAGTGGTTGTTAAGCAGTAATTTTAATCa                                                                                                                               >  2:755164/1‑103 (MQ=255)                                       cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCa                                                                                                                               <  1:755164/103‑1 (MQ=255)                                       cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                                   <  1:533922/131‑1 (MQ=255)                                       cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                                   >  2:533922/1‑131 (MQ=255)                                           aaaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACa                                                                                       >  2:849171/1‑139 (MQ=255)                                            aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAg                                                                                      >  1:661126/1‑139 (MQ=255)                                                cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  <  1:1033359/139‑1 (MQ=255)                                                cGTCGAAACCGCAGATAAAGTGCTGAAAGGAGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  <  2:380416/139‑1 (MQ=255)                                                  tCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAAc                                                                                <  1:361991/139‑1 (MQ=255)                                                     aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             <  1:303637/139‑1 (MQ=255)                                                     aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             <  2:661126/139‑1 (MQ=255)                                                     aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             <  2:925466/139‑1 (MQ=255)                                                     aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             <  1:145173/139‑1 (MQ=255)                                                     aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTGGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCGGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                             <  2:202201/139‑1 (MQ=255)                                                                        gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  <  1:383376/115‑1 (MQ=255)                                                                        gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                                                  >  2:383376/1‑115 (MQ=255)                                                                        gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAgg                                                          >  2:811961/1‑139 (MQ=255)                                                                           aGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCa                                                       <  1:553901/139‑1 (MQ=255)                                                                                agaAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                                                  <  2:953098/139‑1 (MQ=255)                                                                                agaAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                                                  <  2:385590/139‑1 (MQ=255)                                                                                          ggCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAaca                                        <  2:854619/139‑1 (MQ=255)                                                                                          ggCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAaca                                        <  1:113273/139‑1 (MQ=255)                                                                                            cTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                            >  1:487782/1‑85 (MQ=255)                                                                                            cTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAt                                                                                            <  2:487782/85‑1 (MQ=255)                                                                                                   tCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcgc                               <  2:38796/139‑1 (MQ=255)                                                                                                     cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                            >  1:599040/1‑60 (MQ=255)                                                                                                     cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                                            <  2:599040/60‑1 (MQ=255)                                                                                                           aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATg                       >  2:284122/1‑139 (MQ=255)                                                                                                           aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATg                       <  2:580502/139‑1 (MQ=255)                                                                                                           aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATg                       <  1:8936/139‑1 (MQ=255)                                                                                                              tGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGAc                    <  1:372299/139‑1 (MQ=255)                                                                                                              tGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGAc                    <  1:849171/139‑1 (MQ=255)                                                                                                                  aCTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAata                                                                        <  2:387164/83‑1 (MQ=255)                                                                                                                  aCTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAata                                                                        >  1:387164/1‑83 (MQ=255)                                                                                                                  aCTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATcc                >  1:458186/1‑139 (MQ=255)                                                                                                                      gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGaa            >  2:300293/1‑139 (MQ=255)                                                                                                                       ttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCa                              >  1:283858/1‑120 (MQ=255)                                                                                                                       ttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCa                              <  2:283858/120‑1 (MQ=255)                                                                                                                          ttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAat           >  1:19241/1‑136 (MQ=255)                                                                                                                          ttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAat           <  2:19241/136‑1 (MQ=255)                                                                                                                                agtaatTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAAc  >  2:448012/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       CTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAAC  >  NC_000913/3937037‑3937302

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