Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,287,932 A→G 65.8% K131K (AAA→AAG argE → acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,287,9320AG65.8% 50.3 / 34.8 67K131K (AAA→AAGargEacetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (9/14);  new base G (21/23);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

AAAGGCTTTTTTGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAGTATGGCCGGAGCGCGT  >  minE/2287879‑2287995
                                                     |                                                               
aaaGGCTTTTTTGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAACTGAAAAAACCGctc                                                  >  1:1408889/1‑69 (MQ=255)
aaaGGCTTTTTTGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAACTGAAAAAACCGctc                                                  >  1:1643826/1‑69 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:2346188/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:2213243/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:2218334/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:3115964/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:2259370/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:3039867/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:229759/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:2332857/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:3146000/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:2397944/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:3032133/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:77610/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:2941056/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:262957/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:692629/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:329679/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:3301215/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:877265/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:3369070/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:897058/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:148317/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:484096/48‑1 (MQ=255)
        ttttGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCt                                                               <  1:1203872/48‑1 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTTGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:359128/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:785130/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:3295825/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:643799/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:3401593/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:2975038/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:691444/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:610336/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:2071676/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:112134/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:130532/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:1007343/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:1534972/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:1621588/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:1919651/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:865099/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:2219164/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:2285237/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:2501863/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:2561551/1‑64 (MQ=255)
                 tATCCTTGATGCGCTACGCGATGTCGACGTCACGAAGCTTGCAAAACCGCTCTACATTCTGGCg                                      >  1:2690324/1‑64 (MQ=255)
                       tGATGCGCTACGCGATGTCGACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCtgatg                            <  1:257457/68‑1 (MQ=255)
                       tGATGCGCTACGCGATGTCGACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCtgatg                            <  1:539621/68‑1 (MQ=255)
                           gcgcTACGCGATGTCGACTTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATgaaga                        <  1:686318/68‑1 (MQ=255)
                           gcgcTACGCGATGTCGACGTCACGAAACTGAAAAAACCGCTTTACATTCTGGCGACTGCTGATgaaga                        <  1:697570/68‑1 (MQ=255)
                           gcgcTACGCGATGTCGACGTCACGAAACTGAAAAAACCGCTGTACATTCTGGCGACTGCTGATgaaga                        <  1:1737918/68‑1 (MQ=255)
                           gcgcTACGCGATGTCGACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATgaaga                        <  1:662729/68‑1 (MQ=255)
                           gcgcTACGCGATGTCGACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATgaaga                        <  1:185267/68‑1 (MQ=255)
                           gcgcTACGCGATGTCGACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATgaaga                        <  1:473352/68‑1 (MQ=255)
                                aCGCGATGTCGACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAg                  <  1:975876/69‑1 (MQ=255)
                                aCGCGATGTCGACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAg                  <  1:3216340/69‑1 (MQ=255)
                                          gACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAGTATGGCCGGa        >  1:1339470/1‑69 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAg                  <  1:3128376/53‑1 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAg                  <  1:2838520/53‑1 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAg                  <  1:2182228/53‑1 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAg                  <  1:1728843/53‑1 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAGTATGGCCGGAGCGCGt  >  1:337984/1‑69 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAGTATGGCCGGAGCGCGt  >  1:2785235/1‑69 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAGTATGGCCGGAGCGCGt  >  1:2765703/1‑69 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAGTATGGCCGGAGCGCGt  >  1:2518734/1‑69 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAGTATGGCCGGAGCGCGt  >  1:2514639/1‑69 (MQ=255)
                                                aCGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAGTATGGCCGGAGCGCGt  >  1:1762134/1‑69 (MQ=255)
                                                     |                                                               
AAAGGCTTTTTTGCGTTTATCCTTGATGCGCTACGCGATGTCGACGTCACGAAACTGAAAAAACCGCTCTACATTCTGGCGACTGCTGATGAAGAAACCAGTATGGCCGGAGCGCGT  >  minE/2287879‑2287995

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