Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,726,854 T→C R190R (CGT→CGC pstC → phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,726,8540TC75.4% 106.6 / 29.5 61R190R (CGT→CGCpstCphosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (0/15);  new base C (46/0);  total (46/15)
Fisher's exact test for biased strand distribution p-value = 1.42e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCATCATGATTATTCCGTACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCGG  >  W3110S.gb/3726814‑3726894
                                        |                                        
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAAcccc                      >  1:467617/1‑61 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2403897/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1924410/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:464419/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2132095/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2193854/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2216131/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:38030/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2236491/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:354675/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2268969/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2398544/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:903372/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2480572/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:332819/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2539267/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:254491/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2554861/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2602307/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2787044/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2792104/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2854933/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:2867021/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:314329/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1833312/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:887152/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1197247/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:840822/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1332728/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1339519/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1365881/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1397824/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1452019/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:14735/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1488706/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:83259/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1571038/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1584321/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:508960/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:481196/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1710830/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1767730/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1776457/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAACAAACCCCg                     >  1:1829973/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATTCGAAAAAACCCCg                     >  1:1914358/1‑62 (MQ=255)
tCGCCATCATGATTATTCCGTACATTGCAGCGGTAATGCGCGATGTATACGAAAAAAccc                       >  1:1024131/1‑60 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTcag  <  1:1585438/61‑3 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:253329/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:2259217/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:2233612/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:409664/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:1961214/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:1795220/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:165388/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:57615/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:609550/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:1491982/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:124949/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:100148/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgg  <  1:977952/61‑1 (MQ=255)
                    tACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCgc  <  1:104882/61‑2 (MQ=255)
                                        |                                        
TCGCCATCATGATTATTCCGTACATTGCGGCGGTAATGCGTGATGTGTTCGAACAAACCCCGGTGATGATGAAAGAGTCGG  >  W3110S.gb/3726814‑3726894

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: