Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,714,265 A→G 53.0% G389G (GGT→GGC yphG ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,714,2650AG53.0% 8.5 / 70.1 83G389G (GGT→GGCyphGconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base G (38/6);  total (55/28)
Fisher's exact test for biased strand distribution p-value = 5.53e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

GATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTC  >  minE/1714211‑1714327
                                                      |                                                              
gATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGAt                                                  >  1:3083748/1‑69 (MQ=255)
gATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGAt                                                  >  1:2635151/1‑69 (MQ=255)
gATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGAt                                                  >  1:2688924/1‑69 (MQ=255)
gATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGAt                                                  >  1:698738/1‑69 (MQ=255)
gATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGAt                                                  >  1:1761525/1‑69 (MQ=255)
gATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGAt                                                  >  1:1388155/1‑69 (MQ=255)
  tGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATc                                                 >  1:2094588/1‑68 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCg                                                             >  1:1813407/1‑51 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:2608216/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:282259/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:2631951/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:167189/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:2950377/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:925380/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:194584/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:3167418/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:2748582/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:325336/1‑69 (MQ=255)
       gTCGCTGGCATCAGTGCCAAGGCATCATCAAGTAACGCGTAGCATTTGCCGATTTCGCGGATCGCCAgg                                           >  1:291503/1‑69 (MQ=255)
                aTCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCg                                  >  1:871734/1‑69 (MQ=255)
                   aGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTc                               >  1:2495027/1‑69 (MQ=255)
                   aGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTc                               >  1:2330607/1‑69 (MQ=255)
                   aGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTc                               >  1:2195038/1‑69 (MQ=255)
                   aGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTc                               >  1:75102/1‑69 (MQ=255)
                         aCGGCATGATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:144312/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTc                                                        >  1:2710799/1‑38 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCgg                                                    >  1:2173363/1‑42 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGt                                       >  1:1111006/1‑55 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2324686/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2753187/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:3429649/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2777530/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:311187/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2861030/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2542195/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2534468/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2886389/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:555129/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:1121546/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:1123664/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:1209598/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:769877/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2171583/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:209427/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:1432519/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:1543212/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:1789886/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAACGCGTTGCATTTGACGATTTCGCGGATc                                                 >  1:1325133/1‑45 (MQ=255)
                         aCGGCATCATCAAGTAACGCGCTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2138640/1‑64 (MQ=255)
                         aCGGCATCATCAAGTAAAGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCa                              >  1:2342894/1‑64 (MQ=255)
                            gCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgaga                       <  1:1076197/68‑1 (MQ=255)
                            gCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgaga                       <  1:2988342/68‑1 (MQ=255)
                            gCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgaga                       <  1:1992002/68‑1 (MQ=255)
                            gCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgaga                       <  1:2070819/68‑1 (MQ=255)
                            gCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgaga                       <  1:2347792/68‑1 (MQ=255)
                            gCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgaga                       <  1:3434902/68‑1 (MQ=255)
                            gCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgaga                       <  1:395512/68‑1 (MQ=255)
                             catcatCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgagaga                     <  1:362030/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:556697/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:2600821/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:487635/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:2367651/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGGTCAACGGAGAGATg                   <  1:778048/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:977492/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:212674/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:637760/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:2309152/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:2542941/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:2723261/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:3310334/69‑1 (MQ=255)
                               tcatcaAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATg                   <  1:2744849/69‑1 (MQ=255)
                                catcaAGTAACGCGTTGCATTTGCCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGAtt                   <  1:1477923/68‑2 (MQ=255)
                                   caAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGtt                                >  1:805895/1‑52 (MQ=255)
                                   caAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGct                <  1:2581680/68‑2 (MQ=255)
                                   caAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCa                >  1:2751586/1‑68 (MQ=255)
                                   caAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCa                >  1:1161102/1‑68 (MQ=255)
                                    aaGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATa              <  1:1055164/69‑1 (MQ=255)
                                    aaGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATa              <  1:455284/69‑1 (MQ=255)
                                    aaGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATa              <  1:2422122/69‑1 (MQ=255)
                                    aaGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATa              <  1:1317778/69‑1 (MQ=255)
                                    aaGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATa              <  1:3032398/69‑1 (MQ=255)
                                                gCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACtg  >  1:3040398/1‑68 (MQ=255)
                                                gCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTc  >  1:1195094/1‑69 (MQ=255)
                                                      |                                                              
GATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTC  >  minE/1714211‑1714327

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: