Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,970 T→C 23.3% K428K (AAA→AAG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9700TC23.3% 292.6 / 87.4 159K428K (AAA→AAG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (47/75);  new base C (21/16);  total (68/91)
Fisher's exact test for biased strand distribution p-value = 5.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

AACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAT  >  minE/608901‑609038
                                                                     |                                                                    
aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                                     <  1:898099/71‑1 (MQ=255)
aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                                     >  1:2669789/1‑71 (MQ=255)
 aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                                     >  1:3459613/1‑70 (MQ=255)
   cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                                                  >  1:3756992/2‑71 (MQ=255)
     cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                                                                <  1:2065457/71‑1 (MQ=255)
     cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                                                 >  1:5186119/1‑70 (MQ=255)
       ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                                                   >  1:5516307/1‑66 (MQ=255)
       ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaat                                                              <  1:2375314/71‑2 (MQ=255)
        cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                                             <  1:2771617/71‑1 (MQ=255)
        cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                                             <  1:775996/71‑1 (MQ=255)
          ttGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                                           >  1:4718237/1‑71 (MQ=255)
            gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:2818379/70‑1 (MQ=255)
            gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:4810051/70‑1 (MQ=255)
             aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         <  1:261160/70‑1 (MQ=255)
             aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         <  1:3114830/70‑1 (MQ=255)
               tGCGTACAGCATGGTTTTCCCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTa                                                                   >  1:3366110/1‑58 (MQ=255)
                   tACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                  <  1:4018357/71‑1 (MQ=255)
                     cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                >  1:2240311/1‑71 (MQ=255)
                     cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                >  1:5572364/1‑71 (MQ=255)
                      aGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAc                                               <  1:196137/71‑1 (MQ=255)
                         aTGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                >  1:3132667/1‑67 (MQ=255)
                          tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag                                           >  1:814231/1‑71 (MQ=255)
                           cgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                           >  1:1143248/2‑70 (MQ=255)
                            gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                        <  1:2207997/56‑1 (MQ=255)
                            gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                                          >  1:1569553/1‑70 (MQ=255)
                            gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                                          <  1:1461504/70‑1 (MQ=255)
                             ttttCGCCAGGTTTGCACGCGCACCGATGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                <  1:4923260/63‑1 (MQ=255)
                             ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         <  1:399597/54‑1 (MQ=255)
                             ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         <  1:576990/54‑1 (MQ=255)
                             ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCtt                                                 <  1:1943784/62‑1 (MQ=255)
                             ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa                                        <  1:5109958/71‑1 (MQ=255)
                              tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          >  1:5134384/1‑52 (MQ=255)
                              tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                                       >  1:5591309/1‑71 (MQ=255)
                               ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                                       >  1:2391923/1‑70 (MQ=255)
                               ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                       >  1:5184458/1‑68 (MQ=255)
                               ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                       >  1:32789/1‑68 (MQ=255)
                               ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                       >  1:5207214/1‑68 (MQ=255)
                               ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                       >  1:5332640/1‑68 (MQ=255)
                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                  <  1:5587630/58‑1 (MQ=255)
                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:2032018/70‑1 (MQ=255)
                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:3043363/70‑1 (MQ=255)
                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:5274704/70‑1 (MQ=255)
                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:5547664/70‑1 (MQ=255)
                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:202592/70‑1 (MQ=255)
                                tCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                                      <  1:4334380/70‑1 (MQ=255)
                                 cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                     <  1:2937922/70‑1 (MQ=255)
                                  gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    <  1:4224293/70‑1 (MQ=255)
                                  gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    <  1:4653418/70‑1 (MQ=255)
                                    cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaatc                                                             <  1:1905524/43‑3 (MQ=19)
                                     aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                >  1:4947648/1‑71 (MQ=255)
                                      ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:4240566/44‑1 (MQ=255)
                                      ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          >  1:1341037/1‑44 (MQ=255)
                                       gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                  <  1:4772410/51‑1 (MQ=255)
                                        tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                                              <  1:1089070/38‑1 (MQ=255)
                                        tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                        <  1:1834487/44‑1 (MQ=255)
                                           gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:2935523/39‑1 (MQ=255)
                                           gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                     <  1:2511788/44‑1 (MQ=255)
                                           gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                     <  1:4919366/44‑1 (MQ=255)
                                            cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:2138586/58‑1 (MQ=255)
                                             aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                        >  1:4914727/1‑39 (MQ=255)
                                             aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                        >  1:1657053/1‑39 (MQ=255)
                                             aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         <  1:252889/38‑1 (MQ=255)
                                               gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg                      <  1:1110034/71‑1 (MQ=255)
                                                 gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                    >  1:3664825/1‑71 (MQ=255)
                                                 gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAgcaggcaa                                    <  1:4796127/55‑1 (MQ=255)
                                                  cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                   <  1:2154030/71‑1 (MQ=255)
                                                  cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctg                   >  1:3185884/1‑71 (MQ=255)
                                                   aCCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                          <  1:5172893/47‑1 (MQ=16)
                                                     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)
                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                               <  1:192719/56‑1 (MQ=255)
                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                <  1:4459885/71‑1 (MQ=255)
                                                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                         >  1:4493100/1‑45 (MQ=255)
                                                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                         >  1:2612485/1‑45 (MQ=255)
                                                      aaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                >  1:4544652/2‑70 (MQ=255)
                                                      aaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                >  1:633249/2‑70 (MQ=255)
                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                  >  1:1469284/1‑51 (MQ=255)
                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                  >  1:4730602/1‑51 (MQ=255)
                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                  >  1:4163056/1‑51 (MQ=255)
                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                <  1:3781532/53‑1 (MQ=255)
                                                       aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                          >  1:1492323/1‑43 (MQ=16)
                                                       aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc                                  >  1:4440321/1‑50 (MQ=255)
                                                       aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                >  1:2232396/1‑69 (MQ=255)
                                                        agaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                                      <  1:3828394/46‑1 (MQ=16)
                                                        agaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                                      <  1:3829407/46‑1 (MQ=16)
                                                         gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt            >  1:24494/1‑71 (MQ=255)
                                                          aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgc                     <  1:3580464/61‑2 (MQ=255)
                                                          aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc           <  1:1569837/71‑1 (MQ=255)
                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    >  1:5648645/1‑44 (MQ=255)
                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                <  1:2567882/64‑1 (MQ=255)
                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa              >  1:5081927/1‑66 (MQ=255)
                                                             tGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTc                          <  1:762743/53‑1 (MQ=255)
                                                             tGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTc                          <  1:2918683/53‑1 (MQ=255)
                                                              gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAAt                                   >  1:1195933/1‑43 (MQ=255)
                                                               cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    >  1:4831019/1‑41 (MQ=255)
                                                                  ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTc                          <  1:3881398/48‑1 (MQ=255)
                                                                     tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                   <  1:172749/52‑1 (MQ=255)
                                                                     tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt  >  1:4548540/1‑69 (MQ=255)
                                                                     |                                                                    
AACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAT  >  minE/608901‑609038

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