Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,461,663 A→G 50.0% G301G (GGA→GGG gidA → glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,461,6630AG50.0% ‑1.0 / 96.7 106G301G (GGA→GGGgidAglucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base A (49/4);  new base G (10/43);  total (59/47)
Fisher's exact test for biased strand distribution p-value = 3.64e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

GATCGAAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAG  >  minE/2461597‑2461723
                                                                  |                                                            
gatcgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACt                                                            >  1:2947022/1‑69 (MQ=255)
gatcgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACt                                                            >  1:2920423/1‑69 (MQ=255)
gatcgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACt                                                            >  1:677003/1‑69 (MQ=255)
gatcgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACt                                                            >  1:2146714/1‑69 (MQ=255)
gatcgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACt                                                            >  1:1244060/1‑69 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGa                                                              <  1:508416/64‑1 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGa                                                              <  1:2860108/64‑1 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGa                                                              <  1:197373/64‑1 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGa                                                              <  1:3172848/64‑1 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGa                                                              <  1:1899397/64‑1 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGa                                                              <  1:1471225/64‑1 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGa                                                              <  1:951855/64‑1 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGa                                                              <  1:1035625/64‑1 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGAc                                                         >  1:2278736/1‑69 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGAc                                                         >  1:1305713/1‑69 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGAc                                                         >  1:1845993/1‑69 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGAc                                                         >  1:3044192/1‑69 (MQ=255)
   cgaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGAc                                                         >  1:2317583/1‑69 (MQ=255)
    gaAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGAc                                                         <  1:3226642/68‑1 (MQ=255)
                   cgcTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTAt                                         >  1:2796926/1‑69 (MQ=255)
                     cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                        >  1:2449035/1‑68 (MQ=255)
                     cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                        >  1:777639/1‑68 (MQ=255)
                     cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                        >  1:2704592/1‑68 (MQ=255)
                     cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                        >  1:1128729/1‑68 (MQ=255)
                     cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                        >  1:838567/1‑68 (MQ=255)
                     cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                        >  1:870113/1‑68 (MQ=255)
                     cTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATc                                        >  1:849215/1‑68 (MQ=255)
                       tCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGa                                     >  1:3409260/1‑69 (MQ=255)
                       tCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGa                                     >  1:1961185/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                    >  1:661500/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                    >  1:1059796/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                    >  1:3119434/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                    >  1:1003799/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                    >  1:857918/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                    >  1:1062831/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                    >  1:1935978/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                    >  1:513842/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGaa                                    >  1:2885619/1‑69 (MQ=255)
                        cGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACATCTAACGAAATTTATCCGaa                                    >  1:877584/1‑69 (MQ=255)
                           cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:2847136/68‑1 (MQ=255)
                           cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:1298133/68‑1 (MQ=255)
                           cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:311100/68‑1 (MQ=255)
                           cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:1488040/68‑1 (MQ=255)
                           cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:771337/68‑1 (MQ=255)
                           cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:1064483/68‑1 (MQ=255)
                           cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:554866/68‑1 (MQ=255)
                           cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:1565415/68‑1 (MQ=255)
                           cGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:3138459/68‑1 (MQ=255)
                            gACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGAAATTTATCCGAACg                                  <  1:407897/67‑1 (MQ=255)
                            gACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGt                                <  1:3093624/69‑1 (MQ=255)
                             aCAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTa                               <  1:43411/69‑1 (MQ=255)
                                            aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                           >  1:3419445/1‑58 (MQ=255)
                                            aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                           >  1:2790551/1‑58 (MQ=255)
                                            aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                           >  1:261406/1‑58 (MQ=255)
                                            aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                           >  1:2369010/1‑58 (MQ=255)
                                            aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                           >  1:3189238/1‑58 (MQ=255)
                                            aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                           >  1:3209019/1‑58 (MQ=255)
                                            aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtctc                           >  1:3328173/1‑58 (MQ=255)
                                            aGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtc                             >  1:3293823/1‑56 (MQ=255)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:421047/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:3037914/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:960816/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:2887380/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:110876/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:949210/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:128514/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:1340552/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:1374013/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:1420104/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:1420283/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:1706440/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:1884817/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:1902937/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:2219123/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:2308313/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:2508339/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:2866280/38‑1 (MQ=38)
                                             gATCTTCCTTGAACCGGAAGGGCTGACCTCTAACGaaa                                              <  1:2887348/38‑1 (MQ=38)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGATCGGTATCTCCACCAGCCTGCCGt             >  1:1643278/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTAtc                             >  1:2943073/1‑52 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCCGCCTGCCGt             >  1:1218204/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCt                  >  1:1489215/1‑63 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:2753029/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:482732/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:470518/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:466382/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:2831568/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:2617141/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:2294326/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:2152802/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:1911685/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:146878/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:113703/1‑68 (MQ=255)
                                                cTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGt             >  1:1004410/1‑68 (MQ=255)
                                                          ccGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg  <  1:3383815/69‑1 (MQ=255)
                                                          ccGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg  <  1:2420923/69‑1 (MQ=255)
                                                          ccGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg  <  1:1045297/69‑1 (MQ=255)
                                                          ccGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCa   <  1:2319861/68‑1 (MQ=255)
                                                           cGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg  <  1:3274176/68‑1 (MQ=255)
                                                           cGGAAGGGCTGACCTCTAACGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAg  <  1:1911886/68‑1 (MQ=255)
                                                                  |                                                            
GATCGAAGACAAAGTCATGCGCTTCGCCGACAGAAATCAGCATCAGATCTTCCTTGAACCGGAAGGACTGACCTCTAATGAAATTTATCCGAACGGTATCTCCACCAGCCTGCCGTTCGATGTGCAG  >  minE/2461597‑2461723

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