Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,863,467 T→C 26.3% R378R (CGT→CGC hflX → predicted GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,863,4670TC26.3% 77.8 / 27.2 72R378R (CGT→CGChflXpredicted GTPase
Reads supporting (aligned to +/- strand):  ref base T (23/30);  new base C (11/8);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

GCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATGGAGGAGGACGGCAGCGTAAG  >  minE/2863399‑2863535
                                                                    |                                                                    
gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                                      >  1:1851809/1‑69 (MQ=255)
gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                                      >  1:1862517/1‑69 (MQ=255)
gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                                      >  1:1351995/1‑69 (MQ=255)
gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                                      >  1:2417473/1‑69 (MQ=255)
         gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             <  1:1172859/69‑1 (MQ=255)
         gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             <  1:1955753/69‑1 (MQ=255)
         gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             <  1:2317449/69‑1 (MQ=255)
         gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             <  1:2690075/69‑1 (MQ=255)
         gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             <  1:1345896/69‑1 (MQ=255)
                        gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             >  1:1370601/1‑54 (MQ=255)
                        gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             >  1:2166537/1‑54 (MQ=255)
                        gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             >  1:161298/1‑54 (MQ=255)
                        gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             >  1:2728995/1‑54 (MQ=255)
                        gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                                             >  1:3131461/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:971609/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                                                    >  1:2437332/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                                                    >  1:2691715/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                                                    >  1:2281596/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                                                    >  1:2739706/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                                                    >  1:526767/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                                                    >  1:1715876/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                                                    >  1:1355996/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                                                    >  1:1043125/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                                                    >  1:1628336/1‑52 (MQ=39)
                                   gcaTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGa                                   <  1:2625580/69‑1 (MQ=255)
                                    caTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGa                                   <  1:589439/68‑1 (MQ=255)
                                    caTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGa                                   <  1:3191018/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:447038/1‑69 (MQ=255)
                                        cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag                              >  1:757783/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)
                                        cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag                              >  1:2931011/1‑69 (MQ=255)
                                        cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag                              >  1:1312595/1‑69 (MQ=255)
                                            gCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:2848947/69‑1 (MQ=255)
                                            gCGTCTGCCACAGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          >  1:1431696/1‑69 (MQ=255)
                                               tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga                        <  1:3261932/68‑1 (MQ=255)
                                               tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga                        <  1:2226142/68‑1 (MQ=255)
                                               tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga                        <  1:2181522/68‑1 (MQ=255)
                                               tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga                        <  1:1943021/68‑1 (MQ=255)
                                               tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga                        <  1:634910/68‑1 (MQ=255)
                                               tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga                        <  1:803378/68‑1 (MQ=255)
                                               tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga                        <  1:849140/68‑1 (MQ=255)
                                               tCTGACGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga                        <  1:229727/68‑1 (MQ=255)
                                                      ccGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATggaggag                <  1:125611/69‑1 (MQ=255)
                                                      ccGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATggaggag                <  1:2790894/69‑1 (MQ=255)
                                                      ccGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATggaggag                <  1:546713/69‑1 (MQ=255)
                                                      ccGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATggaggag                <  1:2578659/69‑1 (MQ=255)
                                                             aaGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATggaggag                <  1:1317458/62‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGaa                                  >  1:1673478/1‑43 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          >  1:1892409/1‑51 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:104012/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:1625643/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:3333273/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:3318150/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          >  1:2212987/1‑51 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:1038482/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:3314114/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          >  1:3154337/1‑51 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          >  1:1183226/1‑51 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          >  1:2965185/1‑51 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:2940876/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:2644809/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:224288/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg                          <  1:2807002/51‑1 (MQ=255)
                                                              aGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATACAAAAAGAGtg                          >  1:1742270/1‑51 (MQ=255)
                                                                 gCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATggaggag                <  1:2491115/58‑1 (MQ=255)
                                                                 gCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATggaggag                <  1:1275195/58‑1 (MQ=255)
                                                                  cGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATggaggag                <  1:2235929/57‑1 (MQ=255)
                                                                    tCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATGGAGGAGGACGGCAGCGTAAg  <  1:2872836/69‑1 (MQ=255)
                                                                    tCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATGGAGGAGGACGGCAGCGTAAg  <  1:541108/69‑1 (MQ=255)
                                                                    tCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATGGAGGAGGACGGCAGCGTAAg  <  1:1621776/69‑1 (MQ=255)
                                                                    tCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATGGAGGAGGACGGCAGCGTAAg  <  1:1668146/69‑1 (MQ=255)
                                                                    tCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATGGAGGAGGACGGCAGCGTAAg  <  1:1556007/69‑1 (MQ=255)
                                                                    tCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATGGAGGAGGACGGCAGCGTAAg  <  1:3341340/69‑1 (MQ=255)
                                                                    |                                                                    
GCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGATGGAGGAGGACGGCAGCGTAAG  >  minE/2863399‑2863535

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: