Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,499,305 T→C 60.4% H51H (CAT→CAC ilvN → acetolactate synthase I, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,499,3050TC60.4% 31.8 / 43.9 68H51H (CAT→CACilvNacetolactate synthase I, small subunit
Reads supporting (aligned to +/- strand):  ref base T (11/16);  new base C (17/24);  total (28/40)
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.81e-01

CCGCGCTTTTAACGTTGAAGGCATTCTTTGTCTGCCGATTCAGGACAGCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATAAGCTGGAA  >  minE/2499245‑2499371
                                                            |                                                                  
ccgcGCTTTTAACGTTGAAGGCATTCTTTGTCTGCCGATTCAGGACAGCGACAAAAGCCATATCTGGCt                                                            <  1:2673054/69‑1 (MQ=255)
ccgcGCTTTTAACGTTGAAGGCATTCTTTGTCTGCCGATTCAGGACAGCGACAAAAGCCATATCTGGCt                                                            <  1:2365219/69‑1 (MQ=255)
ccgcGCTTTTAACGTTGAAGGCATTCTTTGTCTGCCGATTCAGGACAGCGACAAAAGCCATATCTGGCt                                                            <  1:1420857/69‑1 (MQ=255)
ccgcGCTTTTAACGTTGAAGGCATTCTTTGTCTGCCGATTCAGGACAGCGACAAAAGCCATATCTGGCt                                                            <  1:1675384/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:1020739/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:2513609/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:2888197/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:2935262/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:2285439/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:2258635/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:2249797/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:3069841/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:3386407/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:160179/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:827241/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:722486/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:1280046/69‑1 (MQ=255)
     cTTTTAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:1431411/69‑1 (MQ=255)
      ttttAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:383060/68‑1 (MQ=255)
      ttttAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:3025838/68‑1 (MQ=255)
      ttttAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:54471/68‑1 (MQ=255)
      ttttAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:1466063/68‑1 (MQ=255)
      ttttAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:1794634/68‑1 (MQ=255)
      ttttAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:1696039/68‑1 (MQ=255)
      ttttAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:3351464/68‑1 (MQ=255)
         tAACGTTTAAGGCATTCTTTGTCTGCCGATTCAGGACAGCGACAAAAGCCATATCTGGCt                                                            <  1:2864696/60‑1 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:3208669/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:3302426/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:3321198/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:3018972/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:637253/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:918133/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:283826/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:2808826/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:945077/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:2611434/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:2605327/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:2736170/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:1189577/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:22770/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:2235071/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:1083857/1‑68 (MQ=255)
         tAACGTTGAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTGGTCa                                                    >  1:1883640/1‑68 (MQ=255)
                gAAGGCATTCTCTGCCTGCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:2979480/58‑1 (MQ=255)
                        tctcTGCCTGACGATTCAGGACAGCGACAAAAGCCACATCTGGCTACCgg                                                       <  1:1462347/50‑1 (MQ=255)
                              tCTGCCGATTCAGGACAGCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTgg                                  >  1:3280902/1‑65 (MQ=255)
                                 gCCGATTCAGGACAGCGACAAAAGCCACATCTGGCTACTgg                                                       <  1:1669747/41‑1 (MQ=255)
                                    gATTCTGGACAGCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAgatga                         <  1:1553303/68‑1 (MQ=255)
                                    gATTCAGGACAGCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAgatga                         <  1:175962/68‑1 (MQ=255)
                                    gATTCAGGACAGCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAgatga                         <  1:2785319/68‑1 (MQ=255)
                                    gATTCAGGACAGCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAgatga                         <  1:1095656/68‑1 (MQ=255)
                                    gATTCAGGACAGCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAgatga                         <  1:1271585/68‑1 (MQ=255)
                                    gATTCAGGACAGCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAgatga                         <  1:855664/68‑1 (MQ=255)
                                               gCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGc                                        >  1:232922/1‑42 (MQ=255)
                                               gCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATc              >  1:1692380/1‑68 (MQ=255)
                                               gCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATc              >  1:3246365/1‑68 (MQ=255)
                                               gCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATc              >  1:2222860/1‑68 (MQ=255)
                                               gCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATc              >  1:3398034/1‑68 (MQ=255)
                                               gCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATc              >  1:2382471/1‑68 (MQ=255)
                                               gCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATc              >  1:2588256/1‑68 (MQ=255)
                                                   cAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATaa         >  1:1063195/1‑69 (MQ=255)
                                                   cAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATaa         >  1:566993/1‑69 (MQ=255)
                                                   cAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATaa         >  1:452707/1‑69 (MQ=255)
                                                          cATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATAAGCTGGaa  <  1:744977/69‑1 (MQ=255)
                                                          cATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATAAGCTGGaa  <  1:3368991/69‑1 (MQ=255)
                                                          cATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATAAGCTGGaa  <  1:3321115/69‑1 (MQ=255)
                                                          cATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATAAGCTGGaa  <  1:31617/69‑1 (MQ=255)
                                                          cATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATAAGCTGGaa  <  1:2504225/69‑1 (MQ=255)
                                                            |                                                                  
CCGCGCTTTTAACGTTGAAGGCATTCTTTGTCTGCCGATTCAGGACAGCGACAAAAGCCATATCTGGCTACTGGTCAATGACGACCAGCGTCTGGAGCAGATGATAAGCCAAATCGATAAGCTGGAA  >  minE/2499245‑2499371

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: