Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,425,807 T→C 17.0% A95A (GCA→GCG rffG ← dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,425,8070TC17.0% 183.5 / 27.4 82A95A (GCA→GCGrffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base T (32/36);  new base C (4/10);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGAT  >  minE/2425743‑2425874
                                                                |                                                                   
tCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCggg                                                               >  1:820647/1‑71 (MQ=255)
tCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCggg                                                               >  1:5055634/1‑71 (MQ=255)
tCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCggg                                                               >  1:1730489/1‑71 (MQ=255)
 ccAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCggg                                                               <  1:5008626/70‑1 (MQ=255)
    gTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGc                                                                   <  1:5234880/63‑1 (MQ=255)
      aaGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTc                                                          <  1:1561731/70‑1 (MQ=255)
      aaGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCCTc                                                          >  1:4848317/1‑70 (MQ=255)
       aGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                                        <  1:5107199/71‑1 (MQ=255)
       aGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCa                                                         <  1:1187261/70‑1 (MQ=255)
          ccGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAg                                                     <  1:3373636/71‑1 (MQ=255)
           cgcgCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAg                                                     <  1:1520510/70‑1 (MQ=255)
            gcgcCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGaa                                                   >  1:2461800/1‑71 (MQ=255)
            gcgcCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGaa                                                   >  1:672439/1‑71 (MQ=255)
                  cTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGt                                                           <  1:4145212/57‑1 (MQ=255)
                  cTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                                              >  1:4552253/1‑70 (MQ=255)
                  cTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTc                                              >  1:2357644/1‑70 (MQ=255)
                     cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAAc                                           <  1:2728608/70‑1 (MQ=255)
                     cAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAAc                                           <  1:4946705/70‑1 (MQ=255)
                           aTGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTc                                                          <  1:4551543/49‑1 (MQ=255)
                            tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCaa                                                        <  1:3042882/50‑1 (MQ=255)
                            tGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                                   >  1:2624697/1‑71 (MQ=255)
                              tataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGt                                                           >  1:1200551/1‑45 (MQ=255)
                              tataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCt                                     >  1:977172/1‑67 (MQ=255)
                              tataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCg                                 >  1:1051877/1‑71 (MQ=255)
                              tataAGTCCCGACAATGTAGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCt                                     >  1:2583339/1‑67 (MQ=255)
                               ataAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCgg                                >  1:2868128/1‑71 (MQ=255)
                                taAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATg                                        <  1:2085425/62‑1 (MQ=255)
                                 aaGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCaa                                            <  1:2888384/57‑1 (MQ=255)
                                 aaGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAAc                                           <  1:827688/58‑1 (MQ=255)
                                  aGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTg                             >  1:1189486/1‑71 (MQ=255)
                                  aGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGcag                             >  1:3274020/1‑69 (MQ=255)
                                  aGTACCAACAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGca                              <  1:2768039/70‑2 (MQ=255)
                                  aGTACCAACAATATTGGTTTCAATAAATGCCGCACGGCCTGTAATTGAACGGTCAACATGGCTTTCAGca                              <  1:1468395/70‑2 (MQ=255)
                                    tCCCGACAATGTTGGTTTCAATAAACGCTGCCggg                                                               >  1:87077/1‑35 (MQ=255)
                                    tCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc                            >  1:2915128/1‑70 (MQ=255)
                                    tCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGc                            >  1:2165144/1‑70 (MQ=255)
                                      ccaaCAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAg                         <  1:4016703/68‑1 (MQ=255)
                                      ccGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAg                         <  1:1554660/71‑1 (MQ=255)
                                        gACAATGTTGGTTTCAATGAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGa                        <  1:2008773/70‑1 (MQ=255)
                                        gACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGAt                       >  1:2648124/1‑71 (MQ=255)
                                         aCAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGcc                           <  1:423368/66‑1 (MQ=255)
                                         aCAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGAt                       >  1:4480877/1‑70 (MQ=255)
                                          cAATATTGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAGGTGc                     <  1:668410/71‑1 (MQ=255)
                                           aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCttt                                   >  1:4805586/1‑56 (MQ=255)
                                           aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCa                    >  1:3004546/1‑71 (MQ=255)
                                           aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCa                    >  1:2081447/1‑71 (MQ=255)
                                           aaTGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGGCAACATGGCTTTCGGCTGCCAGATGCa                    <  1:4994585/71‑1 (MQ=255)
                                             tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAt                   >  1:1766674/1‑70 (MQ=255)
                                             tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATg                  <  1:1075914/71‑1 (MQ=255)
                                             tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATg                  <  1:1306547/71‑1 (MQ=255)
                                             tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATg                  <  1:1381664/71‑1 (MQ=255)
                                             tGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATg                  <  1:5304892/71‑1 (MQ=255)
                                               ttGGTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAGGTGCATCac                >  1:4258597/1‑71 (MQ=255)
                                               ttGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGc                               >  1:199525/1‑56 (MQ=255)
                                               ttGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCa                          <  1:2187802/61‑1 (MQ=255)
                                                 ggTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCa                    >  1:5260986/1‑65 (MQ=255)
                                                  gTTTCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGcag                             >  1:834593/1‑53 (MQ=255)
                                                   tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCt                                     <  1:2382549/46‑1 (MQ=255)
                                                   tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCt                                     <  1:4754659/46‑1 (MQ=255)
                                                   tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAg            <  1:4314033/71‑1 (MQ=255)
                                                   tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAg            <  1:846269/71‑1 (MQ=255)
                                                   tttCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAg            <  1:3643897/71‑1 (MQ=255)
                                                     tCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAgg                        <  1:4306628/57‑2 (MQ=255)
                                                     tCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAgg                        <  1:4566631/57‑2 (MQ=255)
                                                     tCAATAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAg                         >  1:5231748/1‑56 (MQ=255)
                                                     tCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACa                                          >  1:3626959/1‑39 (MQ=255)
                                                      cAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt         >  1:4457491/1‑71 (MQ=255)
                                                       aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt         <  1:4878499/70‑1 (MQ=255)
                                                       aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt         <  1:3736104/70‑1 (MQ=255)
                                                       aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt         <  1:3665213/70‑1 (MQ=255)
                                                       aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCt         <  1:2843927/70‑1 (MQ=255)
                                                       aaTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTg        >  1:4198487/1‑71 (MQ=255)
                                                        aTAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGcag                             <  1:3045889/49‑3 (MQ=255)
                                                        aTAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTgg       <  1:630480/71‑1 (MQ=255)
                                                         tAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAggt                       <  1:2181131/54‑3 (MQ=12)
                                                         tAAATGCCGCAGGGCCTGTAATTGAACGGTCAACATGGCTTTCAGCAGCCAggt                       <  1:4471808/54‑3 (MQ=12)
                                                          aaaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCa                          >  1:1235619/1‑50 (MQ=255)
                                                          aaaCGCTGCCGGGCCGTCAATAGAACGGTAAACATGGCTTTCGGCTGCCa                          >  1:5329179/1‑50 (MQ=255)
                                                           aaCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCt     <  1:1390005/70‑1 (MQ=255)
                                                            aCGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGa   <  1:2618622/71‑1 (MQ=255)
                                                             cGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGAt  <  1:1752427/71‑1 (MQ=255)
                                                             cGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGAt  >  1:2204662/1‑71 (MQ=255)
                                                                tGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAt                   <  1:2851498/51‑1 (MQ=255)
                                                                tGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCAt                   <  1:2917800/51‑1 (MQ=255)
                                                                |                                                                   
TCCAGTAAGCCCGCGCCGCTTCAAGCAATGTATAAGTCCCGACAATGTTGGTTTCAATAAACGCTGCCGGGCCGTCAATAGAACGGTCAACATGGCTTTCGGCTGCCAGATGCATGACACAGTCTGGCTGAT  >  minE/2425743‑2425874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: