Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,989,743 (G)7→6 intergenic (‑617/+12) yhiS ← / ← arsC hypothetical protein/arsenate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,989,7370G.81.8% 47.2 / 14.0 22intergenic (‑611/+18)yhiS/arsChypothetical protein/arsenate reductase
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base . (5/13);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGGATAAGACGCGCCAGCGTCGCATCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTTCTCGCCATCTTC  >  W3110S.gb/3989673‑3989802
                                                                |                                                                 
ggTCGGATAAGACGCGCCAGCGTCGCATCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTgggggga                                                             >  1:2745512/1‑70 (MQ=255)
  tCGGATAAGACGCGCCAGCGTCGCATCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTggggggat                                                            <  1:1860328/70‑3 (MQ=255)
   cGGATAAGACGCGCCAGCGTCGCATCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGGata                                                          >  1:4077969/1‑71 (MQ=255)
   cGGATAAGACGCGCCAGCGTCGCATCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGGata                                                          >  1:1389939/1‑71 (MQ=255)
   cGGATAAGACGCGCCAGCGTCGCATCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGGata                                                          >  1:3058315/1‑71 (MQ=255)
             gcgcCAGCGTCGCATCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGGATATCGCCGCtt                                                 <  1:1788310/70‑1 (MQ=255)
             gcgcCAGCGTCGCATCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTa                                               <  1:2482940/71‑1 (MQ=255)
                          aTCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTa                                  >  1:890559/1‑71 (MQ=255)
                          aTCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTa                                  >  1:1437080/1‑71 (MQ=255)
                               aCATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGc                             <  1:251991/71‑1 (MQ=255)
                               aCATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGc                             <  1:3218420/71‑1 (MQ=255)
                               aCATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGc                             <  1:3005152/71‑1 (MQ=255)
                               aCATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGc                             <  1:2797247/71‑1 (MQ=255)
                               aCATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGc                             <  1:1121186/71‑1 (MQ=255)
                                 ataaataaAGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCtt                           <  1:2946369/71‑1 (MQ=255)
                                     ataaaGCGCACTTTTCTAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGc                             <  1:4077553/65‑1 (MQ=255)
                                                     tAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTtctc           <  1:584504/67‑1 (MQ=255)
                                                     tAACAACCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTtctc           <  1:1285429/67‑1 (MQ=255)
                                                       acaacCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTtctc           <  1:2451526/65‑1 (MQ=255)
                                                       acaacCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTTCTCGc         >  1:1317548/1‑67 (MQ=255)
                                                       acaacCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTTCTCGCCa       <  1:1767683/69‑1 (MQ=255)
                                                       acaacCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTTCTCGCCa       <  1:1539704/69‑1 (MQ=255)
                                                          acCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTTCTCGCCATCTTc  >  1:2394646/1‑71 (MQ=255)
                                                          acCTGTGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTTCTCGCCATCTTc  >  1:46749/1‑71 (MQ=255)
                                                                |                                                                 
GGTCGGATAAGACGCGCCAGCGTCGCATCCGACATAAATAAAGCGCACTTTTCTAACAACCTGTGGGGGGGATATCGCCGCTTTATTTCAGGCGCTTACCCGCTTCATCAACCACTTTCTCGCCATCTTC  >  W3110S.gb/3989673‑3989802

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: