Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,019 T→C 17.6% N412S (AAC→AGC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,0190TC17.6% 268.9 / 47.5 125N412S (AAC→AGC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (43/60);  new base C (12/10);  total (55/70)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 150 total aligned reads displayed.

CGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAG  >  minE/608949‑609088
                                                                      |                                                                     
cgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgt                                                                       <  1:1980355/71‑1 (MQ=255)
 gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                                      >  1:3664825/1‑71 (MQ=255)
 gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                                      >  1:1677563/1‑71 (MQ=255)
  cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                                                     <  1:2154030/71‑1 (MQ=255)
  cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctg                                                                     >  1:1946567/1‑71 (MQ=255)
  cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctg                                                                     >  1:3185884/1‑71 (MQ=255)
     caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctgtt                                                                   >  1:514048/3‑70 (MQ=255)
     caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                                  <  1:5360228/69‑1 (MQ=255)
     caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                                  <  1:3782391/69‑1 (MQ=255)
     cGAAGATCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                                  <  1:2340665/71‑1 (MQ=255)
     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                                  <  1:4459885/71‑1 (MQ=255)
     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                                  <  1:3985486/71‑1 (MQ=255)
      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGa                                                                 <  1:217997/71‑1 (MQ=255)
      aaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                                  >  1:1616873/2‑70 (MQ=255)
      aaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                                  >  1:4544652/2‑70 (MQ=255)
       aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                                  >  1:2232396/1‑69 (MQ=255)
         gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                                                              >  1:24494/1‑71 (MQ=255)
          aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                                                             <  1:1569837/71‑1 (MQ=255)
                aTTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACg                                                       <  1:1878427/71‑1 (MQ=255)
                  ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                                           <  1:5176578/65‑1 (MQ=255)
                  ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                                           <  1:659474/65‑1 (MQ=255)
                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                                    >  1:268613/1‑71 (MQ=255)
                        aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                                                <  1:3253912/70‑1 (MQ=255)
                           aaCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                            <  1:2096741/71‑1 (MQ=255)
                              gATCATCGGTCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                                         >  1:3346541/1‑71 (MQ=255)
                                tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                              >  1:4006540/1‑64 (MQ=255)
                                tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                              >  1:1007848/1‑64 (MQ=255)
                                tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                              >  1:947131/1‑64 (MQ=255)
                                tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                              >  1:1557215/1‑64 (MQ=255)
                                tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                       <  1:3016969/71‑1 (MQ=255)
                                 catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                      >  1:5542797/1‑71 (MQ=255)
                                 catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                      >  1:1482523/1‑71 (MQ=255)
                                 catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                      >  1:4750025/1‑71 (MQ=255)
                                  atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                                                 <  1:2042020/59‑1 (MQ=255)
                                   ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAg                                                 <  1:5126323/56‑1 (MQ=255)
                                   ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                    >  1:5684195/3‑71 (MQ=255)
                                   ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                    >  1:4556021/3‑71 (MQ=255)
                                   tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                                           <  1:2755966/48‑1 (MQ=255)
                                     gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                                   <  1:4952891/70‑1 (MQ=255)
                                     gggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTca                                                            >  1:642556/1‑44 (MQ=255)
                                     gggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTca                                                            >  1:1596959/1‑44 (MQ=255)
                                         ttACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                               <  1:2953341/70‑1 (MQ=255)
                                           aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                                              <  1:3564746/37‑1 (MQ=255)
                                           aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                               >  1:4535093/1‑68 (MQ=255)
                                             gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg                          <  1:4350307/71‑4 (MQ=255)
                                             gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg                          <  1:3895082/71‑4 (MQ=255)
                                              cagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga                         <  1:2118396/71‑1 (MQ=255)
                                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                         >  1:1982035/1‑38 (MQ=255)
                                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                         >  1:3094628/1‑38 (MQ=255)
                                                gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                    <  1:553177/58‑1 (MQ=255)
                                                gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                    <  1:1710426/58‑1 (MQ=255)
                                                 caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                           <  1:154232/66‑1 (MQ=255)
                                                 caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                           <  1:312338/66‑1 (MQ=255)
                                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                         >  1:3000644/1‑35 (MQ=255)
                                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                         >  1:1423125/1‑35 (MQ=255)
                                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                        >  1:1216229/1‑68 (MQ=255)
                                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                        >  1:3710429/1‑68 (MQ=255)
                                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                     >  1:3078264/1‑71 (MQ=255)
                                                   ggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                    >  1:905889/2‑55 (MQ=255)
                                                   agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                     <  1:3330683/70‑1 (MQ=255)
                                                    gcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgt                   <  1:4835971/71‑1 (MQ=255)
                                                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                   <  1:1435471/71‑1 (MQ=255)
                                                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                   >  1:4337394/1‑71 (MQ=255)
                                                     caaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATc                                              >  1:41501/1‑43 (MQ=12)
                                                      aaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgtga                 <  1:4618897/71‑2 (MQ=255)
                                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCa                                                  <  1:2824134/37‑1 (MQ=255)
                                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCa                                             <  1:4323708/42‑1 (MQ=255)
                                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                <  1:5588921/71‑1 (MQ=255)
                                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                <  1:1082422/71‑1 (MQ=255)
                                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                >  1:2302400/1‑71 (MQ=255)
                                                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                    <  1:1385924/50‑1 (MQ=255)
                                                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                           <  1:1803017/59‑1 (MQ=255)
                                                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                           <  1:1820379/59‑1 (MQ=255)
                                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                                                 >  1:3662086/1‑35 (MQ=255)
                                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                            >  1:1150438/1‑40 (MQ=255)
                                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                            >  1:4247255/1‑40 (MQ=255)
                                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                            >  1:5055751/1‑40 (MQ=255)
                                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                       >  1:2581567/1‑45 (MQ=255)
                                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt             <  1:2751295/71‑1 (MQ=255)
                                                          gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa              <  1:940828/70‑1 (MQ=255)
                                                          gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGAt             >  1:2137280/1‑71 (MQ=255)
                                                             tAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGa          <  1:3301044/71‑1 (MQ=255)
                                                              aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                                 <  1:5126632/47‑1 (MQ=255)
                                                              aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                                 <  1:2223316/47‑1 (MQ=255)
                                                               gTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                                   <  1:2293534/44‑1 (MQ=255)
                                                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                          <  1:5290414/36‑1 (MQ=255)
                                                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                          <  1:3235701/36‑1 (MQ=255)
                                                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                      <  1:235236/40‑1 (MQ=255)
                                                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                          <  1:5667610/52‑1 (MQ=255)
                                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                                  >  1:770698/1‑42 (MQ=255)
                                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa                             <  1:4137198/47‑1 (MQ=255)
                                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt     >  1:3179307/1‑71 (MQ=255)
                                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt     >  1:3379225/1‑71 (MQ=255)
                                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt     >  1:705335/1‑71 (MQ=255)
                                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt     >  1:445539/1‑71 (MQ=255)
                                                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                       <  1:2456467/52‑1 (MQ=255)
                                                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                       <  1:5640415/52‑1 (MQ=255)
                                                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                       <  1:1334793/52‑1 (MQ=255)
                                                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt     >  1:3116851/1‑70 (MQ=255)
                                                                     gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa   <  1:845802/70‑1 (MQ=255)
                                                                      |                                                                     
CGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAG  >  minE/608949‑609088

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: