Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,172,268 T→C 48.3% K1092E (AAG→GAG)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,172,2680TC48.3% ‑5.1 / 40.1 29K1092E (AAG→GAG) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (9/6);  new base C (6/8);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTT  >  W3110S.gb/1172206‑1172327
                                                              |                                                           
gCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCgg                                                        >  1:667295/1‑68 (MQ=255)
      aTGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTc                                                 >  1:1335222/1‑69 (MQ=255)
        gctgcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCg                                                         >  1:3092/1‑59 (MQ=255)
        gctgcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCg                                                         >  1:2531076/1‑59 (MQ=255)
        gctgcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCg                                                         >  1:2052584/1‑59 (MQ=255)
           gcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                            >  1:2826545/1‑69 (MQ=255)
           gcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                            >  1:3108053/1‑69 (MQ=255)
              gCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCAc                                           <  1:709255/67‑1 (MQ=255)
                  tttttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCACCCCGcc                                     >  1:3127759/1‑69 (MQ=255)
                    tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCg                                                >  1:2561289/1‑56 (MQ=255)
                    tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCg                                                >  1:3425654/1‑56 (MQ=255)
                    tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCtt                                   >  1:964283/1‑69 (MQ=255)
                    tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCACCCCGCCt                                    <  1:1714862/68‑1 (MQ=255)
                    tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCACCCCGCCt                                    <  1:2639120/68‑1 (MQ=255)
                       gcagcaAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCACCCCGCCTTtc                                 <  1:3246486/68‑1 (MQ=255)
                         agcaAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCACCCCGCCTTtctc                               <  1:1559790/68‑1 (MQ=255)
                         agcaAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCACCCCGCCTTtctc                               <  1:1334485/68‑1 (MQ=255)
                            aaaCCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCg                                                <  1:3488525/48‑1 (MQ=255)
                                    cAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCACCCCGc                                      >  1:1676330/1‑50 (MQ=255)
                                    cAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCACCCCGCCTTtct                                >  1:2082878/1‑56 (MQ=255)
                                       ccAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCtt                <  1:2917773/69‑1 (MQ=255)
                                       ccAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCtt                <  1:234871/69‑1 (MQ=255)
                                       ccAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATCACCTTCCAGCtt                <  1:59910/69‑1 (MQ=255)
                                        cAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCAc                                           <  1:3013320/41‑1 (MQ=255)
                                          ggCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCt             >  1:145853/1‑69 (MQ=255)
                                                  tAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCa     >  1:2371477/1‑69 (MQ=255)
                                                      gTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGtt  <  1:3351488/68‑1 (MQ=255)
                                                      gTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGtt  <  1:2755109/68‑1 (MQ=255)
                                                      gTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGtt  <  1:3484076/68‑1 (MQ=255)
                                                              |                                                           
GCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTT  >  W3110S.gb/1172206‑1172327

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: