Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,886,205 A→G 40.6% intergenic (‑330/+29) ygcB ← / ← cysH conserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,886,2050AG40.6% 13.7 / 35.9 32intergenic (‑330/+29)ygcB/cysHconserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base A (12/7);  new base G (6/7);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

CTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCA  >  W3110S.gb/2886136‑2886270
                                                                     |                                                                 
cttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGa                                                                   >  1:3725242/1‑70 (MQ=255)
              cACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGa                                                                   >  1:2426940/1‑56 (MQ=255)
              cACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGGCAGAGGCGTGTAAtt                                                    >  1:806943/1‑71 (MQ=255)
                cGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTg                                                  <  1:1103369/71‑1 (MQ=255)
                cGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTg                                                  >  1:121167/1‑71 (MQ=255)
                cGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTAAATTTg                                                  <  1:662/71‑1 (MQ=255)
                    tAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGGCAGAGGCGTTTAATTTGTCCg                                              >  1:3842591/1‑71 (MQ=255)
                         tCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGGCAGAgg                                                             <  1:3894747/51‑1 (MQ=255)
                          cAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAata                                        >  1:804581/1‑71 (MQ=255)
                               aTGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTa                                    <  1:2687511/70‑1 (MQ=255)
                                 gctgctCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAcc                                  >  1:2677358/1‑70 (MQ=255)
                                  ctgctCGCAAAATAAAAACGCCCCCATCATTTCTGGCAGAGGCGttt                                                        <  1:1089104/47‑1 (MQ=255)
                                  ctgctCGCAAAATAAAAACGCCCCCATCATTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                >  1:2183547/1‑71 (MQ=255)
                                  ctgctCGCAAAATAAAAACGCCCCCATCATTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                >  1:505161/1‑71 (MQ=255)
                                  ctgctCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                >  1:2882473/1‑71 (MQ=255)
                                  ctgctCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                >  1:1714354/1‑71 (MQ=255)
                                    gctCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAggg                                                            >  1:3338060/1‑40 (MQ=255)
                                        ccAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAat                                         <  1:464834/55‑1 (MQ=255)
                                         cAAAATAAACACGCCCCCATCATTTCTGGCAGAGGCGTTTAAt                                                     <  1:3941824/43‑1 (MQ=255)
                                         cAAAATAAAAACGCCCCCATCATTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtgt                          <  1:1648281/70‑1 (MQ=255)
                                          aaaaTAAAAACGCCCCCATCATTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTa                         <  1:3149195/70‑1 (MQ=255)
                                           aaaTAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTa                         <  1:2096858/69‑1 (MQ=255)
                                              tAAAAACGCCCCCATCATTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCca                    <  1:1286966/71‑1 (MQ=255)
                                                        cccATCATTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATAttt                                     <  1:1537212/44‑1 (MQ=255)
                                                          catcatTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                   >  1:999879/1‑44 (MQ=255)
                                                          catcatTTCTGACAGAGGCGTTTAATTTTTCCGGCAata                                        >  1:410782/1‑39 (MQ=255)
                                                           atcatTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTAGTTTAACCCACATTCCCtttt         >  1:2235511/1‑69 (MQ=255)
                                                            tcatTTCTGACAGAGGCGTTTAATTTGTCCGGCAat                                         <  1:4070232/36‑1 (MQ=255)
                                                            tcatTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTcc              <  1:3056810/63‑1 (MQ=255)
                                                              atTTCTGACAGAGGCGTTTAATTTGTCCGGCAatat                                       >  1:3898731/1‑36 (MQ=255)
                                                              atTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTa                                    >  1:3619376/1‑39 (MQ=255)
                                                              atTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAg     >  1:3885651/1‑70 (MQ=255)
                                                                ttCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGcc   >  1:3960500/1‑70 (MQ=255)
                                                                 tCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCa  >  1:819058/1‑70 (MQ=255)
                                                                     |                                                                 
CTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCA  >  W3110S.gb/2886136‑2886270

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: