Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,175,236 C→A 66.6% Q102H (CAG→CAT mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,175,2360CA66.6% 22.2 / 59.7 51Q102H (CAG→CATmfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base C (9/8);  new base A (12/22);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAGAAAACTGTGTGGGCAAACGCGCTGCATAAGCGTATTCACCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGAAAAACT  >  W3110S.gb/1175171‑1175304
                                                                 |                                                                    
tGGAGAAAACTGTGTGGGCAAACGCGCTGCATAAGCGTATTCACCGGAACAATCAGTACGCCACGCTGc                                                                   <  1:2918063/69‑1 (MQ=255)
tGGAGAAAACTGTGTGGGCAAACGCGCTGCATAAGCGTATTCACCGGAACAATCAGTACGCCACGCTGc                                                                   <  1:1919938/69‑1 (MQ=255)
              tgGGCAAACGCGCTGCATAAGCGTATTCACCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGc                                                       <  1:518857/67‑1 (MQ=255)
                ggCAAACGCGCTGCATAAGCGTATTCACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTgg                                                    >  1:3556827/1‑68 (MQ=255)
                         gcTGCATAAGCGTATTCACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTgg                                          <  1:2259623/69‑1 (MQ=255)
                         gcTGCATAAGCGTATTCACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTgg                                          <  1:641317/69‑1 (MQ=255)
                          cTGCATAAGCGTATTCACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                         <  1:1956367/69‑1 (MQ=255)
                                       ttCACCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCgaggag                            <  1:1228873/69‑1 (MQ=255)
                                       ttCACCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCgaggag                            <  1:1477340/69‑1 (MQ=255)
                                       ttCACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTaa                                                 >  1:1149317/1‑48 (MQ=255)
                                       ttCACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTaa                                                 <  1:3061003/48‑1 (MQ=255)
                                       ttCACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTaa                                                 >  1:515101/1‑48 (MQ=255)
                                       ttCACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTaa                                                 >  1:580988/1‑48 (MQ=255)
                                       ttCACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTaa                                                 >  1:1522545/1‑48 (MQ=255)
                                         cACCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGgaga                           >  1:3648638/1‑68 (MQ=255)
                                         cACCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGgaga                           >  1:2488541/1‑68 (MQ=255)
                                         cACCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                          <  1:653356/69‑1 (MQ=255)
                                         cACCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                          >  1:29641/1‑69 (MQ=255)
                                         cACCGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGAt                          <  1:2849368/69‑1 (MQ=255)
                                          aCCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATa                         >  1:770408/1‑69 (MQ=255)
                                          aCCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATa                         >  1:2908281/1‑69 (MQ=255)
                                           ccGGAACAATCAGTACGCCACGATGCATCGTCGGTAGc                                                       <  1:638866/38‑1 (MQ=255)
                                           ccGGAACAATCAGTACGCCACGATGCATCGTCGGTAGc                                                       <  1:640186/38‑1 (MQ=255)
                                           ccGGAACAATCAGTACGCCACGATGCATCGTCGGTAGc                                                       <  1:479578/38‑1 (MQ=255)
                                           ccGGAACAATCAGTACGCCACGATGCATCGTCGGTAGc                                                       <  1:744106/38‑1 (MQ=255)
                                           ccGGAACAATCAGTACGCCACGATGCATCGTCGGTAGc                                                       <  1:904913/38‑1 (MQ=255)
                                           ccGGAACAATCAGTACGCCACGATGCATCGTCGGTAGc                                                       <  1:1953410/38‑1 (MQ=255)
                                            cGGAACAATCAGTACGCCACGATGCATCGTCGGTAGc                                                       <  1:3325447/37‑1 (MQ=255)
                                            cGGAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAAt                       >  1:204862/1‑69 (MQ=255)
                                             ggAACAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATg                      <  1:900361/69‑1 (MQ=255)
                                                aCAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGt                                                   >  1:1838571/1‑37 (MQ=255)
                                                aCAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTc                    >  1:600253/1‑68 (MQ=255)
                                                aCAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTc                    >  1:708467/1‑68 (MQ=255)
                                                aCAATCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTc                    >  1:1619696/1‑68 (MQ=255)
                                                   aTCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCtg                 >  1:950588/1‑68 (MQ=255)
                                                    tCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCgag                               >  1:586819/1‑53 (MQ=255)
                                                    tCAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCgag                               >  1:2126173/1‑53 (MQ=255)
                                                     cAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                         <  1:614199/42‑1 (MQ=255)
                                                     cAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                         <  1:1805028/42‑1 (MQ=255)
                                                     cAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                         <  1:2510895/42‑1 (MQ=255)
                                                     cAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                         <  1:1595798/42‑1 (MQ=255)
                                                     cAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGa                                         <  1:889510/42‑1 (MQ=255)
                                                     cAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCt                  <  1:2089628/65‑1 (MQ=255)
                                                     cAGTACGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCt                  <  1:826255/65‑1 (MQ=255)
                                                        tACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAg            >  1:951161/1‑68 (MQ=255)
                                                         aCGCCACCATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAgg                                     <  1:2133059/42‑1 (MQ=38)
                                                          cGCCACGATGCATCGTCGGTAGCTGGTAAAGGGTGGAAAgg                                     <  1:1821845/41‑1 (MQ=255)
                                                            ccACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCtgatga             <  1:2940602/63‑1 (MQ=255)
                                                             cACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGaaa      <  1:1116816/69‑1 (MQ=255)
                                                                gCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGAAAAAc   >  1:2962982/1‑69 (MQ=255)
                                                                gCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGAAAAAc   >  1:200272/1‑69 (MQ=255)
                                                                 cTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGAAAAACt  >  1:142318/1‑69 (MQ=255)
                                                                 cTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGAAAAACt  >  1:959989/1‑69 (MQ=255)
                                                                 |                                                                    
TGGAGAAAACTGTGTGGGCAAACGCGCTGCATAAGCGTATTCACCGGAACAATCAGTACGCCACGCTGCATCGTCGGTAGCTGGTAAAGGGTGGAAAGGCGCGAGGAGATAATGTCCTGATGAGGCGAAAAACT  >  W3110S.gb/1175171‑1175304

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: