Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,469 G→A 6.0% Y595Y (TAC→TAT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,4690GA 6.0% 248.6 / 3.7 83Y595Y (TAC→TATpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (37/41);  new base A (3/2);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

ACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCG  >  minE/608400‑608539
                                                                     |                                                                      
aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                                                        >  1:115970/1‑70 (MQ=255)
    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                                                        >  1:3634530/1‑66 (MQ=255)
    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                                                        >  1:1062633/1‑66 (MQ=255)
    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                                                   >  1:4710746/1‑71 (MQ=255)
    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                                                   >  1:2879597/1‑71 (MQ=255)
    gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                                                   >  1:2541054/1‑71 (MQ=255)
    gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGgt                                                                   >  1:2102100/1‑71 (MQ=255)
       ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgtgt                                                                 <  1:4831545/70‑1 (MQ=255)
       ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg                                                                >  1:1910693/1‑71 (MQ=255)
            aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCa                                                              <  1:3568046/68‑1 (MQ=255)
            aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt                                                            >  1:1920407/1‑70 (MQ=255)
            aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt                                                            <  1:4449503/70‑1 (MQ=255)
             ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATATCGTCACGGTAGGTGTGCAGttt                                                          <  1:82175/71‑1 (MQ=255)
             ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAg                                                             <  1:5049511/68‑1 (MQ=255)
               aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGt                                                                       <  1:3255532/56‑1 (MQ=255)
               aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCt                                                        >  1:1650878/1‑71 (MQ=255)
                  caccacGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAgg                                                                    >  1:471017/1‑56 (MQ=255)
                  cacaACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGaa                                                     >  1:2005261/1‑71 (MQ=255)
                     aaCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  >  1:932886/1‑71 (MQ=255)
                     aaCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  >  1:1961911/1‑71 (MQ=255)
                     aaCGTTAGAAGTGATGGTCAGAACAGACTCAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  <  1:1590691/71‑1 (MQ=255)
                      aCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  >  1:5092046/1‑70 (MQ=255)
                       cGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGGTTCTGAATTttc                                                <  1:4513366/71‑1 (MQ=255)
                         ttAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  <  1:1665726/67‑1 (MQ=255)
                         ttAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  <  1:478686/67‑1 (MQ=255)
                         ttAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTttctt                                              <  1:3191103/71‑1 (MQ=255)
                           aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa                                            >  1:1923857/1‑71 (MQ=255)
                           aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa                                            >  1:2664794/1‑71 (MQ=255)
                           aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa                                            >  1:4144598/1‑71 (MQ=255)
                           aGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCa                                            >  1:3691980/1‑71 (MQ=255)
                            gAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCAt                                           <  1:3737269/71‑1 (MQ=255)
                            gAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTGAATTTTCTTCAt                                           <  1:903448/71‑1 (MQ=255)
                                tGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTttc                                                <  1:1842762/62‑1 (MQ=255)
                                tGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGa                                         <  1:5430284/69‑1 (MQ=255)
                                tGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGa                                         <  1:2555051/69‑1 (MQ=255)
                                 gATGGTCAGAACAGACTGAGTCGGGATAGAGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                                      >  1:1607965/1‑71 (MQ=255)
                                  aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAg                                                                     >  1:1619658/1‑39 (MQ=255)
                                  aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                                      <  1:5169601/70‑1 (MQ=255)
                                  aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                                      <  1:2557183/70‑1 (MQ=255)
                                  aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg                                     >  1:4177682/1‑71 (MQ=255)
                                  aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg                                     <  1:411506/71‑1 (MQ=255)
                                     gTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGGGTGCAGTTTCTGAATTTTCTTCATGAAACGtt                                   <  1:470384/70‑1 (MQ=255)
                                      tCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCt                                 <  1:836843/71‑1 (MQ=255)
                                       cAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  >  1:512016/1‑53 (MQ=255)
                                       cAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  >  1:3275993/1‑53 (MQ=255)
                                        aGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTa                                >  1:457559/1‑70 (MQ=255)
                                         gAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTAcc                              <  1:4320318/71‑1 (MQ=255)
                                          aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                                      <  1:727378/62‑1 (MQ=255)
                                          aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                                      <  1:4152402/62‑1 (MQ=255)
                                          aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                                      <  1:3001829/62‑1 (MQ=255)
                                           aCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAg                                                             >  1:5564359/1‑38 (MQ=255)
                                           aCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTg                                                       >  1:2980204/1‑44 (MQ=255)
                                            cAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAgg                           >  1:2523699/1‑71 (MQ=255)
                                               aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                                      <  1:1551903/57‑1 (MQ=255)
                                               aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCa                        <  1:712001/71‑1 (MQ=255)
                                               aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCa                        <  1:5237242/71‑1 (MQ=255)
                                               aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGATCTACCAGGTc                         >  1:2196096/1‑70 (MQ=255)
                                                cTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  <  1:4267765/44‑1 (MQ=255)
                                                cTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  <  1:4296603/44‑1 (MQ=255)
                                                 tGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  <  1:1389460/43‑1 (MQ=255)
                                                     tCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTc                         <  1:2804514/64‑1 (MQ=255)
                                                      cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                                  <  1:5526022/38‑1 (MQ=255)
                                                      cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGt                                    <  1:3884306/52‑1 (MQ=255)
                                                      cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGt                                    <  1:1309733/52‑1 (MQ=255)
                                                       gggATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTGTACCAGGTCAACAGCCAg                >  1:3342243/1‑71 (MQ=255)
                                                       gggATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg                <  1:4270063/71‑1 (MQ=255)
                                                       gggATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg                >  1:4812152/1‑71 (MQ=255)
                                                       gggATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg                >  1:3315199/1‑71 (MQ=255)
                                                       gggATAGCGTCAAGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg                >  1:1855319/1‑71 (MQ=255)
                                                          aTAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTc             <  1:5285724/71‑1 (MQ=255)
                                                            aGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg                                     >  1:4113556/1‑45 (MQ=255)
                                                             gCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCaa                       >  1:4700207/1‑58 (MQ=255)
                                                             gCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCaa                       >  1:2090656/1‑58 (MQ=255)
                                                             gCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCaa                       >  1:3637429/1‑58 (MQ=255)
                                                             gCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATc          <  1:3679454/71‑1 (MQ=255)
                                                              cGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTa                                <  1:2057071/48‑1 (MQ=255)
                                                              cGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCAt           <  1:418698/69‑1 (MQ=255)
                                                              agTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTa                                <  1:965159/47‑1 (MQ=255)
                                                               gTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTa        >  1:3990721/1‑71 (MQ=255)
                                                               gTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTa        <  1:2031069/71‑1 (MQ=255)
                                                                tCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTa        >  1:1616694/1‑70 (MQ=255)
                                                                 cACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTaca      >  1:2079075/1‑71 (MQ=255)
                                                                  aCGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAg                            >  1:4443120/1‑48 (MQ=255)
                                                                  aCGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTacac     >  1:4155460/1‑71 (MQ=255)
                                                                  aCGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTacac     >  1:4110448/1‑71 (MQ=255)
                                                                   cGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCaa                       <  1:3484473/52‑1 (MQ=255)
                                                                     gTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACAcgcg  >  1:3661614/1‑71 (MQ=255)
                                                                     |                                                                      
ACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCG  >  minE/608400‑608539

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: