Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,733,223	G→A	D317D (GAC→GAT)	clpB ←	protein disaggregation chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,733,223	0	G	A	100.0%	99.9 / NA	30	D317D (GAC→GAT)	clpB	protein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/15);  total (15/15)

GGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATGGTATGTAATTCGTCGATAAATAGGATGACG  >  NC_000913/2733090‑2733358                                                                                                                                      |                                                                                                                                        ggtggtGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAg                                                                                                                                    >  2:608624/1‑139 (MQ=255)      tgCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTGATCAAGcgtcg                                                                                                                               >  2:399217/1‑139 (MQ=255)       gCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGcgtcgt                                                                                                                              <  1:608624/139‑1 (MQ=255)           ttCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGTTATTCATCAAGCGTCGTGGCa                                                                                                                          >  2:547269/1‑139 (MQ=255)            tCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCAc                                                                                                                         <  1:304627/139‑1 (MQ=255)            tCGTAACGTTCTTGCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACTACGTCCCAGCGCAGCATCTCTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCAc                                                                                                                         <  1:406342/139‑1 (MQ=255)                aaCGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTAc                                                                                                                     >  1:399603/1‑139 (MQ=255)                   gTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCa                                                                                                                  <  2:52596/139‑1 (MQ=255)                             gCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACg                                                                                                                    <  2:602437/127‑1 (MQ=255)                             gCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACg                                                                                                                    >  1:602437/1‑127 (MQ=255)                                 cGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac                                                                                                    >  2:327757/1‑139 (MQ=255)                                                gTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGtt                                                                                     >  2:532387/1‑139 (MQ=255)                                                      tCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCa                                                                               >  2:456937/1‑139 (MQ=255)                                                               ggCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCgg                                                                      >  2:354555/1‑139 (MQ=255)                                                                                      gAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTa                                               <  1:547269/139‑1 (MQ=255)                                                                                             cgTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCgcgc                                        <  1:399217/139‑1 (MQ=255)                                                                                                tCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGa                                     <  2:554815/139‑1 (MQ=255)                                                                                                tCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGa                                     <  1:327757/139‑1 (MQ=255)                                                                                                       cagcaTCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATg                                >  1:434707/1‑137 (MQ=255)                                                                                                       cagcaTCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATg                                <  2:434707/137‑1 (MQ=255)                                                                                                          caTCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATGGTATg                           >  2:252937/1‑139 (MQ=255)                                                                                                           aTCTTTTTCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATGGTATGt                          <  1:252937/139‑1 (MQ=255)                                                                                                                 ttCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTg                                                             <  2:195428/98‑1 (MQ=255)                                                                                                                 ttCAATGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTg                                                             >  1:195428/1‑98 (MQ=255)                                                                                                                      tGTACTGGCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATGGTATGTAATTCGTCGAt               >  2:316846/1‑139 (MQ=255)                                                                                                                       gTACTGGCGATGTTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATGGTATGTAATTCGTCGATa              <  2:399603/139‑1 (MQ=255)                                                                                                                             gCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCa                                                                               >  2:296164/1‑68 (MQ=255)                                                                                                                             gCGATATTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCa                                                                               <  1:296164/68‑1 (MQ=255)                                                                                                                                  atTCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATGGTATGTAATTCGTCGATAAATAGGATGAc   >  2:605577/1‑139 (MQ=255)                                                                                                                                   ttCATCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATGGTATGTAATTCGTCGATAAATAGGATGACg  <  1:605577/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        GGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATTGCGCCATCGGCTTTACCCGCGCCGACCATGGTATGTAATTCGTCGATAAATAGGATGACG  >  NC_000913/2733090‑2733358

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