Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,970 T→C 10.9% 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,9700TC10.9% 302.2 / 24.7 110K428K (AAA→AAG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (54/44);  new base C (5/7);  total (59/51)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 114 total aligned reads displayed.

AACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAT  >  minE/608901‑609038
                                                                     |                                                                    
aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                                     >  1:1856739/1‑71 (MQ=255)
 aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt                                                                    <  1:3815239/71‑1 (MQ=255)
  ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                                                   <  1:2023686/69‑1 (MQ=255)
  ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                                                   <  1:725150/69‑1 (MQ=255)
      gccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCaa                                                               >  1:128065/1‑71 (MQ=255)
        cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                                              >  1:953683/1‑70 (MQ=255)
        cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                                              >  1:3075564/1‑70 (MQ=255)
        cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                                             <  1:677044/71‑1 (MQ=255)
           tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                                           >  1:657558/1‑70 (MQ=255)
            gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         >  1:587614/1‑71 (MQ=255)
            gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                         >  1:1852293/1‑71 (MQ=255)
            gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:4005748/70‑1 (MQ=255)
              ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                                       <  1:910646/71‑1 (MQ=255)
              ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                                       <  1:1823578/71‑1 (MQ=255)
                    aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                        >  1:111947/1‑64 (MQ=255)
                    aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAAGatcat                                                        >  1:1450022/1‑64 (MQ=255)
                       gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                  >  1:1288974/1‑67 (MQ=255)
                       gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                                              <  1:43335/71‑1 (MQ=255)
                       gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                                              <  1:2648335/71‑1 (MQ=255)
                          tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                                                                <  1:1360880/50‑1 (MQ=255)
                          tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                     <  1:4123172/61‑1 (MQ=255)
                          tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                     <  1:2636176/61‑1 (MQ=255)
                          tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                     <  1:1040688/61‑1 (MQ=255)
                          tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCGTACg                                              >  1:3760010/1‑68 (MQ=255)
                            gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                         <  1:3498025/71‑1 (MQ=255)
                             ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgca                                            <  1:2918905/67‑1 (MQ=255)
                             ttttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa                                                               <  1:3893971/48‑1 (MQ=25)
                               ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          >  1:3624840/1‑51 (MQ=255)
                               ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                          <  1:1990208/51‑1 (MQ=255)
                                tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:3395187/70‑1 (MQ=255)
                                 cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                     <  1:1670833/70‑1 (MQ=255)
                                 cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                     <  1:11636/70‑1 (MQ=255)
                                 cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                     <  1:1199145/70‑1 (MQ=255)
                                 cGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggca                                     <  1:1900468/70‑1 (MQ=255)
                                    cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                 <  1:1472220/71‑1 (MQ=255)
                                    cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGctt                                                                    <  1:2241829/36‑1 (MQ=16)
                                     aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                >  1:4119712/1‑71 (MQ=255)
                                      ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                 >  1:2729084/1‑69 (MQ=255)
                                      ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                 >  1:4078426/1‑69 (MQ=255)
                                      ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                               <  1:2076321/71‑1 (MQ=255)
                                       gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGt                              <  1:1555201/71‑1 (MQ=255)
                                        tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTa                             <  1:232203/71‑1 (MQ=255)
                                         ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    >  1:1117979/1‑63 (MQ=255)
                                         ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                  >  1:2022989/1‑65 (MQ=255)
                                         ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                            >  1:593652/1‑71 (MQ=255)
                                         ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                            >  1:40416/1‑71 (MQ=255)
                                          tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                            >  1:488209/1‑70 (MQ=255)
                                          tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                           >  1:3324975/1‑71 (MQ=255)
                                          tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                           >  1:509491/1‑71 (MQ=255)
                                          tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                           >  1:2525102/1‑71 (MQ=255)
                                          tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                           >  1:2004154/1‑71 (MQ=255)
                                             aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:3363694/57‑1 (MQ=255)
                                               gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    >  1:282008/1‑57 (MQ=255)
                                               gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    >  1:2363300/1‑57 (MQ=255)
                                               gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg                      <  1:2744558/71‑1 (MQ=255)
                                                tgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTa                             <  1:456201/62‑1 (MQ=255)
                                                tgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATAGCGTa                             <  1:802870/62‑1 (MQ=255)
                                                 gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                    >  1:1211478/1‑71 (MQ=255)
                                                 gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                    >  1:3364338/1‑71 (MQ=255)
                                                 gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct                    >  1:808507/1‑69 (MQ=255)
                                                 gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct                    >  1:1385780/1‑69 (MQ=255)
                                                    ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:2944900/50‑1 (MQ=255)
                                                    ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                  >  1:123820/1‑70 (MQ=255)
                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAAt                                   <  1:1006285/52‑1 (MQ=255)
                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                               <  1:2124285/56‑1 (MQ=255)
                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                 >  1:2056076/1‑70 (MQ=255)
                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                <  1:2539150/71‑1 (MQ=255)
                                                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                <  1:151874/71‑1 (MQ=255)
                                                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                 >  1:1840749/1‑69 (MQ=255)
                                                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                >  1:3036046/1‑70 (MQ=255)
                                                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGAAGCAAGCAATAGCGTAGTCATCgttg                   >  1:2618494/1‑67 (MQ=255)
                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                <  1:618985/69‑1 (MQ=255)
                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                >  1:3905432/1‑69 (MQ=255)
                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                >  1:2744738/1‑69 (MQ=255)
                                                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa              <  1:2852718/71‑1 (MQ=255)
                                                       aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGaa              <  1:1272267/71‑1 (MQ=255)
                                                        agaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAg             >  1:3157389/1‑71 (MQ=255)
                                                         gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt            >  1:1475473/1‑71 (MQ=255)
                                                           aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTcc          <  1:1934522/71‑1 (MQ=255)
                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                      <  1:4002379/42‑1 (MQ=255)
                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    >  1:1850420/1‑44 (MQ=255)
                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg         <  1:3254358/71‑1 (MQ=255)
                                                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg         <  1:421580/71‑1 (MQ=255)
                                                              gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                  >  1:1060756/1‑60 (MQ=255)
                                                              gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                  >  1:765660/1‑60 (MQ=255)
                                                               cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    <  1:1817824/41‑1 (MQ=255)
                                                               cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    <  1:244663/41‑1 (MQ=255)
                                                               cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    <  1:3988942/41‑1 (MQ=255)
                                                               cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                >  1:838121/1‑45 (MQ=255)
                                                               cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGAc      >  1:2272008/1‑71 (MQ=255)
                                                                 tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAg             <  1:643970/62‑1 (MQ=255)
                                                                 tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGc    <  1:3672883/71‑1 (MQ=255)
                                                                   tGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                    >  1:2006482/1‑37 (MQ=255)
                                                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCTTAGTCATCGTTGTTGaa              >  1:1225262/1‑59 (MQ=255)
                                                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa              >  1:1600160/1‑59 (MQ=255)
                                                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTc           >  1:56076/1‑62 (MQ=255)
                                                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt  >  1:259203/1‑71 (MQ=255)
                                                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt  >  1:2190831/1‑71 (MQ=255)
                                                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt  >  1:2830485/1‑71 (MQ=255)
                                                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt  >  1:3799987/1‑71 (MQ=255)
                                                                     |                                                                    
AACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAT  >  minE/608901‑609038

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