Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,871,308 G→A T158I (ACA→ATA)  ispF ← 2C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,871,3080GA83.3% 49.1 / 5.6 24T158I (ACA→ATA) ispF2C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (11/9);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTTAACATCATCCATATGGCAGCCGAGATCTTC  >  NC_000913/2871172‑2871429
                                                                                                                                        |                                                                                                                         
gAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTttatt                                                                                                                         >  2:457658/1‑139 (MQ=255)
   tatGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGcc                                                                                                                      <  1:457658/139‑1 (MQ=255)
      gCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTa                                                                                                                   >  1:552300/1‑139 (MQ=255)
                                cccAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGc                                                                                         <  1:605346/139‑1 (MQ=255)
                                  cAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATg                                                                                                                >  2:331395/1‑114 (MQ=255)
                                  cAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATg                                                                                                                <  1:331395/114‑1 (MQ=255)
                                                            ccGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATccc                                                             <  2:156811/139‑1 (MQ=255)
                                                                    gCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTc                                                     >  2:479057/1‑139 (MQ=255)
                                                                               cccGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTgg                                          >  1:436604/1‑139 (MQ=255)
                                                                                         tGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCgtagta                                              >  1:540257/1‑125 (MQ=255)
                                                                                         tGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCgtagta                                              <  2:540257/125‑1 (MQ=255)
                                                                                              tttACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTTAAcat                           <  2:552300/139‑1 (MQ=255)
                                                                                                        aGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTTAAcatc                          <  2:53432/130‑1 (MQ=255)
                                                                                                        aGGTAAGTGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTTAAcatc                          >  1:53432/1‑130 (MQ=255)
                                                                                                               tGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGt                                                                    >  2:201146/1‑81 (MQ=255)
                                                                                                               tGAGATTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGt                                                                    <  1:201146/81‑1 (MQ=255)
                                                                                                                  gaTTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACg                                  <  2:76529/112‑1 (MQ=255)
                                                                                                                  gaTTATCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACg                                  >  1:76529/1‑112 (MQ=255)
                                                                                                                      aTCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACGGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTTAACATCATCCATATGGCAGCCGAGATCtt   >  2:201987/1‑139 (MQ=255)
                                                                                                                      aTCAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTTAACATCATCCATATGGCAGCCGAGATCtt   >  2:510626/1‑139 (MQ=255)
                                                                                                                       tCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTTAACATCATCCATATGGCAGCCGAGATCTTc  >  1:468349/1‑139 (MQ=255)
                                                                                                                        cAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATccc                                                                                   >  2:438428/1‑57 (MQ=255)
                                                                                                                        cAAACTCAATCATTTTATTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATccc                                                                                   <  1:438428/57‑1 (MQ=255)
                                                                                                                                 tCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTcc                                                    >  2:256839/1‑79 (MQ=255)
                                                                                                                                 tCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTcc                                                    <  1:256839/79‑1 (MQ=255)
                                                                                                                                        |                                                                                                                         
GAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTTAACATCATCCATATGGCAGCCGAGATCTTC  >  NC_000913/2871172‑2871429

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