Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,000,025 A→T V94V (GTT→GTA)  uspA ← universal stress global response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,000,0250AT75.9% 275.3 / 79.2 145V94V (GTT→GTA) uspAuniversal stress global response regulator
Reads supporting (aligned to +/- strand):  ref base A (0/35);  new base T (0/110);  total (0/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

CAAACCACCAAATCCATATCGTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCA  >  W3110S.gb/3999981‑4000061
                                            |                                    
caAACCACCAGATCCATATTGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:559190/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGTCCCAGATcgcc                     <  1:134846/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:476980/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:600586/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:58621/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:569270/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:531008/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:523098/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:518579/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:505558/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:49375/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:608772/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:41589/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:412596/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:387238/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:375861/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:353383/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:259456/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:235718/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:209873/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:194536/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:194083/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:830670/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:988860/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:988638/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:954888/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:947208/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:924730/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:923303/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:89656/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:895985/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:84480/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:833255/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1005586/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:760661/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:730215/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:726200/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:71566/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:691489/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:676013/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:65310/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:637947/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:632583/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:622160/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1162720/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1297982/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:128875/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1270008/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1265000/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1242656/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1237104/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1198324/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1189816/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1184403/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1177848/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:178416/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1159866/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:114733/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1137253/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1129517/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:112754/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1118435/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1074500/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1054957/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1003373/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1317536/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:183220/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1451638/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1445584/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1436952/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1421803/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1406650/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1370325/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1350437/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:134297/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1352353/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1332481/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATagcc                     <  1:1160411/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTCGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:354177/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCCTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:658657/62‑1 (MQ=255)
caAACCACCACATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1356928/62‑1 (MQ=255)
 aaaCCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:618534/61‑1 (MQ=255)
           aTCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:607310/51‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCTAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:811313/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGGTGCCGCTCAGGGTTTCa  <  1:757783/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1016738/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1068313/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1313066/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:237124/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:880810/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:865688/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:285660/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1147437/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:801480/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1148020/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:349845/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1418493/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1406732/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1230383/61‑1 (MQ=255)
                     tATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:220410/60‑1 (MQ=255)
                                            |                                    
CAAACCACCAAATCCATATCGTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCA  >  W3110S.gb/3999981‑4000061

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