Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,869,681 A→G 23.9% G25G (GGT→GGC ygcW ← predicted deoxygluconate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,869,6810AG23.9% 156.3 / 44.0 125G25G (GGT→GGCygcWpredicted deoxygluconate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (30/65);  new base G (8/22);  total (38/87)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 125 total aligned reads displayed.

AAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTC  >  minE/1869619‑1869746
                                                              |                                                                 
aaGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGGCTAAACCGCTATTCCCACCGGt                                                              >  1:116834/1‑68 (MQ=255)
aaGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCgg                                                               >  1:1825929/1‑67 (MQ=255)
aaGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGt                                                              >  1:96967/1‑68 (MQ=255)
aaGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGt                                                              >  1:2922858/1‑68 (MQ=255)
aaGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGt                                                              >  1:2782066/1‑68 (MQ=255)
aaGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGt                                                              >  1:393880/1‑68 (MQ=255)
aaGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACAGCTATTCCCACCGGt                                                              >  1:638579/1‑68 (MQ=255)
     atTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAAt                                                        >  1:1522154/1‑69 (MQ=255)
              aGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                                                <  1:276526/68‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:3108313/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:3178940/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:3034978/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:174440/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:2881763/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:2733804/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:2729739/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:2517570/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:2410810/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:2258661/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:1195339/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:8236/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:3414354/69‑1 (MQ=255)
                cTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:354276/69‑1 (MQ=255)
                     gCCAACCCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                                             <  1:641623/64‑1 (MQ=255)
                       cAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGTTTTACCTTTCa                                       >  1:774826/1‑68 (MQ=255)
                       cAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGTTTTACCTTTCa                                       >  1:1013997/1‑68 (MQ=255)
                       cAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGTTTTACCTTTCa                                       >  1:968747/1‑68 (MQ=255)
                       cAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGTTTTACCTTTCa                                       >  1:3172870/1‑68 (MQ=255)
                       cAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGTTTTACCTTTCa                                       >  1:1285854/1‑68 (MQ=255)
                       cAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGTTTTACCTTTCa                                       >  1:3251363/1‑68 (MQ=255)
                       cAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGTTTTACCTTTCa                                       >  1:1718508/1‑68 (MQ=255)
                                ggCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGaaaaa                             >  1:3276240/1‑69 (MQ=255)
                                ggCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGaaaaa                             >  1:2544904/1‑69 (MQ=255)
                                 gCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGaaaaaa                            <  1:1224590/69‑1 (MQ=255)
                                 gCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGaaaaaa                            <  1:2546252/69‑1 (MQ=255)
                                   aaaTGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATc                          <  1:2681938/69‑1 (MQ=255)
                                   aaaTGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATc                          <  1:1095533/69‑1 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                                                <  1:1153671/44‑1 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                                                <  1:2176993/44‑1 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                                                <  1:3062160/44‑1 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                                                <  1:811719/44‑1 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                                                <  1:1724039/44‑1 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                                                <  1:636841/44‑1 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                                                <  1:3280522/44‑1 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                                                <  1:3380067/44‑1 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCa                        >  1:2198818/1‑68 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCa                        >  1:281005/1‑68 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCa                        >  1:2525530/1‑68 (MQ=255)
                                      tGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCa                        >  1:429511/1‑68 (MQ=255)
                                       gCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCa                        >  1:1827387/1‑67 (MQ=255)
                                          tGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTgag                   <  1:595607/69‑1 (MQ=255)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:1465006/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:552750/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:1711862/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:696060/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:892993/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:2191763/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:2955121/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:2949505/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:300934/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:2785283/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:2024994/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:1891525/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:2439262/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:269612/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:1947173/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:2460976/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:2232915/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:2474533/38‑1 (MQ=38)
                                            gCCTAAACCGCTATTCCCGCCGGTAACAATTGCTGttt                                                <  1:1952804/38‑1 (MQ=38)
                                              cTAAACCGCTATTCCCACCTGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACg               <  1:180449/69‑1 (MQ=255)
                                              cTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACg               <  1:817451/69‑1 (MQ=255)
                                              cTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACg               <  1:1348384/69‑1 (MQ=255)
                                              cTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACg               <  1:545412/69‑1 (MQ=255)
                                              cTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACg               <  1:1982120/69‑1 (MQ=255)
                                              cTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACg               <  1:198034/69‑1 (MQ=255)
                                              cTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACg               <  1:2809803/69‑1 (MQ=255)
                                               tAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACg               <  1:1796456/68‑1 (MQ=255)
                                                   ccGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCAtt           <  1:1593560/68‑1 (MQ=255)
                                                   ccGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCAtt           <  1:159672/68‑1 (MQ=255)
                                                   ccGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCAtt           <  1:2402209/68‑1 (MQ=255)
                                                   ccGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCAtt           <  1:1449806/68‑1 (MQ=255)
                                                   ccGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCAtt           <  1:1862522/68‑1 (MQ=255)
                                                   ccGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCAtt           <  1:130939/68‑1 (MQ=255)
                                                      cTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTgaga       <  1:2551814/69‑1 (MQ=255)
                                                      cTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTgaga       <  1:2326139/69‑1 (MQ=255)
                                                      cTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTgaga       <  1:1274755/69‑1 (MQ=255)
                                                      cTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTgaga       <  1:2592337/69‑1 (MQ=255)
                                                      cTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTgaga       <  1:3084592/69‑1 (MQ=255)
                                                      cTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTgaga       <  1:1251514/69‑1 (MQ=255)
                                                      cTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTgaga       <  1:172693/69‑1 (MQ=255)
                                                           cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTgaga       >  1:1563391/1‑64 (MQ=255)
                                                           cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGAttt  >  1:2505553/1‑68 (MQ=255)
                                                           cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTc  >  1:1208187/1‑69 (MQ=255)
                                                           cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTc  >  1:977015/1‑69 (MQ=255)
                                                           cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTc  >  1:579275/1‑69 (MQ=255)
                                                           cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTc  >  1:1507272/1‑69 (MQ=255)
                                                           cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTc  >  1:3149574/1‑69 (MQ=255)
                                                           cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTc  >  1:2604058/1‑69 (MQ=255)
                                                           cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTc  >  1:1918813/1‑69 (MQ=255)
                                                              |                                                                 
AAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTC  >  minE/1869619‑1869746

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: