Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,346,506 C→T E651K (GAA→AAA)  gmr ← modulator of Rnase II stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,346,5060CT76.0% 41.7 / 12.6 25E651K (GAA→AAA) gmrmodulator of Rnase II stability
Reads supporting (aligned to +/- strand):  ref base C (6/0);  new base T (19/0);  total (25/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

TAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTG  >  W3110S.gb/1346450‑1346531
                                                        |                         
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:1210150/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:842566/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:814891/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:685535/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:623870/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:603543/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:440654/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:391244/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:249946/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:2279914/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:2031857/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:1943244/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:1838519/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:1742036/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:164364/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:1546166/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:1442646/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcgg                    >  1:1003502/1‑61 (MQ=255)
tAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTTGAAgcg                     >  1:2479331/1‑61 (MQ=255)
                  tACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTg  >  1:1937749/1‑64 (MQ=255)
                  tACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTg  >  1:2336855/1‑64 (MQ=255)
                  tACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTg  >  1:1746329/1‑64 (MQ=255)
                  tACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTg  >  1:1462556/1‑64 (MQ=255)
                  tACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTg  >  1:702827/1‑64 (MQ=255)
                  tACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTg  >  1:913711/1‑64 (MQ=255)
                                                        |                         
TAAGTATATATATTCATCTACTTATGCGCGCTTCAGATAGCGTTTATACCAGCGTTCGAAGGCGACGGCGGGCATCGGTTTG  >  W3110S.gb/1346450‑1346531

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: ‑