Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,916,684 T→C 34.8% intergenic (‑670/‑15) nudH ← / → mutH nucleotide hydrolase/methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,916,6840TC34.8% 57.3 / 49.7 86intergenic (‑670/‑15)nudH/mutHnucleotide hydrolase/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (48/8);  new base C (10/20);  total (58/28)
Fisher's exact test for biased strand distribution p-value = 1.59e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACA  >  minE/1916619‑1916745
                                                                 |                                                             
gAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCaa                                                            >  1:3141162/1‑69 (MQ=255)
gAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCaa                                                            >  1:3177256/1‑69 (MQ=255)
gAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCaa                                                            >  1:2151195/1‑69 (MQ=255)
     acGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTAt                                                       >  1:2908822/1‑69 (MQ=255)
     acGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTAt                                                       >  1:2346379/1‑69 (MQ=255)
     acGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTAt                                                       >  1:803924/1‑69 (MQ=255)
     acGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTAt                                                       >  1:3301992/1‑69 (MQ=255)
     acGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTAt                                                       >  1:1341945/1‑69 (MQ=255)
          tgtgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGa                                                  <  1:474802/69‑1 (MQ=255)
            tgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                >  1:92395/1‑69 (MQ=255)
            tgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                >  1:101285/1‑69 (MQ=255)
            tgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                >  1:2445237/1‑69 (MQ=255)
            tgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                >  1:2433666/1‑69 (MQ=255)
            tgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                >  1:2101336/1‑69 (MQ=255)
            tgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                >  1:3418467/1‑69 (MQ=255)
            tgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                >  1:1330235/1‑69 (MQ=255)
                  aaGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTccc                                          <  1:2603994/69‑1 (MQ=255)
                                   ttCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctc                         >  1:37389/1‑69 (MQ=255)
                                   ttCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctc                         >  1:631529/1‑69 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGt                                             >  1:1592006/1‑40 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:1418056/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:1127522/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:860819/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:2985867/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:2214855/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:2163831/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:1292838/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:1297283/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:1690955/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:460846/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                 >  1:1503157/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtccacc                 >  1:2723886/1‑68 (MQ=255)
                                            tAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCCCCCACTGCTCTCtcctcc                 >  1:881156/1‑68 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTccc                >  1:667774/1‑67 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGa              >  1:1210185/1‑69 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGa              >  1:2768986/1‑69 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGa              >  1:2548065/1‑69 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:2730219/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:2880963/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:3010641/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:3128750/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:3138081/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:3152564/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:605753/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:437791/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:43882/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:2624804/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:1241095/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:1304840/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:131458/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:1334172/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:134028/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:2104235/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:216172/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:2613348/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:2530109/52‑1 (MQ=255)
                                              aTTGCTGTTTGTTTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCa                               <  1:2815240/52‑1 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:2380709/1‑67 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:910714/1‑67 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:22995/1‑67 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:3171726/1‑67 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:1192827/1‑67 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:301289/1‑67 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:855861/1‑67 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:1584090/1‑67 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:1569370/1‑67 (MQ=255)
                                                   tgtttgttTTTTAACCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAAc           >  1:3007337/1‑67 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCACCCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:1232807/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGTCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:300075/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAATCTgaag      >  1:592373/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:809776/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:1044380/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:61939/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:532003/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:605133/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:512358/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:31049/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:230377/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaag      >  1:2356256/1‑68 (MQ=255)
                                                       tgtttTTTAATCAAGGTATCATGACATGTCCCAACCTCGACCACTGCTCTCTCCTCCCGAAACTgaag      >  1:1101590/1‑68 (MQ=255)
                                                           ttttAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACa  <  1:548213/68‑1 (MQ=255)
                                                           ttttAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACa  <  1:494823/68‑1 (MQ=255)
                                                           ttttAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACa  <  1:723439/68‑1 (MQ=255)
                                                           ttttAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACa  <  1:15684/68‑1 (MQ=255)
                                                           ttttAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACa  <  1:1801303/68‑1 (MQ=255)
                                                           ttttAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACa  <  1:3360916/68‑1 (MQ=255)
                                                                 |                                                             
GAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACA  >  minE/1916619‑1916745

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: