Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,669,714 T→G 13.6% intergenic (+38/+37) prpD → / ← acnB 2‑methylcitrate dehydratase/aconitate hydratase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,669,7140TG13.6% 136.7 / 17.6 59intergenic (+38/+37)prpD/acnB2‑methylcitrate dehydratase/aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base T (26/25);  new base G (4/4);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

ACCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTC  >  NC_002947/2669569‑2669862
                                                                                                                                                 |                                                                                                                                                    
aCCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCgat                                                                                                                                                  >  2:309955/1‑148 (MQ=255)
 cccGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCGCAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTg                                                                                                                                                 >  1:1102803/1‑150 (MQ=255)
  ccGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGc                                                                                                                                                >  2:174447/1‑150 (MQ=255)
    gCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCg                                                                                                                                              >  2:299367/1‑150 (MQ=255)
        tGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCGGGGc                                                                                                                                          <  2:113383/150‑1 (MQ=255)
          cccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCtt                                                                                                                                        >  1:15375/1‑150 (MQ=255)
           ccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCttt                                                                                                                                       <  2:1032581/150‑1 (MQ=255)
           ccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCGGGGCttt                                                                                                                                       <  1:309955/150‑1 (MQ=255)
                         agcagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTa                                                                                                                         >  2:1109899/1‑150 (MQ=255)
                                tttCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGc                                                                                                                  <  1:1000350/150‑1 (MQ=255)
                                    gAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTAc                                                                                                              <  1:174447/150‑1 (MQ=255)
                                     aaCTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTCCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGACTTTGTTTACAACCTGTAACTCAGCTTACg                                                                                                             >  1:1569678/1‑150 (MQ=255)
                                      aCTGTGCAAGCACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGc                                                                                                            <  1:244296/150‑1 (MQ=255)
                                        tgtgCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCt                                                                                                          >  1:914205/1‑150 (MQ=255)
                                                         aaGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGGTTGTGATATGAAGCGGCCGGCACTTTCCCTGAGATCAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTg                                                                                         <  1:283109/150‑1 (MQ=255)
                                                              ggAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGtt                                                                                    <  2:1068766/150‑1 (MQ=255)
                                                                  ggCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCt                                                                                <  2:491180/150‑1 (MQ=255)
                                                                  ggCACGGCTGTGCAAAGGTTCGTGGACTTGTTTGTGGTATCAGGCGTCCGGCACATGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCt                                                                                <  2:15375/150‑1 (MQ=255)
                                                                   gCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTg                                                                               <  1:819193/150‑1 (MQ=255)
                                                                   gCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCgggg                                                                                                                                           <  2:650671/90‑1 (MQ=255)
                                                                   gCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCgggg                                                                                                                                           >  1:650671/1‑90 (MQ=255)
                                                                       ggCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGc                                                                           >  1:524054/1‑150 (MQ=255)
                                                                               acaGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCggg                                                                                                                                            <  2:437709/77‑1 (MQ=255)
                                                                               acaGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCggg                                                                                                                                            >  1:437709/1‑77 (MQ=255)
                                                                                   gTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTc                                                               >  2:143789/1‑150 (MQ=255)
                                                                                    ttCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCg                                                              <  2:1102803/150‑1 (MQ=255)
                                                                                     tCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCTTTGGCTGCCGCTTCGCGGAACTCgg                                                             <  1:299367/150‑1 (MQ=255)
                                                                                       gTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTc                                                                        >  1:102298/1‑137 (MQ=255)
                                                                                       gTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTc                                                                        <  2:102298/137‑1 (MQ=255)
                                                                                          gACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATc                                                        >  2:150147/1‑150 (MQ=255)
                                                                                           aCTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTg                                                                                   <  1:1447116/122‑1 (MQ=255)
                                                                                           aCTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTg                                                                                   >  2:1447116/1‑122 (MQ=255)
                                                                                            cTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGAt                                                         >  1:1561837/1‑147 (MQ=255)
                                                                                            cTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGAt                                                         <  2:1561837/147‑1 (MQ=255)
                                                                                             ttgtttgtTATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTgg                                                     >  1:1593403/1‑150 (MQ=255)
                                                                                             ttgtttgtGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTgg                                                     >  2:515879/1‑150 (MQ=255)
                                                                                             ttgtttgtGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGACTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTgg                                                     <  2:1569678/150‑1 (MQ=255)
                                                                                                  tgtgATATGAAGCGACCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCgg                                                                                                                                             <  1:795503/57‑1 (MQ=255)
                                                                                                  tgtgATATGAAGCGACCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCgg                                                                                                                                             >  2:795503/1‑57 (MQ=255)
                                                                                                      atatGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTc                                            >  1:922394/1‑150 (MQ=255)
                                                                                                        atGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAg                                          <  2:914205/150‑1 (MQ=255)
                                                                                                             gCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAaccacc                                                                                                  <  2:89338/89‑1 (MQ=255)
                                                                                                             gCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAaccacc                                                                                                  >  1:89338/1‑89 (MQ=255)
                                                                                                              cggccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCg                                    <  2:524054/150‑1 (MQ=255)
                                                                                                                gccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTaa                                                                                                                        >  2:1604767/1‑64 (MQ=255)
                                                                                                                gccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTaa                                                                                                                        <  1:1604767/64‑1 (MQ=255)
                                                                                                                      acaGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGt                            >  1:410279/1‑150 (MQ=255)
                                                                                                                        aGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCg                          <  1:1109899/150‑1 (MQ=255)
                                                                                                                         ggCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCgg                         <  1:143789/150‑1 (MQ=255)
                                                                                                                            ccAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGt                                                    >  1:79357/1‑120 (MQ=255)
                                                                                                                            ccAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGt                                                    <  2:79357/120‑1 (MQ=255)
                                                                                                                                  aacaaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGc                >  2:1403633/1‑150 (MQ=255)
                                                                                                                                  aacaaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGAAGGTCACCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGc                >  1:1252684/1‑150 (MQ=255)
                                                                                                                                    caaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCg                                    >  1:1386440/1‑128 (MQ=255)
                                                                                                                                    caaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCg                                    <  2:1386440/128‑1 (MQ=255)
                                                                                                                                        gCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGa          >  1:976363/1‑150 (MQ=255)
                                                                                                                                         ccccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAg                                               <  1:1160090/112‑1 (MQ=255)
                                                                                                                                         ccccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAg                                               >  2:1160090/1‑112 (MQ=255)
                                                                                                                                                aTCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGCTCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGTTCTATGCTCAGGTAGCGGTACACGTCGGCAGCCATGCTGCCGATGTCTTTc  >  2:754771/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTC  >  NC_002947/2669569‑2669862

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