Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,886,205 A→G 17.4% 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,2050AG17.4% 45.9 / 5.4 23intergenic (‑330/+29)ygcB/cysHconserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base A (14/5);  new base G (2/2);  total (16/7)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

CTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGA  >  W3110S.gb/2886136‑2886274
                                                                     |                                                                     
cttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGa                                                                       >  1:454966/1‑70 (MQ=255)
                cGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGGCAGAGGCGTTTAATTTg                                                      >  1:1867762/1‑71 (MQ=255)
                  tttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTc                                                    >  1:21804/1‑71 (MQ=255)
                   ttAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTcc                                                   >  1:496601/1‑71 (MQ=255)
                        aTCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGGCAGAGGCGTTTAATTTGTCCGGCaa                                              <  1:3842202/71‑1 (MQ=255)
                             gAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGa                                                                       >  1:4014199/1‑41 (MQ=255)
                                  ctgctCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                    >  1:196715/1‑71 (MQ=255)
                                      tCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTc                                  <  1:4542581/69‑1 (MQ=255)
                                          aaaaTAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAat                                             <  1:2472644/54‑1 (MQ=255)
                                             aTAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCa                                               >  1:1469177/1‑49 (MQ=255)
                                                     gCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAat                                             <  1:2467253/43‑1 (MQ=255)
                                                        cccATCATTTCTAACAGAGGCGTTTAATTTGTCCGGCAATAttt                                         <  1:3944949/44‑1 (MQ=255)
                                                          catcatTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                       >  1:717058/1‑44 (MQ=255)
                                                            tcatTTCTGACAGAGGCGTTTAATTTGTCCGGCAATAtt                                          <  1:3820293/39‑1 (MQ=255)
                                                              atTTCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGc        >  1:1076977/1‑71 (MQ=255)
                                                              atTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAc                                       >  1:4511878/1‑40 (MQ=255)
                                                              atTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtgt                              >  1:4008992/1‑49 (MQ=255)
                                                              atTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAtt                    >  1:1036508/1‑59 (MQ=255)
                                                              atTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGc        >  1:3913269/1‑71 (MQ=255)
                                                              atTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGc        >  1:4520072/1‑71 (MQ=255)
                                                              atTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGc        >  1:4605929/1‑71 (MQ=255)
                                                               tttCTGGCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGcc       <  1:1583728/71‑1 (MQ=255)
                                                                 tCTGACAGAGGCGTTTAATTTGTCCTGCAATATTTACCCTTCgt                                >  1:2211520/1‑44 (MQ=255)
                                                                 tCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgt                                >  1:4478861/1‑44 (MQ=255)
                                                                  cTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCa      >  1:1317160/1‑69 (MQ=255)
                                                                     aCAGAGGCGTTTAATTTGTCCGGCGATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGa  <  1:1536127/70‑1 (MQ=255)
                                                                     aCAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGa  <  1:3685983/70‑1 (MQ=255)
                                                                     |                                                                     
CTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGA  >  W3110S.gb/2886136‑2886274

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: