Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,495,966:1 +T 73.2% coding (820/924 nt) ydcI ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,495,9661.T73.2% 78.1 / 37.3 41coding (820/924 nt)ydcIpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (8/3);  new base T (20/10);  total (28/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

AACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAG  >  W3110S.gb/1495899‑1496025
                                                                    |                                                           
aaCTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGt                                                            >  1:1025069/1‑69 (MQ=255)
aaCTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGt                                                            >  1:48457/1‑69 (MQ=255)
  cTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGtt                                                           <  1:2286399/68‑2 (MQ=255)
    ggCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTtcc                                                         >  1:810936/1‑68 (MQ=255)
    ggCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTtcc                                                         >  1:2009459/1‑68 (MQ=255)
    ggCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTtcc                                                         >  1:1934173/1‑68 (MQ=255)
          cAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAAca                                                   <  1:1469807/69‑1 (MQ=255)
          cAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAAca                                                   >  1:1969900/1‑69 (MQ=255)
             aagaGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTcc                                                          <  1:2459518/58‑1 (MQ=255)
               gagaACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAAcaggcag                                              >  1:2576584/1‑69 (MQ=255)
               gagaACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAAcaggcag                                              >  1:2110125/1‑69 (MQ=255)
               gagaACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAAcaggcag                                              >  1:278700/1‑69 (MQ=255)
                   aCGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGcc                                          <  1:549496/69‑1 (MQ=255)
                        gCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGttc                                                          >  1:539062/1‑45 (MQ=255)
                        gCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGttc                                                          >  1:211089/1‑45 (MQ=255)
                             tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGc                               >  1:3700025/1‑69 (MQ=255)
                                  gggTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTgg                                <  1:1464954/64‑1 (MQ=255)
                                         tATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAgg                     >  1:3458923/1‑68 (MQ=255)
                                         tATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgt                    >  1:351065/1‑69 (MQ=255)
                                         tATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgt                    >  1:3030542/1‑69 (MQ=255)
                                         tATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgt                    >  1:3606897/1‑69 (MQ=255)
                                         tATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgt                    >  1:1648064/1‑69 (MQ=255)
                                         tATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtc                   <  1:2439149/69‑1 (MQ=255)
                                         tATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtc                   <  1:1270157/69‑1 (MQ=255)
                                           ttCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc                        >  1:1372728/1‑63 (MQ=255)
                                             ccAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcgt                 <  1:3386295/68‑1 (MQ=255)
                                             ccAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcgt                 <  1:1341854/68‑1 (MQ=255)
                                             ccAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCt               >  1:40751/1‑69 (MQ=255)
                                             ccAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCt               >  1:3277131/1‑69 (MQ=255)
                                                 tCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcg                  >  1:969030/1‑63 (MQ=255)
                                                  cGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt             >  1:1024555/1‑67 (MQ=255)
                                                   ggTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTAc          >  1:732572/1‑69 (MQ=255)
                                                   ggTTCGCCGGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAAGACGCAGgtcg                  >  1:2046036/1‑60 (MQ=255)
                                                      tCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAgc       >  1:2345050/1‑69 (MQ=255)
                                                      tCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAgc       >  1:2299136/1‑69 (MQ=255)
                                                      tCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAgc       >  1:2193757/1‑69 (MQ=255)
                                                      tCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAgc       >  1:1547428/1‑69 (MQ=255)
                                                           tGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAg  <  1:1165793/69‑1 (MQ=255)
                                                           tGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAg  <  1:957238/69‑1 (MQ=255)
                                                             cACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcgtc                >  1:2510581/1‑53 (MQ=255)
                                                             cACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcgtc                >  1:494087/1‑53 (MQ=255)
                                                             cACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCtttt            <  1:2400207/57‑1 (MQ=255)
                                                             cACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCtttt            <  1:2078216/57‑1 (MQ=255)
                                                             cACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCtttt            <  1:611518/57‑1 (MQ=255)
                                                             cACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt             >  1:213015/1‑55 (MQ=255)
                                                               ccATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCtttt            >  1:1556941/1‑55 (MQ=255)
                                                                    |                                                           
AACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAG  >  W3110S.gb/1495899‑1496025

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: