Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,548,960 C→T 13.5% L13L (CTG→TTG)  kbl → glycine C‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,548,9600CT13.5% 165.7 / 15.6 59L13L (CTG→TTG) kblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/21);  new base T (5/3);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

ATGTGCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCGCAGCAAGCAGATAT  >  minE/2548890‑2549030
                                                                      |                                                                      
atGTGCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATc                                                                        >  1:2704255/1‑71 (MQ=255)
atGTGCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAAAGATc                                                                        >  1:2655210/1‑71 (MQ=255)
 tgtgCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCt                                                                       >  1:2971222/1‑71 (MQ=255)
 tgtgCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCt                                                                       >  1:1950573/1‑71 (MQ=255)
 tgtgCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCt                                                                       >  1:3254434/1‑71 (MQ=255)
 tgtgCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCt                                                                       >  1:139842/1‑71 (MQ=255)
 tgtgCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCt                                                                       >  1:2772000/1‑71 (MQ=255)
    gCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTgg                                                                     >  1:2383972/1‑70 (MQ=255)
     cTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGaa                                                                   <  1:1939735/71‑1 (MQ=255)
     cTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGa                                                                    <  1:2406042/70‑1 (MQ=255)
       ggAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTg                                                                      <  1:1507518/66‑1 (MQ=255)
         aaTTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATTTGGAAACCg                                                               <  1:1773210/71‑1 (MQ=255)
             tGCCCTGACTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACg                                                           >  1:3196345/1‑71 (MQ=255)
                    tCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCg                                                               <  1:2212314/60‑1 (MQ=255)
                    tCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGa                                                     <  1:2052152/70‑1 (MQ=255)
                     cTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATTTGGa                                                                    <  1:2494095/54‑1 (MQ=255)
                       ggAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAggg                                                 >  1:932277/1‑71 (MQ=255)
                       ggAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAgg                                                  <  1:2177934/70‑1 (MQ=255)
                       ggAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAgg                                                  <  1:16831/70‑1 (MQ=255)
                        gagaATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCg                                                               <  1:2771855/56‑1 (MQ=255)
                        gagaATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCg                                                               <  1:2704877/56‑1 (MQ=255)
                        gagaATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAggg                                                 <  1:2042279/70‑1 (MQ=255)
                        gagaATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGgt                                                >  1:1643106/1‑71 (MQ=255)
                           aaTCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGaa                                                    >  1:1739854/1‑64 (MQ=255)
                           aaTCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGaa                                                    >  1:903950/1‑64 (MQ=255)
                             tCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGgttgtt                                            >  1:1724850/1‑70 (MQ=255)
                             tCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGgttgtt                                            >  1:68207/1‑70 (MQ=255)
                                  aTGCGTGGAGAATTTTATCAGCAGTTAACCAACGATTTGGAAACCGCACGGGCGGAAGGGTTGTTTAaag                                       >  1:2859208/1‑70 (MQ=255)
                                   tGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAggg                                                 >  1:2091979/1‑59 (MQ=255)
                                   tGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAggg                                                 >  1:1332456/1‑59 (MQ=255)
                                   tGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAaagaa                                     <  1:3354069/71‑1 (MQ=255)
                                   tGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAaagaa                                     <  1:782329/71‑1 (MQ=255)
                                   tGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAaag                                       <  1:2761743/69‑1 (MQ=255)
                                   tGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAaag                                       <  1:97701/69‑1 (MQ=255)
                                    gCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAaag                                       >  1:1637622/1‑68 (MQ=255)
                                       tGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAgc                                 >  1:231704/1‑71 (MQ=255)
                                         gagaATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAgcg                                >  1:2824108/1‑70 (MQ=255)
                                          agaATTTTATCAGCAGTTAACCAACGGTCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCa                              <  1:2527410/71‑1 (MQ=255)
                                          agaATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAgcg                                <  1:2567520/69‑1 (MQ=255)
                                              ttttATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCat                             >  1:3363786/1‑68 (MQ=255)
                                                  aTCAGCAGTTAACCAACGATTTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTAc                       <  1:1131720/70‑1 (MQ=255)
                                                   tCAGCAGTTAACCAACGATTTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCat                             >  1:2712192/1‑63 (MQ=255)
                                                   tCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAgg                                                  <  1:1208119/42‑1 (MQ=255)
                                                    cagcagTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAaag                                       <  1:1363495/52‑1 (MQ=255)
                                                           tAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCgcag             >  1:2428734/1‑71 (MQ=255)
                                                            aaCCAACGATTTGGAAACCGCACGGGCGGAAGGGTTGTTTaaa                                        >  1:1731199/1‑43 (MQ=255)
                                                            aaCCAACGATTTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGt                     >  1:252189/1‑62 (MQ=255)
                                                            aaCCAACGATTTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGt                     >  1:2208913/1‑62 (MQ=255)
                                                            aaCCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCgcag             >  1:1637464/1‑70 (MQ=255)
                                                            aaCCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCgcag             >  1:2440006/1‑70 (MQ=255)
                                                             aCCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCgcagca           <  1:2915963/71‑1 (MQ=255)
                                                               cAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGt                     >  1:559091/1‑59 (MQ=255)
                                                               cAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGt                     >  1:611724/1‑59 (MQ=255)
                                                               cAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGt                     >  1:1158281/1‑59 (MQ=255)
                                                               cAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCgcagca           <  1:545997/69‑1 (MQ=255)
                                                                aaCGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTTCGTCTGCGCagcaag         >  1:1197358/1‑70 (MQ=255)
                                                                aaCGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCGCagcaag         >  1:3336372/1‑70 (MQ=255)
                                                                aaCGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCGCagcaag         >  1:411093/1‑70 (MQ=255)
                                                                 aCGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCGCagcaagca       <  1:1134320/71‑1 (MQ=255)
                                                                 aCGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCGCagcaag         >  1:339436/1‑69 (MQ=255)
                                                                 aCGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCGCagcaag         >  1:1113140/1‑69 (MQ=255)
                                                                      cTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCGCAGCAAGCAGatat  >  1:2168744/1‑71 (MQ=255)
                                                                      |                                                                      
ATGTGCTGGAATTTGCCCTGTCTGGAGAATCGCAATGCGTGGAGAATTTTATCAGCAGTTAACCAACGATCTGGAAACCGCACGGGCGGAAGGGTTGTTTAAAGAAGAGCGCATTATTACGTCTGCGCAGCAAGCAGATAT  >  minE/2548890‑2549030

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: