Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 820,100 C→T 54.0% A189A (GCC→GCT ychA → predicted transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE820,1000CT54.0% 8.6 / 73.8 87A189A (GCC→GCTychApredicted transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (4/36);  new base T (29/18);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 5.55e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAC  >  minE/820038‑820166
                                                              |                                                                  
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAATCCTCGtt                                                               <  1:2242342/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:3339150/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:2371857/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:2421606/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:1839197/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:1799965/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:243856/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:357265/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:2280810/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:3334932/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:2741982/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:286520/68‑1 (MQ=255)
gAAGACCTTGATGAAGCTGATAACATTGAGGTAATCAGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:2489501/68‑1 (MQ=255)
  aGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:2205658/66‑1 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGGTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:3142492/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGAtag                                                          >  1:2296400/1‑51 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGTAAATCAGATGGa                                          >  1:2996701/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATgg                                           >  1:1750673/1‑68 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATgg                                           >  1:452125/1‑68 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:2187946/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:230726/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:3255566/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:2364501/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:2394744/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:3379040/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:2676576/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:2770988/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:2814691/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:302621/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:323228/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:102710/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:12256/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:1682604/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:795963/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:862200/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:577227/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:1841854/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:1944431/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:195043/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:199221/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:2032344/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGa                                          >  1:206552/1‑69 (MQ=255)
                    tAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCTTAGTTGATAGAAGAAAATCAGATGGa                                          >  1:482524/1‑69 (MQ=255)
                       cATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:117105/45‑1 (MQ=255)
                             ggTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGtt                                 >  1:3285550/1‑69 (MQ=255)
                             ggTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGtt                                 >  1:1173825/1‑69 (MQ=255)
                             ggTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGtt                                 >  1:2589083/1‑69 (MQ=255)
                             ggTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGtt                                 >  1:1565761/1‑69 (MQ=255)
                                 aTCCGCAAATTGCTGGATACACTCAAAGCCTCGtt                                                               <  1:1969502/35‑1 (MQ=255)
                                         aTTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGa                     <  1:3127987/69‑1 (MQ=255)
                                         aTTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGa                     <  1:2838062/69‑1 (MQ=255)
                                         aTTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGa                     <  1:476679/69‑1 (MQ=255)
                                         aTTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGa                     <  1:836538/69‑1 (MQ=255)
                                         aTTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGa                     <  1:2168116/69‑1 (MQ=255)
                                               ggATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCtt                <  1:3375240/68‑1 (MQ=255)
                                               ggATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCtt                <  1:577813/68‑1 (MQ=255)
                                               ggATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCtt                <  1:1780894/68‑1 (MQ=255)
                                               ggATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCtt                <  1:3074229/68‑1 (MQ=255)
                                               ggATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCtt                <  1:1911357/68‑1 (MQ=255)
                                               ggATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCtt                <  1:1146641/68‑1 (MQ=255)
                                                  tACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTttat             <  1:1635508/68‑1 (MQ=255)
                                                  tACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTttat             <  1:79652/68‑1 (MQ=255)
                                                  tACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTttat             <  1:2889277/68‑1 (MQ=255)
                                                  tACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTttat             <  1:366349/68‑1 (MQ=255)
                                                  tACACTCAAAGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTttat             <  1:2324733/68‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:3331161/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:3290489/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:1119139/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:2836885/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:2697050/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:1538599/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:2634684/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:2621205/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:2617990/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:1671353/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:2201644/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:1736004/67‑1 (MQ=255)
                                                           aGCTTCGTTGATAGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTc     <  1:1756516/67‑1 (MQ=255)
                                                            gCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAc  <  1:1187624/69‑1 (MQ=255)
                                                            gCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAc  <  1:3117468/69‑1 (MQ=255)
                                                            gCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAc  <  1:1750337/69‑1 (MQ=255)
                                                            gCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAc  <  1:542390/69‑1 (MQ=255)
                                                            gCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAc  <  1:184906/69‑1 (MQ=255)
                                                            gCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAc  <  1:2096251/69‑1 (MQ=255)
                                                            gCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAc  <  1:247604/69‑1 (MQ=255)
                                                            gCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAc  <  1:2668380/69‑1 (MQ=255)
                                                            gCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAc  <  1:2976890/69‑1 (MQ=255)
                                                              |                                                                  
GAAGACCTTGATGAAGCTGATAACATTGAGGTAATCCGCAAATTGCTGGATACACTCAAAGCCTCGTTGATGGAAGAAAATCAGATGGAGCTGGCGTTACGCACCAGCGAAGCTTTATTACAATTCAAC  >  minE/820038‑820166

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: