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%	156.3 / NA	49	*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 (21/28);  total (21/28)

CAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGA  >  NC_000913/3937033‑3937283                                                                                                                                        |                                                                                                                    cAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAAttt                                                                                                                  >  2:684309/1‑139 (MQ=255)   aaCTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTa                                                                                                                >  2:1032219/1‑139 (MQ=255)     cTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAAt                                                                                                              <  2:224374/139‑1 (MQ=255)     cTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAAt                                                                                                              <  1:585300/139‑1 (MQ=255)         cagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAgg                                                                                                          <  2:112683/139‑1 (MQ=255)              aCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAACGTTCAGGCTAAGTATCCGGTCGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAg                                                                                                           >  2:737786/1‑133 (MQ=255)              aCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTa                                                                                                     <  1:956408/139‑1 (MQ=255)              aCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATCTTAATCAg                                                                                                           <  1:737786/133‑1 (MQ=255)                             cccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATg                                                                                           >  2:366949/1‑134 (MQ=255)                             cccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATg                                                                                           <  1:366949/134‑1 (MQ=255)                              ccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAgg                                                                                                          <  1:413125/118‑1 (MQ=255)                              ccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAgg                                                                                                          >  2:413125/1‑118 (MQ=255)                               cGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACAt                                                                                    <  1:965077/139‑1 (MQ=255)                                       ttGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                         <  2:825/126‑1 (MQ=255)                                       ttGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                         >  1:825/1‑126 (MQ=255)                                           cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGAc                                                                                      <  1:591779/125‑1 (MQ=255)                                           cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGAc                                                                                      >  2:591779/1‑125 (MQ=255)                                                aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                         <  1:39005/117‑1 (MQ=255)                                                aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                         >  2:39005/1‑117 (MQ=255)                                                aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                         >  1:211727/1‑117 (MQ=255)                                                aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGt                                                                                         <  2:211727/117‑1 (MQ=255)                                                aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAg                                                                   >  1:66318/1‑139 (MQ=255)                                                    cGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGa                                                               <  1:774298/139‑1 (MQ=255)                                                     gTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                <  1:60102/121‑1 (MQ=255)                                                     gTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                >  2:60102/1‑121 (MQ=255)                                                         aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                          <  1:684309/139‑1 (MQ=255)                                                           aCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTa                                                        <  2:409170/139‑1 (MQ=255)                                                             cGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGa                                                                  <  2:707268/127‑1 (MQ=255)                                                             cGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGa                                                                  >  1:707268/1‑127 (MQ=255)                                                              gCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAata                                                     <  1:1032219/139‑1 (MQ=255)                                                                                  cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATc                                                                                                             <  1:354797/63‑1 (MQ=255)                                                                                  cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATc                                                                                                             >  2:354797/1‑63 (MQ=255)                                                                                  cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCg                                                                       <  2:829538/101‑1 (MQ=255)                                                                                  cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCg                                                                       >  1:829538/1‑101 (MQ=255)                                                                                  cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgc                                 >  1:809633/1‑139 (MQ=255)                                                                                    agaAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                               <  2:809633/139‑1 (MQ=255)                                                                                      aaaGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCa                             <  2:66318/139‑1 (MQ=255)                                                                                                cTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAg                                           >  2:484415/1‑115 (MQ=255)                                                                                                cTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAg                                           <  1:484415/115‑1 (MQ=255)                                                                                                         cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                          >  2:335705/1‑91 (MQ=255)                                                                                                         cGGTTGATCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                          <  1:335705/91‑1 (MQ=255)                                                                                                               aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg             >  2:548256/1‑130 (MQ=255)                                                                                                               aTCTGAAACTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg             <  1:548256/130‑1 (MQ=255)                                                                                                                      aCTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                          >  1:1129707/1‑78 (MQ=255)                                                                                                                      aCTGGTTGTTAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                          <  2:1129707/78‑1 (MQ=255)                                                                                                                          gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGGCATAAATACGTCATCGACAGATGAAc                                                             >  2:800230/1‑71 (MQ=255)                                                                                                                          gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAAc                                                             <  1:800230/71‑1 (MQ=255)                                                                                                                          gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGa  >  2:140325/1‑130 (MQ=255)                                                                                                                          gttgttAAGCAGTAATTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGa  <  1:140325/130‑1 (MQ=255)                                                                                                                                    agtaatTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg             <  1:319348/109‑1 (MQ=255)                                                                                                                                    agtaatTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg             >  2:319348/1‑109 (MQ=255)                                                                                                                                        |                                                                                                                    CAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGA  >  NC_000913/3937033‑3937283

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