Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,072 G→T 24.9% Q728K (CAG→AAG) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,0720GT24.9% 158.5 / 29.9 64Q728K (CAG→AAG) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (27/21);  new base T (7/9);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

GTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAG  >  minE/608004‑608138
                                                                    |                                                                  
gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                                  >  1:5527442/1‑71 (MQ=255)
gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                                  >  1:3563115/1‑71 (MQ=255)
   ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                                >  1:2966344/1‑70 (MQ=255)
   ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                                >  1:3375138/1‑70 (MQ=255)
   ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                                >  1:1178530/1‑70 (MQ=255)
   ctgctgcagTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                                >  1:562899/1‑70 (MQ=255)
         ctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                                                         <  1:1119205/71‑1 (MQ=255)
            ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                                       <  1:1488244/70‑1 (MQ=255)
             tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                                     >  1:5385737/1‑71 (MQ=255)
             tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                                     >  1:2144533/1‑71 (MQ=255)
             tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                                     >  1:3291109/1‑71 (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)
                        cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATcg                                          >  1:2599153/1‑71 (MQ=255)
                          cgTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTcc                                                       <  1:5500462/55‑1 (MQ=12)
                            ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                                                        <  1:1169838/53‑1 (MQ=255)
                            ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                                                        <  1:1855915/53‑1 (MQ=255)
                              gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCa                                             >  1:3090896/1‑62 (MQ=255)
                              gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                     <  1:1841906/70‑1 (MQ=255)
                                aaCGTACTGCGTAGCCAGATACAGGGATGGTCAGCTGCGGATATTTTTcc                                                       <  1:4784597/50‑1 (MQ=255)
                                 aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                                 <  1:2574704/71‑1 (MQ=255)
                                 aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                   >  1:4723853/1‑69 (MQ=255)
                                    tACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttt                                                           <  1:1550897/42‑1 (MQ=255)
                                      cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                            >  1:2496566/3‑71 (MQ=12)
                                      cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAt                             >  1:2145640/3‑70 (MQ=12)
                                        gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                                       <  1:4827450/42‑1 (MQ=255)
                                        gCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTc                          >  1:3894100/1‑71 (MQ=16)
                                         cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                                            >  1:4987806/1‑52 (MQ=255)
                                         cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                           <  1:2534990/69‑1 (MQ=255)
                                         cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                         >  1:686846/1‑71 (MQ=255)
                                            aGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga                      >  1:611247/1‑70 (MQ=12)
                                             gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCgat                     >  1:204410/1‑69 (MQ=12)
                                                agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtc                   <  1:3476824/70‑1 (MQ=9)
                                                aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                  >  1:4707475/1‑71 (MQ=255)
                                                 gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                 <  1:3047823/68‑1 (MQ=9)
                                                 gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                 <  1:1208217/68‑1 (MQ=9)
                                                 gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                 <  1:1890052/68‑1 (MQ=9)
                                                   tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                <  1:3986411/70‑1 (MQ=255)
                                                   tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                <  1:4459557/70‑1 (MQ=255)
                                                   gacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcatta               >  1:2583400/2‑69 (MQ=9)
                                                    acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtc                   >  1:2752349/1‑66 (MQ=255)
                                                    acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc              >  1:125761/1‑71 (MQ=255)
                                                    acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc              >  1:378263/1‑71 (MQ=255)
                                                     cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg             <  1:5515929/71‑1 (MQ=12)
                                                     cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg             <  1:2610403/71‑1 (MQ=12)
                                                     cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg             <  1:2878476/71‑1 (MQ=12)
                                                     cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg             <  1:3994210/71‑1 (MQ=12)
                                                     cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACg             >  1:1900190/1‑71 (MQ=255)
                                                       cGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt           <  1:724340/71‑1 (MQ=255)
                                                        ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt           >  1:2256425/1‑70 (MQ=255)
                                                         gATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa         >  1:3215964/1‑71 (MQ=255)
                                                             gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa         <  1:388756/67‑1 (MQ=255)
                                                             gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGt      <  1:623050/70‑1 (MQ=255)
                                                               cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                >  1:3268764/1‑58 (MQ=255)
                                                               cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt     >  1:672858/1‑69 (MQ=255)
                                                                aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa         >  1:5346504/1‑64 (MQ=255)
                                                                aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg  >  1:4238926/1‑71 (MQ=255)
                                                                    |                                                                  
GTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAG  >  minE/608004‑608138

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: