Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,262,619	C→T	G87R (GGA→AGA)	ispE ←	4‑diphosphocytidyl‑2‑C‑methylerythritol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,262,619	0	C	T	100.0%	175.4 / NA	45	G87R (GGA→AGA)	ispE	4‑diphosphocytidyl‑2‑C‑methylerythritol kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/25);  total (20/25)

GCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCCCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACGAAGCTCAATGCTGATGGTGTCGCC  >  NC_000913/1262472‑1262766                                                                                                                                                    |                                                                                                                                                    gCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGt                                                                                                                                                  >  2:806658/1‑151 (MQ=255)  ccccATTTCCCCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTAAATGCCCCCAGTCCCGTCGCGCCTTTTAATGAACCCCCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCg                                                                                                                                                   <  2:626728/149‑1 (MQ=255)   cccATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCtcgtcg                                                                                                                                                >  1:111130/1‑151 (MQ=255)      aTTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGa                                                                                                                                              >  1:222583/1‑150 (MQ=255)      aTTGCCGCCAGCTCACCCATGCTTAGCCCGTTTTGCCAGAGACGCTTTATTTCCACCCGACCCGTCGCGCCATTGGATGACCCGCCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTAGGa                                                                                                                                              <  2:629752/150‑1 (MQ=255)                 cccaaccATGCTTTGCCGGCATTGCCAGAGATTATTTAATGCCCCCAGGACCGCCGCGGCTTTGGATTAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTg                                                                                                                                  <  1:14582/146‑1 (MQ=255)                    aTCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCt                                                                                                                               >  1:714702/1‑151 (MQ=255)                      ccATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCt                                                                                                                               <  1:802714/149‑1 (MQ=255)                          gCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGc                                                                                                                          <  1:893664/150‑1 (MQ=255)                                 ccGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTtc                                                                                                                   >  1:30045/1‑150 (MQ=255)                                 ccGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGAGGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTtca                                                                                                                  <  2:89014/151‑1 (MQ=255)                                         cgagagATGATTTAATGCCACCAGGACCGTCGCGGTTTTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTgcatca                                                                                                               <  2:690946/144‑1 (MQ=255)                                         ccAGAGATGATTTAATGCCCCCAGGACAGTTGCGGCATTGGATGAACCAACGCCGAGACCGCCGCCCATCGGCAAAGGCTTGTCAATGCTGATATTCGCACAGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGGTTTCATCAACAAt                                                                                                          >  2:870981/1‑151 (MQ=255)                                           agagATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATcg                                                                                                        >  2:402284/1‑151 (MQ=255)                                           agagATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATc                                                                                                         <  1:923436/150‑1 (MQ=255)                                            gagaTGATTTAATGCCACCAGGACCGTCGCGGCTTTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATcaaca                                                                                                            <  1:985302/146‑1 (MQ=255)                                                tgatTTAATGCCACCAGGACCGTCGCGGCATTGGATGACCCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCt                                                                                                    <  1:407615/150‑1 (MQ=255)                                                   tttAATGCCCCCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTgc                                                                                                  <  1:148902/149‑1 (MQ=255)                                                   tttAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTgc                                                                                                  <  1:795460/149‑1 (MQ=255)                                                     tAATCCCACAAGATCGGTGGCGGCATTGGATGACCCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCGgcgc                                                                                                <  2:477116/149‑1 (MQ=255)                                                      aaTGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGa                                                                                              >  1:548192/1‑150 (MQ=255)                                                      aaTCCCACAGGGCCCGTCGCGGAATGGGAGGACCCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGaa                                                                                             <  2:803423/151‑1 (MQ=255)                                                          cccccaGGCCCGTCGCGGCATTGGATGACCCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGa                                                                                              <  2:611128/143‑1 (MQ=255)                                                             ccaGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCa                                                                                       <  1:240892/150‑1 (MQ=255)                                                               aGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGt                                                                                    >  1:273541/1‑151 (MQ=255)                                                                               tGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCa                                                                     <  1:101107/150‑1 (MQ=255)                                                                                    gAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCtt                                                               >  2:569334/1‑151 (MQ=255)                                                                                          ccaccgAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGACCAGGTTATCTTCATGTTCCACGCCTTCAACgg                                                         <  2:955604/148‑1 (MQ=255)                                                                                           cgccgAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACggg                                                        >  1:200815/1‑151 (MQ=255)                                                                                             ccgAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCg                                                      >  2:931596/1‑151 (MQ=255)                                                                                             ccgAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCg                                                      >  1:899272/1‑151 (MQ=255)                                                                                                agaCCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTa                                                   <  1:595475/151‑1 (MQ=255)                                                                                                  aCCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAAc                                                 >  2:103019/1‑151 (MQ=255)                                                                                                         cccATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGa                                           >  1:292854/1‑150 (MQ=255)                                                                                                              cGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATcc                                     >  1:342356/1‑151 (MQ=255)                                                                                                                 cAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAGGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCa                                  <  2:104959/151‑1 (MQ=255)                                                                                                                 cAAACGCTTGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCa                                  <  1:1018794/151‑1 (MQ=255)                                                                                                                         tGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCa                            <  1:565590/149‑1 (MQ=255)                                                                                                                         tGTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACg                          >  2:101740/1‑151 (MQ=255)                                                                                                                          gTCAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACg                          <  1:5978/150‑1 (MQ=255)                                                                                                                            cAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACGaa                        <  1:806704/150‑1 (MQ=255)                                                                                                                            cAATGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACGaa                        <  1:653022/150‑1 (MQ=255)                                                                                                                               tGCTGATATTCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACGAAGCTc                    >  1:539258/1‑151 (MQ=255)                                                                                                                                        tCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACGAAGCTCAATGCTGAt           >  2:152944/1‑151 (MQ=255)                                                                                                                                        tCGCACCGCTTCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACGAAGCTCAATGCTGAt           >  1:839588/1‑151 (MQ=255)                                                                                                                                                 ttCTCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACGAAGCTCAATGCTGATGGTGTCGcc  >  1:961317/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAAACGCTTGTCAATGCTGATATTCGCACCGCTTCCCGTCGGAAGACGCCCGCTGTCTGCCGCAGTTTTCATCAACAATCGCGCTGCGCGAACGATCAGGTTATCTTCATGTTCCACGCCTTCAACGGGCGTTAACAGACGAATATCCCCATCGTCACGAAGCTCAATGCTGATGGTGTCGCC  >  NC_000913/1262472‑1262766

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