Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,318 T→C 22.9% K312K (AAA→AAG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,3180TC22.9% 146.3 / 41.9 83K312K (AAA→AAGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (40/24);  new base C (9/10);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGT  >  minE/609250‑609387
                                                                    |                                                                     
gCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGtt                                                                      >  1:1803539/1‑70 (MQ=255)
 cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                                                                    <  1:3545860/71‑1 (MQ=255)
 cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                                                                    <  1:5186223/71‑1 (MQ=255)
  ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                                                   >  1:3512965/1‑71 (MQ=255)
   cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                                                                  <  1:3130415/71‑1 (MQ=255)
   cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                                                                  <  1:2266684/71‑1 (MQ=255)
    aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                                                                  >  1:3455217/1‑70 (MQ=255)
    aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                                  >  1:4891136/1‑70 (MQ=255)
    aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                                  >  1:123839/1‑70 (MQ=255)
    aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                                  >  1:3794498/1‑70 (MQ=255)
    aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                                  >  1:3589419/1‑70 (MQ=255)
    aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                                  >  1:2193704/1‑70 (MQ=255)
    aGATCGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                                                  >  1:4897799/1‑70 (MQ=255)
     gATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                                                >  1:3407185/1‑71 (MQ=255)
      aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                                                               <  1:3607582/71‑1 (MQ=255)
       ttgggTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                                                >  1:1024587/3‑69 (MQ=255)
       tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                                              <  1:3888221/71‑1 (MQ=255)
        cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                                     >  1:722803/1‑63 (MQ=255)
        cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                                                   >  1:5626028/1‑65 (MQ=255)
        cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                                              <  1:5527717/70‑1 (MQ=255)
        cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                                              <  1:2468755/70‑1 (MQ=255)
             cGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtg                                                         >  1:4023666/1‑70 (MQ=255)
              gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                                                               >  1:4516847/1‑63 (MQ=255)
              gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                                                       >  1:3271539/1‑71 (MQ=255)
              gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                                                       >  1:5303577/1‑71 (MQ=255)
                cAGAGAACAGTTCATCGTATTCCTGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                     <  1:5031927/71‑1 (MQ=255)
                cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGt                                                                       <  1:228596/53‑1 (MQ=255)
                cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                     <  1:278814/71‑1 (MQ=255)
                 agagAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAg                                                            >  1:1972603/1‑63 (MQ=255)
                  gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                                    >  1:2931471/1‑69 (MQ=255)
                  gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                                    >  1:878281/1‑69 (MQ=255)
                    gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa                                                 <  1:1885619/71‑1 (MQ=255)
                      aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt                                               >  1:1719359/1‑71 (MQ=255)
                      aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt                                               >  1:4089161/1‑71 (MQ=255)
                       cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                               >  1:2999037/1‑70 (MQ=255)
                          ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                                     <  1:4923374/45‑1 (MQ=255)
                          ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                                     <  1:4537443/45‑1 (MQ=255)
                          ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                                     <  1:64744/45‑1 (MQ=255)
                          ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                                                     <  1:2220146/45‑1 (MQ=255)
                            cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                                                   >  1:5633186/1‑45 (MQ=255)
                            cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                                                   >  1:2829690/1‑45 (MQ=255)
                            cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCGACCATTTCCtac                                         >  1:3684893/1‑69 (MQ=255)
                            cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                         >  1:1981534/1‑71 (MQ=255)
                            cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                         >  1:3693472/1‑71 (MQ=255)
                            cATCGTATTCCGGAGTAAGCAGGAAGCGAACCATACGCAGTTTCa                                                                   >  1:5174622/1‑45 (MQ=255)
                               caaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgc                                       >  1:1560493/4‑70 (MQ=255)
                                gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGctt                                     >  1:4107251/1‑71 (MQ=255)
                                aaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcc                                      <  1:939224/68‑2 (MQ=255)
                                  aTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                                                >  1:4940931/1‑42 (MQ=255)
                                   ttCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGatg                                                              >  1:4958715/1‑41 (MQ=255)
                                   ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTcta                                                    >  1:958715/1‑51 (MQ=255)
                                   ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGcttctt                                  >  1:5283728/1‑71 (MQ=255)
                                    tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                               <  1:5318468/57‑1 (MQ=255)
                                    tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                               <  1:4064226/57‑1 (MQ=255)
                                     ccGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCaa                                                    >  1:1607554/1‑51 (MQ=255)
                                      cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                               <  1:236672/71‑1 (MQ=255)
                                      cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                               <  1:4578514/71‑1 (MQ=255)
                                        gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                          >  1:4909761/3‑58 (MQ=255)
                                        gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                          >  1:1729511/3‑58 (MQ=255)
                                        gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                     <  1:3155037/47‑1 (MQ=255)
                                        gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAttt                                              <  1:4530147/54‑1 (MQ=255)
                                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                     <  1:1205687/60‑3 (MQ=255)
                                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                     <  1:4567928/60‑3 (MQ=255)
                                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                     <  1:2205378/60‑3 (MQ=255)
                                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                     <  1:4251922/60‑3 (MQ=255)
                                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                     <  1:339874/60‑3 (MQ=255)
                                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                     <  1:3122764/60‑3 (MQ=255)
                                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                     <  1:1301603/60‑3 (MQ=255)
                                             cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                             >  1:299515/1‑50 (MQ=255)
                                             cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                             >  1:1787213/1‑50 (MQ=255)
                                                aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                             >  1:4362599/1‑47 (MQ=255)
                                                aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                             >  1:1682521/1‑47 (MQ=255)
                                                aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                             >  1:2102882/1‑47 (MQ=255)
                                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                               <  1:242077/44‑1 (MQ=255)
                                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                               <  1:811417/44‑1 (MQ=255)
                                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                    >  1:1523225/1‑71 (MQ=255)
                                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                    >  1:299446/1‑71 (MQ=255)
                                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                    <  1:3169799/71‑1 (MQ=255)
                                                  gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCt                      >  1:4128862/1‑68 (MQ=255)
                                                  gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                   >  1:5601422/1‑71 (MQ=255)
                                                    aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                    >  1:5283690/1‑68 (MQ=255)
                                                    aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                   >  1:2070553/1‑69 (MQ=255)
                                                       gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATc                       <  1:4481714/62‑1 (MQ=255)
                                                       gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt              <  1:1970934/71‑1 (MQ=255)
                                                       gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt              >  1:990197/1‑71 (MQ=255)
                                                          ccATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                             >  1:3614193/1‑37 (MQ=255)
                                                            aTACGCAGTTTCTTGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                               <  1:30087/49‑1 (MQ=255)
                                                               cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc      <  1:297722/71‑1 (MQ=255)
                                                               cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc      <  1:4365751/71‑1 (MQ=255)
                                                                   gTTTCATGACCAGGTGGTCTACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGt  <  1:5111201/71‑1 (MQ=255)
                                                                    tttCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCttt            <  1:960013/60‑1 (MQ=255)
                                                                    |                                                                     
GCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGT  >  minE/609250‑609387

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