Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,052,195 C→T 46.9% R92R (CGC→CGT htrG → predicted signal transduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,052,1950CT46.9% 9.3 / 69.2 83R92R (CGC→CGThtrGpredicted signal transduction protein
Reads supporting (aligned to +/- strand):  ref base C (29/15);  new base T (20/19);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

GAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACA  >  minE/2052129‑2052262
                                                                  |                                                                   
gAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTc                                                                   >  1:1854520/1‑69 (MQ=255)
gAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTc                                                                   >  1:552795/1‑69 (MQ=255)
gAAAGACAGCTCTGGCCGTACCGCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTcc                                                                  >  1:994272/1‑69 (MQ=255)
    gACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCgt                                                               >  1:176999/1‑69 (MQ=255)
    gACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCgt                                                               >  1:2093513/1‑69 (MQ=255)
    gACAGCTCTGGCCGTACCGCCTGGATCCCGGTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCgt                                                               >  1:474715/1‑69 (MQ=255)
     aCAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCgtg                                                              >  1:2047982/1‑69 (MQ=255)
     aCAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCgt                                                               >  1:921722/1‑68 (MQ=255)
          tctGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCgtgt                                                             >  1:401198/1‑65 (MQ=255)
          tctGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCgtgt                                                             >  1:1751437/1‑65 (MQ=255)
                 gTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTg                                                   >  1:837957/1‑68 (MQ=255)
                   aCCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaa                                                <  1:422412/69‑1 (MQ=255)
                   aCCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaa                                                <  1:1358895/69‑1 (MQ=255)
                   aCCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGaa                                                <  1:2466156/69‑1 (MQ=255)
                              cccGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCaa                                     >  1:387603/1‑69 (MQ=255)
                              cccGTTGAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCaa                                     <  1:610702/69‑1 (MQ=255)
                              cccGTTGAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCaa                                     <  1:2908005/69‑1 (MQ=255)
                              cccGTTGAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCaa                                     <  1:2206296/69‑1 (MQ=255)
                              cccGTTGAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCaa                                     <  1:177103/69‑1 (MQ=255)
                              cccGTTGAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCaa                                     <  1:2856554/69‑1 (MQ=255)
                              cccGTTGAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCaa                                     <  1:1290993/69‑1 (MQ=255)
                                 gTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCaaaa                                   <  1:1156252/68‑1 (MQ=255)
                                 gTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCaaaa                                   <  1:578531/68‑1 (MQ=255)
                                 gTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCaaa                                    >  1:1545274/1‑67 (MQ=255)
                                 gTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAc                                  >  1:1838190/1‑69 (MQ=255)
                                 gTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAc                                  >  1:889198/1‑69 (MQ=255)
                                 gTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAc                                  >  1:2760166/1‑69 (MQ=255)
                                 gTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAc                                  >  1:3346484/1‑69 (MQ=255)
                                 gTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAc                                  >  1:3369797/1‑69 (MQ=255)
                                 gTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAAc                                  >  1:51066/1‑69 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:3062858/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:3271823/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:3314006/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:358589/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:2410528/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:546675/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:1166628/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:1608085/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:649651/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:1963080/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:992142/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGGAACCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:3139735/45‑1 (MQ=255)
                                    gAAACAACTTAGCACGCAGCCAAGCCTGCGTTCCCGTGTGCCAGa                                                       <  1:1727733/45‑1 (MQ=255)
                                                 aCTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATaaa                     <  1:445900/66‑1 (MQ=255)
                                                 aCTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   <  1:472195/68‑1 (MQ=255)
                                                 aCTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   <  1:3286648/68‑1 (MQ=255)
                                                 aCTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   <  1:1993609/68‑1 (MQ=255)
                                                 aCTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   <  1:209767/68‑1 (MQ=255)
                                                 aCTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   <  1:2118096/68‑1 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAAt                                             >  1:1916895/2‑40 (MQ=38)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTTACCGATAAACTCa                 >  1:144118/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTc                  >  1:2690511/2‑67 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:447361/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:1074637/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:1091214/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:885077/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:733284/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:2072938/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:584302/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:1228668/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:472019/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:125483/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:1398410/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:2547223/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:3434548/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:3427648/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:3373493/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:2737829/2‑68 (MQ=255)
                                                   ggAGCCAAGCCTGCGTTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACGCTGACCGATAAACTCa                 >  1:3301968/2‑68 (MQ=255)
                                                    gAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCAcc               <  1:3003128/69‑1 (MQ=255)
                                                          aGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATc         >  1:839432/1‑69 (MQ=255)
                                                          aGCCTGCGCTCCCGTGTGCC‑‑ATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCg        >  1:1218929/1‑68 (MQ=255)
                                                            ccTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   >  1:3044187/1‑57 (MQ=255)
                                                            ccTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   >  1:1682222/1‑57 (MQ=255)
                                                            ccTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   >  1:3098864/1‑57 (MQ=255)
                                                            ccTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   >  1:3121298/1‑57 (MQ=255)
                                                            ccTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   >  1:85713/1‑57 (MQ=255)
                                                            ccTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   >  1:2847302/1‑57 (MQ=255)
                                                            ccTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   >  1:2017275/1‑57 (MQ=255)
                                                            ccTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACt                   >  1:1470837/1‑57 (MQ=255)
                                                             cTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGAt      <  1:2742606/69‑1 (MQ=255)
                                                             cTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGAt      <  1:1497412/69‑1 (MQ=255)
                                                             cTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGAt      <  1:457634/69‑1 (MQ=255)
                                                                  cTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAaca  <  1:1014060/68‑1 (MQ=255)
                                                                  cTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAaca  <  1:293917/68‑1 (MQ=255)
                                                                  cTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAaca  <  1:2498641/68‑1 (MQ=255)
                                                                  cTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAaca  <  1:2544401/68‑1 (MQ=255)
                                                                  |                                                                   
GAAAGACAGCTCTGGCCGTACCGCCTGGATCCCGTTGAAACAACTTAGCACTGAGCCAAGCCTGCGCTCCCGTGTGCCAGATCTGGAAAATCAGGTCAAAACCCTGACCGATAAACTCACCAATATCGATAACA  >  minE/2052129‑2052262

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: