Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,933,833	C→A	A19D (GCC→GAC)	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,833	0	C	A	100.0%	95.5 / NA	29	A19D (GCC→GAC)	rbsA	ribose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/11);  total (18/11)

GTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCTATACTCGCGATGCCGGTACGCTTTTATGG  >  NC_000913/3933705‑3933966                                                                                                                                 |                                                                                                                                      gTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGgcg                                                                                                                            <  1:409685/140‑1 (MQ=255)        gCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGtt                                                                                                                     <  2:297393/140‑1 (MQ=255)               agaATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTc                                                                                                              >  2:127900/1‑140 (MQ=255)                           gTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCcgc                                                                                                  >  1:413578/1‑140 (MQ=255)                                        aTTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCtgg                                                                                     <  2:194492/140‑1 (MQ=255)                                             ctGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGc                                                                                >  2:282236/1‑140 (MQ=255)                                               gtgCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGa                                                                              >  2:337191/1‑140 (MQ=255)                                               gtgCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGa                                                                              >  2:11178/1‑140 (MQ=255)                                                tgCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGaa                                                                             >  2:33395/1‑140 (MQ=255)                                                       gTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGgcg                                                                      <  1:29549/140‑1 (MQ=255)                                                          aCGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCgg                                                                   >  1:238339/1‑140 (MQ=255)                                                          aCGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCgg                                                                   >  2:14220/1‑140 (MQ=255)                                                          aCGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCgg                                                                   >  1:268572/1‑140 (MQ=255)                                                                        tCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCatg                                                     >  2:156378/1‑140 (MQ=255)                                                                          aTGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCatgat                                                   >  2:255570/1‑140 (MQ=255)                                                                                         cAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTgg                                    <  2:286359/140‑1 (MQ=255)                                                                                          aGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGc                                   <  1:255570/140‑1 (MQ=255)                                                                                            cTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCAt                                 >  1:205310/1‑140 (MQ=255)                                                                                            cTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCAt                                 <  1:179918/140‑1 (MQ=255)                                                                                              tAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCt                               >  2:127733/1‑140 (MQ=255)                                                                                                aaGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCtat                             >  1:188660/1‑140 (MQ=255)                                                                                                aaGGCATCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCtat                             >  2:407490/1‑140 (MQ=255)                                                                                                      tCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCTATACTcgc                       >  1:69487/1‑140 (MQ=255)                                                                                                      tCGATAAAGCCTTCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCTATACTcgc                       >  2:58625/1‑140 (MQ=255)                                                                                                                 ttCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCTATACTCGCGATGCCGGTAc            <  1:282236/140‑1 (MQ=255)                                                                                                                 ttCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCTATACTCGCGATGCCGGTAc            <  2:413578/140‑1 (MQ=255)                                                                                                                  tCCCGGGCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCTATACTCGCGATGCCGGTACg           >  2:30845/1‑140 (MQ=255)                                                                                                                        gCGTAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCTATACTCGCGATGCCGGTACGCTTTTa     <  1:407490/140‑1 (MQ=255)                                                                                                                           tAAAAGACCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCTATACTCGCGATGCCGGTACGCTTTTATgg  <  1:14220/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                      GTAATTCGCAGCGGAGAATGTTCTCCGTATGCGAATATCATTCTCTGTGCTGGCGTGACGTTCTGAGGCCGTCATGGAAGCATTACTTCAGCTTAAAGGCATCGATAAAGCCTTCCCGGGCGTAAAAGCCCTCTCGGGCGCAGCGTTAAATGTCTATCCGGGCCGCGTGATGGCGCTGGTGGGCGAAAACGGCGCGGGTAAATCCACCATGATGAAAGTGCTTACTGGCATCTATACTCGCGATGCCGGTACGCTTTTATGG  >  NC_000913/3933705‑3933966

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