Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 389,226 T→C 45.3% C248C (TGT→TGC ybdL → methionine aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE389,2260TC45.3% 12.0 / 71.0 86C248C (TGT→TGCybdLmethionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base C (9/30);  total (34/52)
Fisher's exact test for biased strand distribution p-value = 7.44e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

AGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAAT  >  minE/389159‑389292
                                                                   |                                                                  
aGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATtgtg                                                                   >  1:2753837/1‑69 (MQ=255)
    ggCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTgc                                                               <  1:1114604/69‑1 (MQ=255)
    ggCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTgc                                                               <  1:3221641/69‑1 (MQ=255)
     gCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTgc                                                               <  1:3430252/68‑1 (MQ=255)
           gCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAgcg                                                        >  1:1339321/1‑69 (MQ=255)
             ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                      >  1:3142824/1‑69 (MQ=255)
             ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                      >  1:3303420/1‑69 (MQ=255)
             ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                      >  1:799519/1‑69 (MQ=255)
             ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                      >  1:1982728/1‑69 (MQ=255)
             ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                      >  1:2809399/1‑69 (MQ=255)
             ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                      >  1:1683570/1‑69 (MQ=255)
               tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:48634/1‑68 (MQ=255)
               tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:1491399/1‑68 (MQ=255)
               tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:3259853/1‑68 (MQ=255)
               tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:3388225/1‑68 (MQ=255)
               tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:2011179/1‑68 (MQ=255)
               tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:971379/1‑68 (MQ=255)
               tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:97426/1‑68 (MQ=255)
               tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:210468/1‑68 (MQ=255)
                ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:1386122/1‑67 (MQ=255)
                ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:282809/1‑67 (MQ=255)
                ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:1470549/1‑67 (MQ=255)
                ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:516525/1‑67 (MQ=255)
                ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:2637388/1‑67 (MQ=255)
                ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                     >  1:2547129/1‑67 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2512579/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2660547/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2826495/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2520722/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:284910/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:1070505/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2880696/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:289473/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2954733/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2994785/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:3063319/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:3402217/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:795639/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:827283/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:956851/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:988601/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:1053493/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2006572/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2420899/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:1375073/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2224989/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2238681/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:1331200/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:1246128/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2319101/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:1118673/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2398030/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:2400174/68‑2 (MQ=255)
                           gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:1082285/68‑2 (MQ=255)
                            cAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                         <  1:1198923/67‑2 (MQ=255)
                                     tCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAgg                               >  1:687168/1‑68 (MQ=255)
                                       atatGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTa                             >  1:2093868/1‑68 (MQ=255)
                                       atatGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTAc                            >  1:3281154/1‑69 (MQ=255)
                                       atatGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTAc                            >  1:2185521/1‑69 (MQ=255)
                                       atatGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTAc                            >  1:2280349/1‑69 (MQ=255)
                                          tGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATc                         >  1:1456229/1‑69 (MQ=255)
                                           gACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGAGATCCGCAAAGTGCATc                         >  1:1828158/1‑68 (MQ=255)
                                           gACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGAGATCCGCAAAGTGCATc                         >  1:1642729/1‑68 (MQ=255)
                                           gACAGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCATCAGCGCCGAGATCCGCAAAGTGCATCAg                       >  1:490343/1‑69 (MQ=255)
                                              cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                     <  1:2461463/69‑1 (MQ=255)
                                              cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                     <  1:2611535/69‑1 (MQ=255)
                                              cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                     <  1:1192918/69‑1 (MQ=255)
                                              cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                     <  1:122172/69‑1 (MQ=255)
                                              cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                     <  1:777118/69‑1 (MQ=255)
                                              cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                     <  1:1307378/69‑1 (MQ=255)
                                              cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                     <  1:3372185/69‑1 (MQ=255)
                                               ggCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                     <  1:1420467/68‑1 (MQ=255)
                                                         tGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCtttt          >  1:1386311/1‑69 (MQ=255)
                                                         tGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCtttt          >  1:2311586/1‑69 (MQ=255)
                                                           ggTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTc                                     >  1:987397/1‑40 (MQ=255)
                                                           ggTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCg        >  1:971663/1‑69 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:2041407/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:1285544/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:527047/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:92415/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:3360409/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:783210/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:2623186/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:1624180/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:3125334/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:914414/69‑1 (MQ=255)
                                                              tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTAGCAAGGTACATCAGTATCTGACCTTTTCGGTg     <  1:1177302/69‑1 (MQ=255)
                                                                 tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGaa   <  1:2403468/68‑1 (MQ=255)
                                                                 tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGaa   <  1:2846773/68‑1 (MQ=255)
                                                                 tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGaa   <  1:2343094/68‑1 (MQ=255)
                                                                 tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGaa   <  1:2847338/68‑1 (MQ=255)
                                                                 tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAAt  >  1:841760/1‑69 (MQ=255)
                                                                 tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAAt  >  1:3148235/1‑69 (MQ=255)
                                                                 tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAAt  >  1:706978/1‑69 (MQ=255)
                                                                   |                                                                  
AGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAAT  >  minE/389159‑389292

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: