Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,030 A→C 9.6% S742A (TCG→GCG) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,0300AC 9.6% 225.7 / 18.3 94S742A (TCG→GCG) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (56/29);  new base C (4/5);  total (60/34)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

AAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCG  >  minE/607963‑608098
                                                                    |                                                                    
aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                                   >  1:973397/1‑71 (MQ=255)
aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                                   >  1:4352840/1‑71 (MQ=255)
aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                                   <  1:645509/71‑1 (MQ=255)
aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                                   >  1:2591745/1‑71 (MQ=255)
aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                                   >  1:229397/1‑71 (MQ=255)
aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGtt                                                                   <  1:2151974/71‑1 (MQ=255)
   tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGa                                                                      <  1:1309261/65‑1 (MQ=255)
    cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                                               >  1:982748/1‑71 (MQ=255)
    cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                                               >  1:4477486/1‑71 (MQ=255)
    cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                                               >  1:2170650/1‑71 (MQ=255)
       aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGaaa                                                               <  1:4419698/68‑1 (MQ=255)
       aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                             >  1:2298041/1‑70 (MQ=255)
       aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                             <  1:1962975/70‑1 (MQ=255)
       aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                            <  1:3026259/71‑1 (MQ=255)
        ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                            >  1:4396549/1‑70 (MQ=255)
        ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                           >  1:2300988/1‑71 (MQ=255)
        ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGATGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                            >  1:3494523/1‑70 (MQ=255)
           cATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACt                                                         <  1:4815956/70‑1 (MQ=255)
             tAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGGCAGCGAGTTGAAACGTACTGc                                                       <  1:3749529/70‑1 (MQ=255)
              aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                                    >  1:3197352/1‑56 (MQ=255)
              aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                                    >  1:71531/1‑56 (MQ=255)
              aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                                    >  1:2244981/1‑56 (MQ=255)
              aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                                    >  1:210154/1‑56 (MQ=255)
              aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                                    >  1:176854/1‑56 (MQ=255)
              aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                                    >  1:2338224/1‑56 (MQ=255)
              aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGt                                                                    >  1:2561674/1‑56 (MQ=255)
              aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCg                                                      >  1:3536577/1‑70 (MQ=255)
                 ttGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                                                  <  1:2617571/71‑1 (MQ=255)
                  tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                             >  1:1801633/1‑59 (MQ=255)
                  tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                             >  1:1953076/1‑59 (MQ=255)
                   gAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGAGAGTTGAAACGTACTGCGTAGcc                                                 <  1:83408/70‑1 (MQ=255)
                     gtaaaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGa                                              <  1:350587/68‑1 (MQ=255)
                      tGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGt                                                            <  1:1294776/56‑1 (MQ=255)
                        aaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                                            <  1:404675/70‑3 (MQ=255)
                        aaGGTACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                                            <  1:2969174/70‑3 (MQ=255)
                        aaGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                            >  1:5315911/1‑70 (MQ=255)
                        aaGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                            >  1:227011/1‑70 (MQ=255)
                         aGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                             >  1:805418/1‑52 (MQ=255)
                         aGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACg                                                             >  1:5010636/1‑52 (MQ=255)
                            tACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac                                         >  1:3279435/1‑69 (MQ=255)
                            tACGTGAAATAACATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                        >  1:398507/1‑70 (MQ=255)
                                aGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACTGATgg                                   <  1:87811/71‑1 (MQ=255)
                                aGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg                                   <  1:3814349/71‑1 (MQ=255)
                                  taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                       <  1:3326983/65‑1 (MQ=255)
                                  taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                       <  1:5158368/65‑1 (MQ=255)
                                  taataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                                        <  1:366473/64‑1 (MQ=255)
                                   aataaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGCATTGTCa                                <  1:569736/71‑1 (MQ=255)
                                    ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                     >  1:3593189/1‑65 (MQ=255)
                                    ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                     >  1:2715948/1‑65 (MQ=255)
                                    ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                     >  1:531529/1‑65 (MQ=255)
                                    ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                     >  1:5538430/1‑65 (MQ=255)
                                    ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                               <  1:472346/71‑1 (MQ=255)
                                    ataaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                               >  1:1064665/1‑71 (MQ=255)
                                    ataaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAg                               >  1:4040379/1‑71 (MQ=255)
                                     taaCATCCTGTTGCT‑GTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                <  1:3450180/69‑1 (MQ=255)
                                      aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                           >  1:2849456/1‑41 (MQ=255)
                                      aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                           >  1:3235411/1‑41 (MQ=255)
                                      aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTa                                                           >  1:739180/1‑41 (MQ=255)
                                      aaCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAg                                                   <  1:1961296/49‑1 (MQ=255)
                                       aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGt                                                     >  1:1115134/1‑46 (MQ=255)
                                       aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                       >  1:2975706/1‑60 (MQ=255)
                                       aCGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                >  1:2340301/1‑67 (MQ=255)
                                         gTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                          >  1:3563115/1‑71 (MQ=255)
                                         gTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                          >  1:5527442/1‑71 (MQ=255)
                                          tCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                  <  1:5225139/62‑1 (MQ=255)
                                           ccTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTAc                                                          >  1:2986568/1‑37 (MQ=255)
                                            ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                >  1:753132/1‑62 (MQ=255)
                                            ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                        >  1:2966344/1‑70 (MQ=255)
                                            ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                        >  1:1178530/1‑70 (MQ=255)
                                            ctgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                        >  1:3375138/1‑70 (MQ=255)
                                            ctgctgca‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                        >  1:562899/1‑70 (MQ=255)
                                             tgctgct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                             <  1:1232043/64‑1 (MQ=255)
                                               ctgctgGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                                       >  1:4867945/1‑37 (MQ=38)
                                                  ct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                                 <  1:1052531/39‑1 (MQ=255)
                                                  ct‑gTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                 <  1:1119205/71‑1 (MQ=255)
                                                      ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc               <  1:1488244/70‑1 (MQ=255)
                                                       tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg             >  1:3291109/1‑71 (MQ=255)
                                                       tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg             >  1:2144533/1‑71 (MQ=255)
                                                       tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg             >  1:5385737/1‑71 (MQ=255)
                                                        cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                              >  1:5045808/1‑37 (MQ=255)
                                                        cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                              >  1:4501310/1‑37 (MQ=255)
                                                        cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                            >  1:2000940/1‑39 (MQ=255)
                                                        cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                            >  1:2851405/1‑39 (MQ=255)
                                                        cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                               >  1:4205704/1‑52 (MQ=255)
                                                            aGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                 <  1:5332446/62‑1 (MQ=255)
                                                             gTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCgg             >  1:3904211/1‑65 (MQ=9)
                                                             gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc       >  1:4713787/1‑71 (MQ=255)
                                                             gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc       <  1:2381112/71‑1 (MQ=255)
                                                              tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc      <  1:4341197/71‑1 (MQ=255)
                                                              tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc      <  1:1067665/71‑1 (MQ=255)
                                                                aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt          <  1:5478999/65‑1 (MQ=255)
                                                                aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt          <  1:862411/65‑1 (MQ=255)
                                                                aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt    >  1:4465790/1‑71 (MQ=255)
                                                                 gCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                  <  1:2426693/40‑1 (MQ=255)
                                                                  cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATcg  >  1:2599153/1‑71 (MQ=255)
                                                                    |                                                                    
AAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCT‑GTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCG  >  minE/607963‑608098

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: ‑