Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,936,381	C→T	T35M (ACG→ATG)	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,936,381	0	C	T	100.0%	160.5 / NA	47	T35M (ACG→ATG)	rbsB	D‑ribose ABC transporter periplasmic binding protein, ribose chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/23);  total (24/23)

AGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCA  >  NC_000913/3936246‑3936504                                                                                                                                        |                                                                                                                            aGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCt                                                                                                                          >  1:61163/1‑139 (MQ=255) aGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCt                                                                                                                          >  1:459896/1‑139 (MQ=255)      taataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaaca                                                                                                                     >  1:1000556/1‑139 (MQ=255)      taataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaaca                                                                                                                     >  2:1086295/1‑139 (MQ=255)      taataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaaca                                                                                                                     >  2:1013209/1‑139 (MQ=255)      taataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaaca                                                                                                                     >  1:759309/1‑139 (MQ=255)        ataaCGACTACAGGACATCTTGAATATGAACATGAAATAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGCCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGGCTCCATGCTTaacagc                                                                                                                   <  1:315914/139‑3 (MQ=255)        ataaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                                                   <  1:1013209/139‑1 (MQ=255)          aaCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAAccc                                                                                                                 <  1:875419/139‑1 (MQ=255)           aCGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCg                                                                                                                <  2:985308/139‑1 (MQ=255)                    ggACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTa                                                                                                       <  1:1086295/139‑1 (MQ=255)                     gACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTAt                                                                                                      <  1:387134/139‑1 (MQ=255)                          cTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTg                                                                                                 >  1:928178/1‑139 (MQ=255)                                  tGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATgg                                                                                         <  2:1000556/139‑1 (MQ=255)                                   gAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgc                                                                                        >  1:873374/1‑139 (MQ=255)                                    aaCATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgcg                                                                                       >  2:1040052/1‑139 (MQ=255)                                     aCATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGgcgc                                                                                      >  2:600342/1‑139 (MQ=255)                                        tGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCa                                                                                   >  2:859430/1‑139 (MQ=255)                                               aCTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAagag                                                                            <  1:1018455/139‑1 (MQ=255)                                                   gCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAAccc                                                                                                                 >  1:216352/1‑98 (MQ=255)                                                   gCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAAccc                                                                                                                 <  2:216352/98‑1 (MQ=255)                                                        ccTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTAt                                                                                                      >  1:242374/1‑104 (MQ=255)                                                        ccTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTAt                                                                                                      <  2:242374/104‑1 (MQ=255)                                                                ccGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCta                                                           <  1:850309/139‑1 (MQ=255)                                                                 cGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCgg                                                                        >  1:767780/1‑125 (MQ=255)                                                                 cGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCgg                                                                        <  2:767780/125‑1 (MQ=255)                                                                      ttGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                  <  2:438888/126‑1 (MQ=255)                                                                      ttGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                  >  1:438888/1‑126 (MQ=255)                                                                          gcTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCtggtg                                                 <  1:466765/139‑1 (MQ=255)                                                                              aGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTc                                                                                                             >  1:850727/1‑75 (MQ=255)                                                                              aGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTc                                                                                                             <  2:850727/75‑1 (MQ=255)                                                                                  ccACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGa                                         <  1:844770/139‑1 (MQ=255)                                                                                   cACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAAccc                                                                                                                 <  2:817448/66‑1 (MQ=255)                                                                                   cACCGTCAGTGCGAATGCGATGGCAAAAGACAACATCGCGCTGGTGGTCTCCATGCTTAACAAccc                                                                                                                 >  1:817448/1‑66 (MQ=255)                                                                                    aCCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACt                                       <  2:459896/139‑1 (MQ=255)                                                                                                aaTGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                                                   <  2:322660/51‑1 (MQ=255)                                                                                                aaTGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTaacaac                                                                                                                   >  1:322660/1‑51 (MQ=255)                                                                                                aaTGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                  <  2:128169/100‑1 (MQ=255)                                                                                                aaTGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAc                                                                  >  1:128169/1‑100 (MQ=255)                                                                                                  tGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCg                         >  2:807915/1‑139 (MQ=255)                                                                                                      aTGGCAAAAGACACCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGa                     <  2:928178/139‑1 (MQ=255)                                                                                                                  acCATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGa         <  1:726518/139‑1 (MQ=255)                                                                                                                   ccATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGaa        >  2:887670/1‑139 (MQ=255)                                                                                                                    cATCGCGCTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAAc       >  2:559001/1‑139 (MQ=255)                                                                                                                         cgcTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACt                                                                 >  2:799138/1‑76 (MQ=255)                                                                                                                         cgcTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACt                                                                 <  1:799138/76‑1 (MQ=255)                                                                                                                         cgcTGGTGGTCTCCATGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGAGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCa  >  1:405628/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                            AGCAGTAATAACGACTACAGGACATCTTGAATATGAACATGAAAAAACTGGCTACCCTGGTTTCCGCTGTTGCGCTAAGCGCCACCGTCAGTGCGAATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCA  >  NC_000913/3936246‑3936504

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