Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,872,022 G→A 29.2% R157C (CGC→TGC)  ispD ← 2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,872,0220GA29.2% 38.3 / 18.8 34R157C (CGC→TGC) ispD2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base A (6/4);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAGGTCATCCTGATGCAAACA  >  NC_000913/2871882‑2872160
                                                                                                                                             |                                                                                                                                          
cTGAGGATGGAATCCGCAATAGTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGAt                                                                                                                                          <  1:420874/143‑1 (MQ=255)
     gATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTCCGATCAACg                                                                                                                                     >  1:91983/1‑143 (MQ=255)
          tctCCGCAATATGCCAGCGCCGAGGCTTCGCCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATg                                                                                                                                <  1:7761/141‑1 (MQ=255)
              cGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCa                                                                                                                            <  1:617410/143‑1 (MQ=255)
               gCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGACAGACAGTCATGTAGCAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGACGTTGCGTTCAACGGTGTGGGcga                                                                                                                           >  1:337453/1‑141 (MQ=255)
                 aaTATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGGGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATg                                                                                                                         <  2:91088/143‑1 (MQ=255)
                           cgcgGAGGTTTGCTCGGTAATAGTCGCGCCTGCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGc                                                                                                               <  1:644561/143‑1 (MQ=255)
                            gcggAGGCTTCGTCGGTAGTAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGcc                                                                                                              <  2:564289/140‑1 (MQ=255)
                                   cTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAACAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCCATGGCATTTTTGCCCGGTTCg                                                                                                       >  2:68804/1‑143 (MQ=255)
                                     tcgtcgGTAGTAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAGATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGc                                                                                                     <  2:270687/143‑1 (MQ=255)
                                        tcgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTGTTGCCCGGTTCGGAACg                                                                                                  >  2:51571/1‑143 (MQ=255)
                                        tcgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGAGCCAGCGCGTGACACCCAGCCCGGCAATCGACGGTATGAGCACTGGCATTTTTGCCCGGTTCGGCACg                                                                                                  >  1:123014/1‑143 (MQ=255)
                                          ggTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCAATCAACGGTATGAGCAATGGCATATTTGCCCGGTAAGGAACGtt                                                                                                >  2:48010/1‑143 (MQ=255)
                                             aatctgcgcgACTTGATTGAGAGCGCGCGTCAGCGAGTCATGTACCAGCTCACGAGGGAAACATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACCTTTCa                                                                                             <  2:546210/137‑1 (MQ=255)
                                             aaTAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTCCAATCAACGGTATGAGCAATGGCATTTTTGCCGGGTTCGGCACGTTTCa                                                                                             >  1:621966/1‑143 (MQ=255)
                                                  tCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGACCCATGGCATCTTTGACCGCTTCCCCACGGCTCCTAGTa                                                                                        >  2:256812/1‑143 (MQ=255)
                                                    gcgcCTTCATTGTGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAATATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATc                                                                                      <  2:448504/143‑1 (MQ=255)
                                                        cTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCgcgc                                                                                  >  1:337145/1‑143 (MQ=255)
                                                                         ggcAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGAt                                                                 <  2:169361/141‑1 (MQ=255)
                                                                           cAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATcc                                                               <  2:647114/143‑1 (MQ=255)
                                                                           cAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATcc                                                               >  2:513347/1‑143 (MQ=255)
                                                                                         aCAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGGGCGGCGTGGATCCCCCCCGTGCGGCTg                                                 >  1:508887/1‑143 (MQ=255)
                                                                                           aGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGACACGTTTCCTAGTATCGCGCACTCGTGCGGCGAGGATCCCCCCCGTGCGGCTggg                                               >  1:471998/1‑142 (MQ=255)
                                                                                           aGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGACACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTggg                                               >  1:40890/1‑142 (MQ=255)
                                                                                             cccACGAGGGAAAAATTGGGGCGTCAGCGCGTGCCATAAGCCGTTGCAAGCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGttt                                             <  1:220484/141‑1 (MQ=255)
                                                                                                aCGAGGGAAAGAGGGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTGTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGc                                          <  1:200259/143‑1 (MQ=255)
                                                                                                    gggAAAGATTGCGGCGTCAGCGCGTGCCATCAGCCGGTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCaa                                      <  2:417642/143‑1 (MQ=255)
                                                                                                       aaaaaTTGCGGCGTCAGCGGGTGCCATACGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGc                                   <  2:445433/143‑1 (MQ=255)
                                                                                                          aaTTGCGGCGTCAGCGCGTGCCATACGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaa                                <  2:183114/143‑1 (MQ=255)
                                                                                                           aTTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTCTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac                               >  2:150668/1‑143 (MQ=255)
                                                                                                                 gcgTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCCCTCAACGCCGACAAgcgc                         >  1:75598/1‑143 (MQ=255)
                                                                                                                      aGCGCGTTCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGGAGCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAg                    <  2:356673/143‑1 (MQ=255)
                                                                                                                                        gTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAGGTCATCCTGATGCAAAc   >  1:556983/1‑143 (MQ=255)
                                                                                                                                         ttGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAGGTCATCCTGATGCAAACa  >  1:35989/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                          
CTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCATAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAGGTCATCCTGATGCAAACA  >  NC_000913/2871882‑2872160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: