Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,863,452 A→G 29.6% P373P (CCA→CCG hflX → predicted GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,863,4520AG29.6% 58.1 / 27.7 64P373P (CCA→CCGhflXpredicted GTPase
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (11/8);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

GCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGA  >  minE/2863389‑2863513
                                                               |                                                             
gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:899222/69‑1 (MQ=255)
gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:734596/69‑1 (MQ=255)
gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:514537/69‑1 (MQ=255)
gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:3166002/69‑1 (MQ=255)
gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:187785/69‑1 (MQ=255)
 cTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCATCTGCCACCGCa                                                          <  1:1055125/68‑1 (MQ=255)
  tttttCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAgg                                                        <  1:2205923/69‑1 (MQ=255)
  tttttCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAgg                                                        <  1:2172608/69‑1 (MQ=255)
  tttttCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAgg                                                        <  1:176796/69‑1 (MQ=255)
   ttttCAGGCTTTGACGGGGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAgg                                                        <  1:746617/68‑1 (MQ=255)
    tttCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:2385807/65‑1 (MQ=255)
          gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                >  1:2417473/1‑69 (MQ=255)
          gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                >  1:1862517/1‑69 (MQ=255)
          gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                >  1:1851809/1‑69 (MQ=255)
          gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                >  1:1351995/1‑69 (MQ=255)
                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:2690075/69‑1 (MQ=255)
                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:1172859/69‑1 (MQ=255)
                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:1955753/69‑1 (MQ=255)
                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:1345896/69‑1 (MQ=255)
                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:2317449/69‑1 (MQ=255)
                                gCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:2298353/37‑1 (MQ=255)
                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAg                                                         >  1:480932/1‑36 (MQ=255)
                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:161298/1‑54 (MQ=255)
                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:2728995/1‑54 (MQ=255)
                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:2166537/1‑54 (MQ=255)
                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:3131461/1‑54 (MQ=255)
                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:1370601/1‑54 (MQ=255)
                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAACGGCGTCTgaga                                          >  1:203482/1‑51 (MQ=255)
                                          gcagcaTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCaa                 <  1:1430466/68‑1 (MQ=255)
                                          gcagcaTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCaa                 <  1:2084463/68‑1 (MQ=255)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTa                               >  1:350162/1‑51 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:1355996/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:1628336/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:526767/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:2739706/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:971609/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:2691715/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:2437332/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:1715876/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:2281596/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:1043125/1‑52 (MQ=39)
                                             gcaTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGa             <  1:2625580/69‑1 (MQ=255)
                                              caTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGa             <  1:3191018/68‑1 (MQ=255)
                                              caTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGa             <  1:589439/68‑1 (MQ=255)
                                                tACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGaaaa          <  1:2019209/69‑1 (MQ=255)
                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGaaaaa         >  1:2894471/1‑68 (MQ=255)
                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:459042/1‑69 (MQ=255)
                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:757783/1‑69 (MQ=255)
                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:1312595/1‑69 (MQ=255)
                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:2931011/1‑69 (MQ=255)
                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:447038/1‑69 (MQ=255)
                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:3346106/1‑69 (MQ=255)
                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:3317135/1‑69 (MQ=255)
                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:318484/1‑69 (MQ=255)
                                                      gCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg    <  1:2848947/69‑1 (MQ=255)
                                                      gCGTCTGCCACAGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg    >  1:1431696/1‑69 (MQ=255)
                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:2181522/68‑1 (MQ=255)
                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:849140/68‑1 (MQ=255)
                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:803378/68‑1 (MQ=255)
                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:1943021/68‑1 (MQ=255)
                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:634910/68‑1 (MQ=255)
                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:2226142/68‑1 (MQ=255)
                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:3261932/68‑1 (MQ=255)
                                                         tCTGACGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:229727/68‑1 (MQ=255)
                                                               |                                                             
GCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGA  >  minE/2863389‑2863513

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