Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,282 C→T 17.1% A658T (GCT→ACT) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,2820CT17.1% 256.6 / 38.4 111A658T (GCT→ACT) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (48/44);  new base T (16/3);  total (64/47)
Fisher's exact test for biased strand distribution p-value = 1.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 111 total aligned reads displayed.

CTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGT  >  minE/608213‑608350
                                                                     |                                                                    
cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:3577571/71‑1 (MQ=255)
cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:2603901/71‑1 (MQ=255)
cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:3001695/71‑1 (MQ=255)
cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:601071/71‑1 (MQ=255)
cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:3063704/71‑1 (MQ=255)
cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:4424500/71‑1 (MQ=255)
   cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                                   <  1:747403/70‑1 (MQ=255)
    gtcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACg                                                                 >  1:4515711/1‑71 (MQ=255)
     tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACgg                                                                >  1:384793/1‑71 (MQ=255)
     tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACgg                                                                >  1:1411211/1‑71 (MQ=255)
     tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTATAGCGTAAGCAAACgg                                                                <  1:1381964/71‑1 (MQ=255)
            tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                                   >  1:2424067/1‑61 (MQ=255)
            tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                                   >  1:2353965/1‑61 (MQ=255)
             aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGt                                                            <  1:1894733/67‑1 (MQ=255)
             aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                         <  1:5430160/70‑1 (MQ=255)
               ccAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGc                                                               <  1:4068465/62‑1 (MQ=255)
                  gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                   >  1:5103641/1‑71 (MQ=255)
                   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                  <  1:666583/71‑1 (MQ=255)
                   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                   >  1:2595823/1‑70 (MQ=255)
                   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                   >  1:3162435/1‑70 (MQ=255)
                     cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                  >  1:726501/1‑69 (MQ=255)
                     cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                <  1:184186/71‑1 (MQ=255)
                          ggAACGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                            >  1:3755348/1‑70 (MQ=255)
                          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           <  1:5653632/71‑1 (MQ=255)
                          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           <  1:1390292/71‑1 (MQ=255)
                             aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCa                                                              >  1:3865439/1‑49 (MQ=255)
                             aCGATTGAGAAGGTGTACGAGATCCCATATTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                          >  1:4234740/1‑67 (MQ=255)
                              cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                                             >  1:4943854/1‑49 (MQ=255)
                              cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                                             >  1:2378580/1‑49 (MQ=255)
                              cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                                             >  1:4597779/1‑49 (MQ=255)
                               gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                <  1:4060079/61‑1 (MQ=255)
                               gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                <  1:5255274/61‑1 (MQ=255)
                                aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     <  1:1327759/71‑1 (MQ=255)
                                 ttgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                     >  1:5678493/3‑70 (MQ=255)
                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:2004201/38‑1 (MQ=255)
                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     <  1:4135293/70‑1 (MQ=255)
                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     >  1:1762916/1‑70 (MQ=255)
                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     >  1:3348759/1‑70 (MQ=255)
                                  agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                   >  1:5573347/1‑71 (MQ=255)
                                  agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                   >  1:2969125/1‑71 (MQ=255)
                                  agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                   >  1:4378990/1‑71 (MQ=255)
                                   gagaAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                                   >  1:2964523/1‑38 (MQ=255)
                                   gagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGc                                                       <  1:4579883/50‑1 (MQ=255)
                                     gaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                <  1:920718/71‑1 (MQ=255)
                                      aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                               <  1:3033685/71‑1 (MQ=255)
                                      aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                               <  1:3026919/71‑1 (MQ=255)
                                        ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                    >  1:3370390/1‑48 (MQ=255)
                                        ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                    >  1:1242477/1‑48 (MQ=255)
                                        ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                    >  1:1318659/1‑48 (MQ=255)
                                        ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                    >  1:3013444/1‑48 (MQ=255)
                                         gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCAcc                            >  1:397857/1‑71 (MQ=255)
                                          tgtAGTAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:4083921/39‑1 (MQ=255)
                                          tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                            <  1:2288181/54‑1 (MQ=255)
                                          tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                            <  1:5068583/54‑1 (MQ=255)
                                           gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                          >  1:267150/1‑71 (MQ=255)
                                            tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                          >  1:324205/1‑54 (MQ=255)
                                             aGGAGATACCATCTTAAGCGAAAGCAAACGGCAGTTTAGc                                                       <  1:177471/40‑1 (MQ=38)
                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:1787219/35‑1 (MQ=255)
                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:3020729/35‑1 (MQ=255)
                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:1009061/35‑1 (MQ=255)
                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:5603934/35‑1 (MQ=255)
                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                         >  1:1595621/1‑37 (MQ=255)
                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                        >  1:5311233/1‑54 (MQ=255)
                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                       >  1:415944/1‑71 (MQ=255)
                                              cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGttt                                                          >  1:3768728/2‑36 (MQ=12)
                                              cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGttt                                                          >  1:1605421/2‑36 (MQ=12)
                                              cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                            <  1:813430/49‑1 (MQ=255)
                                               gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                      >  1:4114886/1‑68 (MQ=255)
                                               gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                      >  1:1693457/1‑68 (MQ=255)
                                               gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                      <  1:834402/71‑4 (MQ=255)
                                               gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                   >  1:4826911/1‑58 (MQ=255)
                                               gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                   >  1:3596066/1‑58 (MQ=255)
                                                  aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                              <  1:3959556/44‑1 (MQ=255)
                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                     >  1:498555/1‑35 (MQ=255)
                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                     >  1:2090346/1‑35 (MQ=255)
                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                     >  1:2515299/1‑35 (MQ=255)
                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGa                                                 >  1:585966/1‑39 (MQ=255)
                                                     ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                 >  1:3762863/1‑70 (MQ=255)
                                                     ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                 >  1:3321072/1‑70 (MQ=255)
                                                     ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                     <  1:1633976/66‑1 (MQ=255)
                                                     ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                 >  1:526394/1‑70 (MQ=255)
                                                      cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa               >  1:489970/1‑71 (MQ=255)
                                                       aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc              <  1:1258203/71‑1 (MQ=255)
                                                        tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGATCACGAc              >  1:432341/1‑70 (MQ=255)
                                                         cTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAcca            >  1:5093733/1‑70 (MQ=255)
                                                          tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                          <  1:1066885/40‑1 (MQ=255)
                                                          tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                          >  1:4230043/1‑40 (MQ=255)
                                                          tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           >  1:3917764/1‑39 (MQ=255)
                                                           ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt           >  1:3836571/1‑70 (MQ=255)
                                                           ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg          >  1:1192396/1‑71 (MQ=255)
                                                           ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCCTg          >  1:4201025/1‑71 (MQ=255)
                                                            tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc             <  1:635675/67‑1 (MQ=255)
                                                            tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt           <  1:2513253/69‑1 (MQ=255)
                                                            tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg          >  1:3385579/1‑70 (MQ=255)
                                                               cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc      <  1:3595044/71‑1 (MQ=255)
                                                                gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg     >  1:5471197/1‑71 (MQ=255)
                                                                 tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg     <  1:2509113/70‑1 (MQ=255)
                                                                   aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc              >  1:3140718/1‑59 (MQ=255)
                                                                   aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc              >  1:4945308/1‑59 (MQ=255)
                                                                   aGCAAACGGCAGTTTAGAAACGGAAGTCAGAGAGGCTACTCCACCTTTCTGGTAACGACCGTGCATCGGGt  >  1:4693906/1‑71 (MQ=255)
                                                                     |                                                                    
CTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGT  >  minE/608213‑608350

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: