Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,262,453 A→G V142A (GTT→GCT)  ispE ← 4‑(cytidine 5'‑diphospho)‑2‑C‑methyl‑D‑erythritol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,262,4530AG96.9% 104.3 / ‑4.1 32V142A (GTT→GCT) ispE4‑(cytidine 5'‑diphospho)‑2‑C‑methyl‑D‑erythritol kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (22/9);  total (23/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

CCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAACATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAA  >  NC_000913/1262331‑1262586
                                                                                                                          |                                                                                                                                     
ccGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCGCCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTc                                                                                                                      <  1:597315/140‑1 (MQ=255)
                    ccAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTc                                                                                                  >  1:542123/1‑140 (MQ=255)
                       gggTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATc                                                                                               <  2:97643/140‑1 (MQ=255)
                       gggTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATc                                                                                               >  1:357579/1‑140 (MQ=255)
                          tGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCGGCACCATTTCCGCCAGCTCATCCAt                                                                                            >  1:220182/1‑140 (MQ=255)
                            cgcCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTTCGCGCAGTGTCAGCCCCATTTCCGCCAGCTCATCCATGc                                                                                          >  2:297120/1‑140 (MQ=255)
                            cgcCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGc                                                                                          >  2:119295/1‑140 (MQ=255)
                             gcCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCt                                                                                         >  1:39003/1‑140 (MQ=255)
                             gcCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAACATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCt                                                                                         >  1:502356/1‑140 (MQ=255)
                                     acaCCGCTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCATCGGCATGCGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGCCAGCCCCATTTCCTCCAGCTCATCCATGCTTAGCCCGc                                                                                 <  2:546072/138‑1 (MQ=255)
                                          ccTTCTCTGGCGGATCCCCCGGAGTTAGTATTTCACCAACGCCTTCGGTAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGc                                                                            <  2:211019/139‑1 (MQ=255)
                                           cTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGcc                                                                           <  1:119295/140‑1 (MQ=255)
                                           cTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGcc                                                                           <  2:64511/140‑1 (MQ=255)
                                           cTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGcc                                                                           >  2:535094/1‑140 (MQ=255)
                                           cTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGcc                                                                           >  1:574866/1‑140 (MQ=255)
                                              ctctGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCTTCAGCCCCATTTTCGCGAGCTCATCGATGCGTAGCCCGCATTGCCaga                                                                        >  2:179545/1‑140 (MQ=255)
                                              ctctGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCaga                                                                        >  2:345087/1‑140 (MQ=255)
                                                      gATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCCTGCTTAGCCCGCACTGCCAGAGATGAttt                                                                >  1:234682/1‑140 (MQ=255)
                                                         ccACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAAt                                                             <  2:357579/140‑1 (MQ=255)
                                                             cGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGcca                                                         >  1:382854/1‑140 (MQ=255)
                                                             cGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGcca                                                         >  1:470575/1‑140 (MQ=255)
                                                                    aGTATTTCACCAACGTCTTCGGAAAGCGCAGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGAc                                                  <  1:323286/140‑1 (MQ=255)
                                                                              cAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCa                                        <  1:345087/140‑1 (MQ=255)
                                                                                aCGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCAtt                                      >  2:529594/1‑140 (MQ=255)
                                                                                aCGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCAtt                                      >  1:655364/1‑140 (MQ=255)
                                                                                 cGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTg                                     >  2:456983/1‑140 (MQ=255)
                                                                                       cGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGaa                               >  1:305249/1‑140 (MQ=255)
                                                                                       cGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCA‑TTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCa                                        >  1:89902/1‑130 (MQ=255)
                                                                                       cGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCA‑TTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCa                                        <  2:89902/130‑1 (MQ=255)
                                                                                                cgGCATGCCCCCGAACAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCAccgccg                      <  2:655364/140‑1 (MQ=255)
                                                                                                    aTGCCCCCGAAAAAAGGCAGGACCATCTGCCGCAAGCGTGAGCCCCATTCCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCCCCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGAc                  <  2:305249/140‑1 (MQ=255)
                                                                                                             aaCAAAGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCa         >  2:392634/1‑140 (MQ=255)
                                                                                                                  aGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGc    >  1:547086/1‑140 (MQ=255)
                                                                                                                  aGACAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGc    >  1:506654/1‑140 (MQ=255)
                                                                                                                    aCAGGAGCATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCaa  >  1:625622/1‑140 (MQ=255)
                                                                                                                          |                                                                                                                                     
CCGGAGTCGGAATACTTACACCAGGGTGCGCCACCAGATACCACTTCTCTGGCGGATCCACCGGCGTTAGTATTTCACCAACGCCTTCGGCAAACGCGGCATGCCCCCGAACAAAGACAGGAACATCTGCGCCCAGCGTCAGCCCCATTTCCGCCAGCTCATCCATGCTTAGCCCGCATTGCCAGAGATGATTTAATGCCACCAGGACCGTCGCGGCATTGGATGAACCACCGCCGAGACCGCCGCCCATCGGCAA  >  NC_000913/1262331‑1262586

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