Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 389,235 A→T 37.0% P251P (CCA→CCT ybdL → methionine aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE389,2350AT37.0% 68.7 / 74.6 89P251P (CCA→CCTybdLmethionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base T (3/30);  total (31/58)
Fisher's exact test for biased strand distribution p-value = 7.04e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

GCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCAC  >  minE/389170‑389302
                                                                 |                                                                   
gCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAgcg                                                                  >  1:1339321/1‑69 (MQ=255)
  ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                                >  1:3142824/1‑69 (MQ=255)
  ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                                >  1:3303420/1‑69 (MQ=255)
  ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                                >  1:1982728/1‑69 (MQ=255)
  ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                                >  1:799519/1‑69 (MQ=255)
  ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                                >  1:1683570/1‑69 (MQ=255)
  ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc                                                                >  1:2809399/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:210468/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)
     ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                               >  1:1386122/1‑67 (MQ=255)
     ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                               >  1:2637388/1‑67 (MQ=255)
     ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                               >  1:282809/1‑67 (MQ=255)
     ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                               >  1:1470549/1‑67 (MQ=255)
     ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                               >  1:2547129/1‑67 (MQ=255)
     ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc                                                               >  1:516525/1‑67 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2994785/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2660547/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:1053493/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2512579/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2420899/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2826495/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2400174/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2398030/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:284910/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:2520722/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:1070505/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:1331200/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:1375073/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:1246128/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2006572/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:1118673/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2224989/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2238681/68‑2 (MQ=255)
                gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTTGCGCCTGCGCCCATCAGCGCCGag                                                   <  1:2319101/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:2280349/1‑69 (MQ=255)
                            atatGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTAc                                      >  1:3281154/1‑69 (MQ=255)
                            atatGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTAc                                      >  1:2185521/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:122172/69‑1 (MQ=255)
                                   cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                               <  1:2611535/69‑1 (MQ=255)
                                   cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                               <  1:777118/69‑1 (MQ=255)
                                   cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                               <  1:1192918/69‑1 (MQ=255)
                                   cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                               <  1:1307378/69‑1 (MQ=255)
                                   cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                               <  1:3372185/69‑1 (MQ=255)
                                   cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                               <  1:2461463/69‑1 (MQ=255)
                                    ggCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTa                               <  1:1420467/68‑1 (MQ=255)
                                              tGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCtttt                    >  1:2311586/1‑69 (MQ=255)
                                              tGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCtttt                    >  1:1386311/1‑69 (MQ=255)
                                                ggTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTc                                               >  1:987397/1‑40 (MQ=255)
                                                ggTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCg                  >  1:971663/1‑69 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:1285544/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:527047/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:783210/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:3360409/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:3125334/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:914414/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:1624180/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:92415/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:2041407/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:2623186/69‑1 (MQ=255)
                                                   tATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTAGCAAGGTACATCAGTATCTGACCTTTTCGGTg               <  1:1177302/69‑1 (MQ=255)
                                                      tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGaa             <  1:2343094/68‑1 (MQ=255)
                                                      tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGaa             <  1:2847338/68‑1 (MQ=255)
                                                      tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGaa             <  1:2846773/68‑1 (MQ=255)
                                                      tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGaa             <  1:2403468/68‑1 (MQ=255)
                                                      tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAAt            >  1:706978/1‑69 (MQ=255)
                                                      tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAAt            >  1:3148235/1‑69 (MQ=255)
                                                      tgtgTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAAt            >  1:841760/1‑69 (MQ=255)
                                                              gcCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCgg     <  1:1528674/68‑1 (MQ=255)
                                                              gcCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCgg     <  1:2162635/68‑1 (MQ=255)
                                                               ccAGCGCCCATCAGCGCCAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGcac  <  1:2387007/69‑1 (MQ=255)
                                                                cAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGcac  <  1:2287666/69‑1 (MQ=255)
                                                                cAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGcac  <  1:2434807/69‑1 (MQ=255)
                                                                 aGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGcac  <  1:3193055/68‑1 (MQ=255)
                                                                 |                                                                   
GCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCAC  >  minE/389170‑389302

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: