Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,460,079 A→G 37.6% S47S (TCA→TCG mioC → FMN‑binding protein MioC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,460,0790AG37.6% 50.0 / 60.2 93S47S (TCA→TCGmioCFMN‑binding protein MioC
Reads supporting (aligned to +/- strand):  ref base A (33/25);  new base G (16/19);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CTGGAAGAGGCGGGTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAACCTTTCTCCTTTC  >  minE/2460014‑2460145
                                                                 |                                                                  
cTGGAAGAGGCGGGTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 <  1:1072327/69‑1 (MQ=255)
cTGGAAGAGGCGGGTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 <  1:1235290/69‑1 (MQ=255)
cTGGAAGAGGCGGGTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 <  1:488111/69‑1 (MQ=255)
cTGGAAGAGGCGGGTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 <  1:794282/69‑1 (MQ=255)
         gCGGGTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAGGGATctggctg                                                        >  1:3271272/1‑69 (MQ=255)
           gggTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 >  1:1012700/1‑58 (MQ=255)
           gggTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 >  1:229348/1‑58 (MQ=255)
           gggTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 >  1:2194494/1‑58 (MQ=255)
           gggTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 >  1:216380/1‑58 (MQ=255)
           gggTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 >  1:1394532/1‑58 (MQ=255)
           gggTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 >  1:1945088/1‑58 (MQ=255)
           gggTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 >  1:1903086/1‑58 (MQ=255)
           gggTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAggg                                                                 >  1:1520518/1‑58 (MQ=255)
               ttaaccaccGAAACGCTGCACGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATc                                                  <  1:103579/66‑1 (MQ=255)
                     accGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTc                                             <  1:785043/68‑1 (MQ=255)
                     accGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTc                                             <  1:2671143/68‑1 (MQ=255)
                     accGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTc                                             <  1:2482238/68‑1 (MQ=255)
                     accGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTc                                             <  1:1806349/68‑1 (MQ=255)
                     accGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTc                                             <  1:1421007/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:3423160/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:3023412/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:2727910/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:252808/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:329879/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:2188496/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:2184651/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:2034629/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:1969910/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:370396/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:508620/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:564386/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:1647397/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:1582226/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:597756/68‑1 (MQ=255)
                                   cGGTCCGCTGTTGGAAGATCTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:335690/68‑1 (MQ=255)
                                   cGGTCCGCAGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:1815621/68‑1 (MQ=255)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGTGGATCTGGCTGGTTATCa                                                 >  1:1201062/1‑49 (MQ=255)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGtt                                                     >  1:2090618/1‑45 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:2713680/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:3114574/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:640735/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:1701365/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:3402655/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:712068/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:3422513/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:2271351/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:2771431/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:835257/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:712686/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:2151655/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTATCa                                                 >  1:368019/1‑49 (MQ=38)
                                    ggTCCGCTGTTGGAAGATTTATCTGCCTCGGGGATCTGGCTGGTTACCa                                                 >  1:1457571/1‑49 (MQ=37)
                                       ccGCTGTTAGAAGATTTACCTGCCACAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGgaga                           <  1:2580758/68‑1 (MQ=255)
                                             ttAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTc                                             <  1:88701/44‑1 (MQ=255)
                                              tAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCg                    >  1:1382401/1‑68 (MQ=255)
                                              tAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCg                    >  1:581947/1‑68 (MQ=255)
                                              tAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCg                    >  1:2930548/1‑68 (MQ=255)
                                               agaagaTTTACCTTCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCgg                   >  1:97500/1‑68 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCgg                   >  1:2481422/1‑68 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCgg                   >  1:2161013/1‑68 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCgg                   >  1:1807040/1‑68 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCgg                   >  1:3398232/1‑68 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCgg                   >  1:3307243/1‑68 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCgg                   >  1:2699147/1‑68 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCgg                   >  1:1470569/1‑68 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCgg                   >  1:280796/1‑68 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCg                    >  1:3129765/1‑67 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCg                    >  1:2463078/1‑67 (MQ=255)
                                               agaagaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACACACGGTGCCGGAGATATTCCgg                   >  1:2092786/1‑68 (MQ=255)
                                                   gaTTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAc              <  1:2284836/69‑1 (MQ=255)
                                                      ttACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctt           >  1:671422/1‑69 (MQ=255)
                                                      ttACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctt           >  1:2529484/1‑69 (MQ=255)
                                                      ttACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctt           >  1:1670049/1‑69 (MQ=255)
                                                      ttACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctt           >  1:846657/1‑69 (MQ=255)
                                                        aCCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGat                          >  1:1065190/1‑52 (MQ=255)
                                                        aCCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttc         >  1:2922818/1‑69 (MQ=255)
                                                        aCCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttc         >  1:1454220/1‑69 (MQ=255)
                                                        aCCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttc         >  1:3310688/1‑69 (MQ=255)
                                                            gCCTCGGGGATCTGGCTGGTTATCAGCTCAACCCACGGTGCCg                               <  1:1954940/43‑1 (MQ=38)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       <  1:3143178/67‑1 (MQ=255)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       <  1:1192587/67‑1 (MQ=255)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       <  1:2769248/67‑1 (MQ=255)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       <  1:1503907/67‑1 (MQ=255)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       <  1:1556103/67‑1 (MQ=255)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       <  1:362995/67‑1 (MQ=255)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       <  1:2388681/67‑1 (MQ=255)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       >  1:1878505/1‑67 (MQ=255)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       <  1:1945410/67‑1 (MQ=255)
                                                            gCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctc       <  1:3410578/67‑1 (MQ=255)
                                                               tCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctcctttc  <  1:3037246/69‑1 (MQ=255)
                                                               tCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctcctttc  <  1:190118/69‑1 (MQ=255)
                                                               tCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAAcctttctcctttc  <  1:1956288/69‑1 (MQ=255)
                                                               tCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGAGCCGGAGATATTCCGGACAAcctttctcctttc  <  1:1304350/69‑1 (MQ=255)
                                                                 |                                                                  
CTGGAAGAGGCGGGTTTTACCACCGAAACGCTGCACGGTCCGCTGTTAGAAGATTTACCTGCCTCAGGGATCTGGCTGGTTATCAGCTCCACCCACGGTGCCGGAGATATTCCGGACAACCTTTCTCCTTTC  >  minE/2460014‑2460145

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: