Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,126,349 A→G 77.6% S95S (TCT→TCC pheS ← phenylalanine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,126,3490AG77.6% 50.1 / 9.1 40S95S (TCT→TCCpheSphenylalanine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (2/7);  new base G (20/11);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 5.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCGCTTTCCAGTTCCGCTTTACGCG  >  minE/1126292‑1126416
                                                         |                                                                   
gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:1332349/68‑1 (MQ=255)
gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:1379135/68‑1 (MQ=255)
gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:879628/68‑1 (MQ=255)
gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:1544558/68‑1 (MQ=255)
gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:171666/68‑1 (MQ=255)
gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:75833/68‑1 (MQ=255)
gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:2194106/68‑1 (MQ=255)
gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:561889/68‑1 (MQ=255)
 tCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:1073963/67‑1 (MQ=255)
 tCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:882632/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:2726898/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:1262032/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:3329/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:3403063/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:526925/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:740394/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:812173/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:2604776/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:1489509/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:1643913/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:1840568/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:1926931/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:1998917/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:2267239/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:2498043/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:2599318/1‑68 (MQ=255)
           gggTCACCGGATGCAGCCAGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc                                                >  1:3399510/1‑68 (MQ=255)
              tcaCCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa                                                           <  1:2658158/52‑1 (MQ=255)
                aCCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGa                                          <  1:3216484/69‑1 (MQ=255)
                aCCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGa                                          <  1:3052475/69‑1 (MQ=255)
                aCCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGa                                          <  1:1344406/69‑1 (MQ=255)
                                    ttCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCGCt                      >  1:707778/1‑69 (MQ=255)
                                    ttCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCGCt                      >  1:317231/1‑69 (MQ=255)
                                         tGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCGCTTTcc                  <  1:246495/68‑1 (MQ=255)
                                         tGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCGCTTTcc                  <  1:688593/68‑1 (MQ=255)
                                         tGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCGCTTTcc                  <  1:823230/68‑1 (MQ=255)
                                         tGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCGCTTTcc                  <  1:852092/68‑1 (MQ=255)
                                                         agagaCATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCGCTTTCCAGTTCCGCTTTAcgcg  <  1:2329357/68‑1 (MQ=255)
                                                         |                                                                   
GTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGACGCGCATTCAGTGCAGCGCTTTCCAGTTCCGCTTTACGCG  >  minE/1126292‑1126416

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: