Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,978 C→T 18.4% V426I (GTT→ATT)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,9780CT18.4% 335.7 / 60.4 152V426I (GTT→ATT) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (46/78);  new base T (15/13);  total (61/91)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 155 total aligned reads displayed.

CGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCG  >  minE/608909‑609048
                                                                     |                                                                      
cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                                                       <  1:775996/71‑1 (MQ=255)
cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                                                       <  1:2771617/71‑1 (MQ=255)
  ttGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                                                     >  1:4718237/1‑71 (MQ=255)
     aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                                   <  1:3114830/70‑1 (MQ=255)
     aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                                   <  1:261160/70‑1 (MQ=255)
        gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                               >  1:3865652/1‑71 (MQ=255)
        gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                               >  1:5340523/1‑71 (MQ=255)
         cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                               <  1:493001/70‑1 (MQ=255)
             cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                          >  1:5663115/1‑71 (MQ=255)
              aGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAc                                                         <  1:196137/71‑1 (MQ=255)
                cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                          <  1:4089256/68‑1 (MQ=255)
                 aTGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                          >  1:3132667/1‑67 (MQ=255)
                  tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag                                                     >  1:814231/1‑71 (MQ=255)
                  tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACAATCATCGGGCTTACg                                                        <  1:508464/68‑1 (MQ=255)
                   cgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                                     >  1:1143248/2‑70 (MQ=255)
                    gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                                                    <  1:1461504/70‑1 (MQ=255)
                    gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                                                    >  1:1569553/1‑70 (MQ=255)
                     ttttCGCCAGGTTTGCACGCGCACCGATGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                                          <  1:4923260/63‑1 (MQ=255)
                     ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                                   <  1:576990/54‑1 (MQ=255)
                     ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                                   <  1:399597/54‑1 (MQ=255)
                     ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa                                                  <  1:4812571/71‑1 (MQ=255)
                      tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                                                 >  1:5117703/1‑71 (MQ=255)
                      tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                                                 >  1:5591309/1‑71 (MQ=255)
                       ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                                                 >  1:2391923/1‑70 (MQ=255)
                       ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                                 >  1:5207214/1‑68 (MQ=255)
                       ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                                 >  1:5184458/1‑68 (MQ=255)
                       ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                                 >  1:5332640/1‑68 (MQ=255)
                       ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                                                 >  1:32789/1‑68 (MQ=255)
                        tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                <  1:2032018/70‑1 (MQ=255)
                        tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                <  1:202592/70‑1 (MQ=255)
                        tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                <  1:5547664/70‑1 (MQ=255)
                        tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                <  1:3043363/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)
                              ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                                    >  1:1341037/1‑44 (MQ=255)
                              ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                                    <  1:4240566/44‑1 (MQ=255)
                              ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                >  1:4642752/1‑64 (MQ=255)
                              ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                         <  1:224252/71‑1 (MQ=255)
                               gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                            <  1:4766030/51‑1 (MQ=255)
                               gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                            <  1:4772410/51‑1 (MQ=255)
                                tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                                                        <  1:1089070/38‑1 (MQ=255)
                                   gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                                    <  1:2935523/39‑1 (MQ=255)
                                   gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                                    <  1:5612882/39‑1 (MQ=255)
                                   gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                               <  1:2511788/44‑1 (MQ=255)
                                   gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                               <  1:4919366/44‑1 (MQ=255)
                                     aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                                  >  1:4914727/1‑39 (MQ=255)
                                     aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                                  >  1:827949/1‑39 (MQ=255)
                                     aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                                   <  1:252889/38‑1 (MQ=255)
                                     aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                                   <  1:5496042/38‑1 (MQ=255)
                                       gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg                                <  1:3806412/71‑1 (MQ=255)
                                       gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg                                <  1:1110034/71‑1 (MQ=255)
                                         gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                              >  1:1677563/1‑71 (MQ=255)
                                         gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                              >  1:3664825/1‑71 (MQ=255)
                                         gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                              >  1:4528867/1‑71 (MQ=255)
                                         gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAgcaggcaa                                              <  1:4796127/55‑1 (MQ=255)
                                          cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                          <  1:3759162/58‑1 (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)
                                             cGAAGATCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                          <  1:2340665/71‑1 (MQ=255)
                                             cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgca                                                      <  1:3427250/43‑1 (MQ=255)
                                             cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                         <  1:192719/56‑1 (MQ=255)
                                             cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                          <  1:4459885/71‑1 (MQ=255)
                                             cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                          <  1:3985486/71‑1 (MQ=255)
                                             cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                          <  1:2209975/71‑1 (MQ=255)
                                              gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                                   >  1:4493100/1‑45 (MQ=255)
                                              gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                                   >  1:2836770/1‑45 (MQ=255)
                                              gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                                   >  1:2612485/1‑45 (MQ=255)
                                              gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGa                         <  1:217997/71‑1 (MQ=255)
                                               aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                            >  1:1469284/1‑51 (MQ=255)
                                               aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                          <  1:3781532/53‑1 (MQ=255)
                                               aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                          >  1:4654861/1‑69 (MQ=255)
                                               aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc                                            >  1:4440321/1‑50 (MQ=255)
                                               aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                          >  1:2232396/1‑69 (MQ=255)
                                                agaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                                                <  1:3829407/46‑1 (MQ=16)
                                                 gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                      >  1:24494/1‑71 (MQ=255)
                                                  aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                <  1:5574477/44‑1 (MQ=255)
                                                  aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgc                               <  1:3580464/61‑2 (MQ=255)
                                                  aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                     <  1:1569837/71‑1 (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)
                                                      gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAAt                                             >  1:515316/1‑43 (MQ=255)
                                                       cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                              >  1:4831019/1‑41 (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)
                                                            gTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                 >  1:5244188/1‑65 (MQ=255)
                                                             tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                             <  1:172749/52‑1 (MQ=255)
                                                             tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt            >  1:4548540/1‑69 (MQ=255)
                                                               tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg         <  1:1626106/70‑1 (MQ=255)
                                                                aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg        <  1:3253912/70‑1 (MQ=255)
                                                                 ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt       <  1:1329495/70‑1 (MQ=255)
                                                                   aaCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt    <  1:2096741/71‑1 (MQ=255)
                                                                     cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg  <  1:4051844/71‑1 (MQ=255)
                                                                     |                                                                      
CGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCG  >  minE/608909‑609048

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: