Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 389,181 A→G 35.0% S233S (TCA→TCG ybdL → methionine aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE389,1810AG35.0% 37.1 / 33.5 60S233S (TCA→TCGybdLmethionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (27/12);  new base G (7/14);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCC  >  minE/389119‑389241
                                                              |                                                            
caGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGc                                                        <  1:3085969/69‑1 (MQ=255)
caGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGc                                                        <  1:225380/69‑1 (MQ=255)
    gCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGa                                                    >  1:2313353/1‑69 (MQ=255)
      cATGCCAGTGTTCTGGCGCATCCGCAGCTGCGTGAGCGGTCAGTGGCGGTTTCTTCATTTGGCAAGAcc                                                  <  1:1129021/69‑1 (MQ=255)
      cATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGAcc                                                  <  1:3338323/69‑1 (MQ=255)
      cATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGAcc                                                  <  1:338402/69‑1 (MQ=255)
       aTGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGAcc                                                  <  1:2296395/68‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:1439648/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:3175842/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:3098355/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:1265173/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:2823280/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:391184/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:2330389/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:605961/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:1890535/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:1756766/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:325763/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:682184/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTGAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:1500850/69‑1 (MQ=255)
               gtgCTGGCGCATCCGCAATTGCGTCAGCGAGCGGTGGCGGTGTCATCGTTTGGCAAGACCTATCATATg                                         <  1:1932213/69‑1 (MQ=255)
                tgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCAtt                                                            >  1:2859832/1‑49 (MQ=255)
                tgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGAcc                                                  <  1:3340606/59‑1 (MQ=255)
                tgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGAcc                                                  <  1:1471997/59‑1 (MQ=255)
                tgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCat                                            >  1:967852/1‑65 (MQ=255)
                tgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCat                                            <  1:1509782/65‑1 (MQ=255)
                tgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCat                                            >  1:2703159/1‑65 (MQ=255)
                tgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGa                                        >  1:137912/1‑69 (MQ=255)
                tgCTGGCGCATCCGCAGCTGCGTGAGCCGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCat                                            >  1:488446/1‑65 (MQ=255)
                                gCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTgg                         >  1:1607110/1‑68 (MQ=255)
                                gCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTgg                         >  1:1847565/1‑68 (MQ=255)
                                gCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTgg                         >  1:1276176/1‑68 (MQ=255)
                                gCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTgg                         >  1:420974/1‑68 (MQ=255)
                                gCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTgg                         >  1:615155/1‑68 (MQ=255)
                                        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)
                                                     ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTgg                         >  1:3010036/1‑47 (MQ=255)
                                                     ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc   >  1:799519/1‑69 (MQ=255)
                                                     ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc   >  1:2809399/1‑69 (MQ=255)
                                                     ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc   >  1:3142824/1‑69 (MQ=255)
                                                     ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc   >  1:1683570/1‑69 (MQ=255)
                                                     ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc   >  1:1982728/1‑69 (MQ=255)
                                                     ggTGTCATCGTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGCGTtgcgcctgcgcc   >  1:3303420/1‑69 (MQ=255)
                                                       tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:2011179/1‑68 (MQ=255)
                                                       tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:97426/1‑68 (MQ=255)
                                                       tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:971379/1‑68 (MQ=255)
                                                       tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:1491399/1‑68 (MQ=255)
                                                       tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:48634/1‑68 (MQ=255)
                                                       tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:210468/1‑68 (MQ=255)
                                                       tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:3388225/1‑68 (MQ=255)
                                                       tttCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:3259853/1‑68 (MQ=255)
                                                        ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:516525/1‑67 (MQ=255)
                                                        ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:2547129/1‑67 (MQ=255)
                                                        ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:2637388/1‑67 (MQ=255)
                                                        ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:1470549/1‑67 (MQ=255)
                                                        ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:1386122/1‑67 (MQ=255)
                                                        ttcttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGccc  >  1:282809/1‑67 (MQ=255)
                                                              |                                                            
CAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCC  >  minE/389119‑389241

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: