Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,655,057	T→C	T29A (ACG→GCG)	rspA ←	bifunctional D‑altronate/D‑mannonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,655,057	0	T	C	100.0%	140.7 / NA	40	T29A (ACG→GCG)	rspA	bifunctional D‑altronate/D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/20);  total (20/20)

GTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAA  >  NC_000913/1654959‑1655187                                                                                                   |                                                                                                                                   gTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTa                                                                                            >  1:1180579/1‑139 (MQ=255)  tGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTAc                                                                                           >  1:1328087/1‑139 (MQ=255)  tGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTAc                                                                                           >  2:622351/1‑139 (MQ=255)  tGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTAc                                                                                           >  1:202205/1‑139 (MQ=255)         cgcgACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCgg                                                                                    <  1:622351/139‑1 (MQ=255)                  tAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAAt                                                                                                        <  1:247391/110‑1 (MQ=255)                  tAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGCGATTTTTAAt                                                                                                        >  2:247391/1‑110 (MQ=255)                            aCAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCa                                                                 <  1:1194382/139‑1 (MQ=255)                                aGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCt                                                                                                                         >  2:595320/1‑79 (MQ=255)                                aGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCt                                                                                                                         <  1:595320/79‑1 (MQ=255)                                        cTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAg                                                                <  1:375208/128‑1 (MQ=255)                                        cTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAg                                                                >  2:375208/1‑128 (MQ=255)                                         tGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCt                                                    <  1:301417/139‑1 (MQ=255)                                          gCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCtt                                                   >  1:33803/1‑139 (MQ=255)                                           cAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTc                                                  >  2:167736/1‑139 (MQ=255)                                                 agagGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTc                                                                                                                           <  1:138268/60‑1 (MQ=255)                                                 agagGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTc                                                                                                                           >  2:138268/1‑60 (MQ=255)                                                     gCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGAtt                                        >  1:722302/1‑139 (MQ=255)                                                     gCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGAtt                                        >  1:302610/1‑139 (MQ=255)                                                         cGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGt                                                                                                                            >  2:510653/1‑51 (MQ=255)                                                         cGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGt                                                                                                                            <  1:510653/51‑1 (MQ=255)                                                                cTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGc                             <  1:167736/139‑1 (MQ=255)                                                                cTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGc                             >  2:731726/1‑139 (MQ=255)                                                                         aTTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGa                    <  2:1060722/139‑1 (MQ=255)                                                                          ttGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                             >  1:200419/1‑65 (MQ=255)                                                                          ttGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAAtt                                                                                             <  2:200419/65‑1 (MQ=255)                                                                          ttGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGt                                                                              <  1:266705/80‑1 (MQ=255)                                                                          ttGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGt                                                                              >  2:266705/1‑80 (MQ=255)                                                                               ggTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGa                    >  2:1130952/1‑133 (MQ=255)                                                                               ggTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGa                    <  1:1130952/133‑1 (MQ=255)                                                                                 tGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGaa            <  1:731726/139‑1 (MQ=255)                                                                                   gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTGCGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAAt          >  1:1110366/1‑139 (MQ=255)                                                                                   gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGc                             >  2:889890/1‑120 (MQ=255)                                                                                   gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGc                             <  1:889890/120‑1 (MQ=255)                                                                                        cccAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACg     <  2:1328087/139‑1 (MQ=255)                                                                                         ccAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGc    >  2:246653/1‑139 (MQ=255)                                                                                          cgaagCCCGCAATACCGTCCTCAGTGGTGATTTTTAACGTGACGAAATTATGCCCCGGACAGGTAACAAAAACTTGAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCa   <  2:837288/135‑1 (MQ=255)                                                                                           aaGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGTATGAAATCTACGCaa  <  2:899755/139‑1 (MQ=255)                                                                                           aaGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCaa  <  2:722302/139‑1 (MQ=255)                                                                                           aaGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCaa  <  2:117047/139‑1 (MQ=255)                                                                                                   |                                                                                                                                   GTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAA  >  NC_000913/1654959‑1655187

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