Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,933,826	Δ1 bp	coding (49/1506 nt)	rbsA →	ribose ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,933,826	0	G	.	100.0%	128.9 / NA	31	coding (49/1506 nt)	rbsA	ribose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (17/14);  total (17/14)

CTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAA  >  NC_000913/3933754‑3933892                                                                         |                                                                   cTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGccc                                                           <  2:783282/81‑1 (MQ=255) cTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGcc                                                            <  1:783281/80‑1 (MQ=255)    gCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCctct                                                        <  1:783283/81‑1 (MQ=255)    gCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCct                                                          <  2:783284/79‑1 (MQ=255)     cGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCtctc                                                       <  2:783286/81‑1 (MQ=255)     cGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCtctc                                                       <  2:783285/81‑1 (MQ=255)        gACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCggg                                                    >  2:783335/1‑81 (MQ=255)        gACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCggg                                                    >  1:783331/1‑81 (MQ=255)           gTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGgcgc                                                 <  1:783287/81‑1 (MQ=255)            ttCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGgcg                                                  >  1:783349/1‑79 (MQ=255)                gAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGt                                            >  1:783347/1‑81 (MQ=255)                    ccGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTaaa                                        >  1:783351/1‑81 (MQ=255)                    ccGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTaaa                                        <  1:783288/81‑1 (MQ=255)                     cGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAAt                                       >  2:783341/1‑81 (MQ=255)                       tCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGt                                     >  1:783340/1‑81 (MQ=255)                         aTGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTc                                    >  1:783348/1‑80 (MQ=255)                                cATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCggg                            >  1:783350/1‑81 (MQ=255)                                cATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCggg                            <  2:783289/81‑1 (MQ=255)                                   tACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCcg                         >  2:783366/1‑81 (MQ=255)                                     cTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCcgcg                       >  1:783333/1‑81 (MQ=255)                                        cAGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGa                    >  2:783355/1‑81 (MQ=255)                                         aGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGAt                   <  1:783290/81‑1 (MQ=255)                                         aGCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGAt                   <  2:783291/81‑1 (MQ=255)                                          gCTTAAAGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATg                  >  1:783327/1‑81 (MQ=255)                                              aaaGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGgcgc              >  2:783356/1‑81 (MQ=255)                                               aaGGCATCGATAAAGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCt             <  1:783292/81‑1 (MQ=255)                                                          aaaGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaa   >  1:783361/1‑80 (MQ=255)                                                          aaaGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaa   >  1:783357/1‑80 (MQ=255)                                                          aaaGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaa   >  1:783346/1‑80 (MQ=255)                                                           aaGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaaa  <  2:783293/80‑1 (MQ=255)                                                           aaGCCTTCCCGGGC‑TAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaaa  <  2:783294/80‑1 (MQ=255)                                                                         |                                                                   CTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAA  >  NC_000913/3933754‑3933892

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