Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,184,950 G→A 40.5% E305E (GAG→GAA yeaG → conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,184,9500GA40.5% 46.0 / 70.5 102E305E (GAG→GAAyeaGconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base A (20/21);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 102 total aligned reads displayed.

TACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACA  >  minE/1184885‑1185011
                                                                 |                                                             
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:1810067/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:2440961/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:251209/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:10002/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:2717705/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:2179299/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:2743781/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:2844698/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:1977447/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:28702/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:2900732/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:3320269/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:1689681/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:3329638/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:3358513/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:8789/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:1542345/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:705320/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATc                                                             >  1:707172/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTc                                                             >  1:2338780/1‑68 (MQ=255)
tACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTc                                                             >  1:1045779/1‑68 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:1232/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:101381/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:2345941/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:263091/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:2675613/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:2252624/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:1120019/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:2147796/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:2789876/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:3384422/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:304005/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:3065703/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:1898560/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:2868074/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:1795405/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:3267920/68‑1 (MQ=255)
    aCGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:1510141/68‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:2607589/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:3357720/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:3353731/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:847954/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:491937/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:78149/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:2809033/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:77319/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:2906404/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1565782/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1648191/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:2478477/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1651484/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1043014/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1197374/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1412742/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1324747/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1286716/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1936202/69‑1 (MQ=255)
              aGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:1955880/69‑1 (MQ=255)
                 tATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTcc                                            >  1:2370141/1‑68 (MQ=255)
                 tATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTcc                                            >  1:3329072/1‑68 (MQ=255)
                 tATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTcc                                            >  1:764484/1‑68 (MQ=255)
                 tATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTcc                                            >  1:1948155/1‑68 (MQ=255)
                 tATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTcc                                            >  1:2039434/1‑68 (MQ=255)
                 tATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTcc                                            >  1:2216135/1‑68 (MQ=255)
                 tATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTc                                             >  1:1928173/1‑67 (MQ=255)
                 tATCTCAGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTcc                                            >  1:435550/1‑68 (MQ=255)
                     tcCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGtaat                                       >  1:196967/1‑69 (MQ=255)
                     tcCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGtaat                                       >  1:215545/1‑69 (MQ=255)
                     tcCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGtaat                                       >  1:1738319/1‑69 (MQ=255)
                           cTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:2635563/45‑1 (MQ=255)
                                     aCGGGATTATTCTCGCACACTCGAACGAGTCCGaa                                                         <  1:2701793/35‑1 (MQ=255)
                                            tattCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTcc                 >  1:1987046/1‑68 (MQ=255)
                                            tattCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTcc                 >  1:3376656/1‑68 (MQ=255)
                                            tattCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTcc                 >  1:3371534/1‑68 (MQ=255)
                                            tattCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTcc                 >  1:2305504/1‑68 (MQ=255)
                                            tattCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTcc                 >  1:638514/1‑68 (MQ=255)
                                            tattCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTcc                 >  1:3071976/1‑68 (MQ=255)
                                            tattCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTcc                 >  1:1211096/1‑68 (MQ=255)
                                            tattCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTc                  >  1:253470/1‑67 (MQ=255)
                                            tattCTCGCACACTCGAACGAGTCCGAATCGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTcc                 >  1:460433/1‑68 (MQ=255)
                                               tctcGCCCACTCGAACGAATCCGAATGGGTCACttt                                              <  1:38401/36‑1 (MQ=255)
                                               tctcGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCg              <  1:607930/68‑1 (MQ=255)
                                               tctcGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCg              <  1:2616770/68‑1 (MQ=255)
                                               tctcGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCg              <  1:3358334/68‑1 (MQ=255)
                                               tctcGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCg              <  1:208610/68‑1 (MQ=255)
                                               tctcGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCg              <  1:151118/68‑1 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGttt     >  1:1012753/1‑68 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTa    >  1:3357219/1‑69 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTa    >  1:17789/1‑69 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTa    >  1:3222750/1‑69 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTa    >  1:3078213/1‑69 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTa    >  1:19315/1‑69 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTa    >  1:2968543/1‑69 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTa    >  1:2651605/1‑69 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTa    >  1:2318745/1‑69 (MQ=255)
                                                        cTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTa    >  1:23408/1‑69 (MQ=255)
                                                          cGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACa  <  1:1621696/69‑1 (MQ=255)
                                                          cGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACa  <  1:3128893/69‑1 (MQ=255)
                                                          cGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACa  <  1:3093157/69‑1 (MQ=255)
                                                          cGAACGACTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACa  <  1:677684/69‑1 (MQ=255)
                                                                 |                                                             
TACAACGGGACGGAAGGTATCTCCGCCCTGCCGTTCAACGGGATTATTCTCGCACACTCGAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACA  >  minE/1184885‑1185011

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: