Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,126,307 A→C 71.5% V109V (GTT→GTG pheS ← phenylalanine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,126,3070AC71.5% 55.9 / 20.5 42V109V (GTT→GTGpheSphenylalanine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (3/9);  new base C (20/10);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCC  >  minE/1126243‑1126370
                                                                |                                                               
gCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                                             <  1:3252901/69‑1 (MQ=255)
gCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCg                                                             <  1:1888421/69‑1 (MQ=255)
      ttGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCa                                                       <  1:2857263/69‑1 (MQ=255)
      ttGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCa                                                       <  1:2338267/69‑1 (MQ=255)
      ttGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCa                                                       <  1:22968/69‑1 (MQ=255)
              gTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGAccgccg                                               <  1:321473/69‑1 (MQ=255)
               tAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGAccgccg                                               <  1:1397440/68‑1 (MQ=255)
                   gCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTc                                          <  1:1860737/69‑1 (MQ=255)
                    ccAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCa                                         >  1:2132630/1‑69 (MQ=255)
                    ccAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCa                                         >  1:3205970/1‑69 (MQ=255)
                         cTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGc                                     >  1:2066082/1‑68 (MQ=255)
                              cGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGcgacga                                <  1:2864777/68‑1 (MQ=255)
                                                 gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:171666/68‑1 (MQ=255)
                                                 gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:75833/68‑1 (MQ=255)
                                                 gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:1332349/68‑1 (MQ=255)
                                                 gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:1379135/68‑1 (MQ=255)
                                                 gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:2194106/68‑1 (MQ=255)
                                                 gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:561889/68‑1 (MQ=255)
                                                 gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:1544558/68‑1 (MQ=255)
                                                 gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:879628/68‑1 (MQ=255)
                                                  tCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:882632/67‑1 (MQ=255)
                                                  tCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa             <  1:1073963/67‑1 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGCCGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:2966534/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATc               >  1:731398/1‑55 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGttt       >  1:2339030/1‑63 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:526925/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:3403063/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:3329/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:812173/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:2498043/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:740394/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:2726898/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:2604776/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:2599318/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:2267239/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:1998917/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:1926931/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:1840568/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:1643913/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:1489509/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:1262032/1‑68 (MQ=255)
                                                            gggTCACCGGATGCAGCCAGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc  >  1:3399510/1‑68 (MQ=255)
                                                                |                                                               
GCCCGGTTGCCACGGTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCC  >  minE/1126243‑1126370

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: