Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,172,268 T→C 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,2680TC86.4% 117.5 / 10.1 44K1092E (AAG→GAG) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base C (20/18);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

ATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCT  >  W3110S.gb/1172212‑1172330
                                                        |                                                              
aTGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCg                                                   <  1:889531/70‑1 (MQ=255)
aTGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCg                                                   <  1:1204297/70‑1 (MQ=255)
aTGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCg                                                   <  1:3611895/70‑1 (MQ=255)
aTGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCg                                                   <  1:2529005/70‑1 (MQ=255)
  gctgcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCg                                                   <  1:2053378/68‑1 (MQ=255)
   ctgcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                               <  1:1203460/71‑1 (MQ=255)
   ctgcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                               <  1:2566830/71‑1 (MQ=255)
      cGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAAt                                                      <  1:1198123/61‑1 (MQ=255)
      cGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAAt                                                      <  1:5828/61‑1 (MQ=255)
      cGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAAt                                                      <  1:4138585/61‑1 (MQ=255)
      cGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAAt                                                      <  1:2806077/61‑1 (MQ=255)
      cGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAAt                                                      <  1:2615454/61‑1 (MQ=255)
      cGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAAt                                                      <  1:1376318/61‑1 (MQ=255)
      cGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCCGCAAAt                                                      <  1:977/61‑1 (MQ=255)
       ggCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAAt                                                      <  1:1565212/60‑1 (MQ=255)
            tttttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCg                                                   <  1:389113/58‑1 (MQ=255)
            tttttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCg                                                   <  1:3577553/58‑1 (MQ=255)
            tttttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCg                                                   <  1:1276817/58‑1 (MQ=255)
            tttttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCg                                                   <  1:2928743/58‑1 (MQ=255)
            tttttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCg                                                   <  1:424210/58‑1 (MQ=255)
            tttttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCg                                                   <  1:2593283/58‑1 (MQ=255)
            tttttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCg                                                   <  1:763475/58‑1 (MQ=255)
            tttttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCg                                                   <  1:1093950/58‑1 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCTATCa                                               >  1:3316452/1‑60 (MQ=39)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:753403/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:875839/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:90771/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:4127529/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:958347/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:1878558/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:132009/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:1512699/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:1520437/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:2304538/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:2141375/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:2077998/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGAt                                                 >  1:1681144/1‑58 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                               >  1:3968544/1‑60 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                               >  1:3387034/1‑60 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                               >  1:1243671/1‑60 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                               >  1:2574669/1‑60 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                               >  1:1655604/1‑60 (MQ=255)
              tttGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTCCTCGGCAAATTCGATCa                                               >  1:978612/1‑60 (MQ=255)
                                                gTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCt  >  1:3109265/1‑71 (MQ=255)
                                                        |                                                              
ATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCT  >  W3110S.gb/1172212‑1172330

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: