Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,526,757 G→A 27.4% N75N (AAC→AAT ttk ← division inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,526,7570GA27.4% 124.7 / 43.7 109N75N (AAC→AATttkdivision inhibitor
Reads supporting (aligned to +/- strand):  ref base G (21/58);  new base A (14/16);  total (35/74)
Fisher's exact test for biased strand distribution p-value = 6.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

Only 100 of 109 total aligned reads displayed.

ACCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATGCGGGTCTT  >  minE/2526691‑2526819
                                                                  |                                                              
aCCGAGAAGCAGCAACACAATCAGACGCAGGCTCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGtt                                                              <  1:821380/69‑1 (MQ=255)
aCCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGtt                                                              <  1:1483525/69‑1 (MQ=255)
aCCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGtt                                                              <  1:109499/69‑1 (MQ=255)
aCCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGtt                                                              <  1:7420/69‑1 (MQ=255)
aCCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGtt                                                              <  1:946605/69‑1 (MQ=255)
aCCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGtt                                                              <  1:155282/69‑1 (MQ=255)
aCCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGtt                                                              <  1:355946/69‑1 (MQ=255)
aCCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGtt                                                              <  1:2995452/69‑1 (MQ=255)
aCCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGtt                                                              <  1:2399082/69‑1 (MQ=255)
 ccGAGAATCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTg                                                             <  1:197582/69‑1 (MQ=255)
 ccGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTg                                                             <  1:877751/69‑1 (MQ=255)
         cagcaACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGt                                                     <  1:1859593/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:1709694/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:2691941/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:2684577/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:2184032/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:364284/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:1493476/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:141632/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:85310/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:299291/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:2190975/69‑1 (MQ=255)
                  aaTCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:3346965/69‑1 (MQ=255)
                                        gtgtCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCa                        >  1:2849008/1‑67 (MQ=255)
                                        gtgtCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAg                       >  1:1080901/1‑68 (MQ=255)
                                        gtgtCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAg                       >  1:2192496/1‑68 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGc                     <  1:2833793/67‑1 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGc                     <  1:2562792/67‑1 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGc                     <  1:676214/67‑1 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGc                     <  1:2627120/67‑1 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGc                     <  1:2799574/67‑1 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGc                     <  1:1145260/67‑1 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGc                     <  1:3043231/67‑1 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGc                     <  1:337888/67‑1 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTa                   >  1:1928442/1‑69 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTa                   >  1:2219135/1‑69 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTa                   >  1:1175859/1‑69 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTa                   >  1:2015233/1‑69 (MQ=255)
                                           tCTTTCTCATCTTTCAGAATCAGGTTGATGCAAGTAATCAGGCTATCTTCGATAAACTCAATCAGGc                     <  1:2437383/67‑1 (MQ=255)
                                             tttCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCt                                            <  1:94947/42‑1 (MQ=255)
                                             tttCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTAt                  >  1:3154099/1‑68 (MQ=255)
                                             tttCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTAt                  >  1:2553054/1‑68 (MQ=255)
                                             tttCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTAt                  >  1:1848293/1‑68 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:1804767/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:3151369/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:114685/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:11983/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:1539970/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:1661995/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:557800/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:3295062/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:1844357/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:1128286/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:1915101/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:2250324/44‑1 (MQ=255)
                                                ctcATCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCtt                                       <  1:2035486/44‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCg                <  1:509375/66‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCg                <  1:2153097/66‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCg                <  1:2709916/66‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCg                <  1:620939/66‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCg                <  1:174758/66‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCg                <  1:3070499/66‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:2121063/68‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:1878102/68‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:1138309/68‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:2175938/68‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:3225219/68‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:1629439/68‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:3354323/68‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:2824174/68‑1 (MQ=255)
                                                 tcATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:3299431/68‑1 (MQ=255)
                                                  cATCTTTCAGAATCATATTGATGCGAGTAATCAGGCAATCtt                                       <  1:2863257/42‑1 (MQ=37)
                                                  cATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGaa              <  1:571941/67‑1 (MQ=255)
                                                   aTCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCaa                           >  1:374521/1‑53 (MQ=255)
                                                   aTCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCaa                           >  1:1366366/1‑53 (MQ=255)
                                                   aTCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCaa                           >  1:72533/1‑53 (MQ=255)
                                                   aTCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCaa                           >  1:2517820/1‑53 (MQ=255)
                                                   aTCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCaa                           >  1:1996350/1‑53 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:3435567/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:750602/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:233904/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:2979964/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:1207502/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:1347813/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:1780174/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:1826170/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:1839546/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:1971033/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:2037147/1‑69 (MQ=255)
                                                    tCTTTCAGAATCAGATTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATg          >  1:2068606/1‑69 (MQ=255)
                                                      tttCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATa                                  <  1:287209/43‑1 (MQ=255)
                                                      tttCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATa                                  <  1:3112686/43‑1 (MQ=255)
                                                      tttCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATa                                  <  1:1134856/43‑1 (MQ=255)
                                                      tttCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATa                                  <  1:1087587/43‑1 (MQ=255)
                                                      tttCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATGCg        >  1:1838185/1‑69 (MQ=255)
                                                      tttCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATGCg        >  1:271073/1‑69 (MQ=255)
                                                      tttCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATGCg        >  1:360080/1‑69 (MQ=255)
                                                      tttCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATGCg        >  1:876636/1‑69 (MQ=255)
                                                            aaTCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATGCGGGTCtt  <  1:2258197/69‑1 (MQ=255)
                                                            aaTCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATGCGGGTCtt  <  1:1194539/69‑1 (MQ=255)
                                                                  |                                                              
ACCGAGAAGCAGCAACACAATCAGACGCAGGCGCGCTGTGGTGTCTTTCTCATCTTTCAGAATCAGGTTGATGCGAGTAATCAGGCTATCTTCGATAAACTCAATCAGGCTATCGAACATGCGGGTCTT  >  minE/2526691‑2526819

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