Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,000,037 G→A D90D (GAC→GAT 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,0370GA75.3% 198.8 / 65.0 146D90D (GAC→GATuspAuniversal stress global response regulator
Reads supporting (aligned to +/- strand):  ref base G (0/36);  new base A (0/110);  total (0/146)
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.31e-01

Only 100 of 146 total aligned reads displayed.

CAAACCACCAAATCCATATCGTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCA  >  W3110S.gb/3999981‑4000061
                                                        |                        
caAACCACGAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:568829/62‑1 (MQ=255)
caAACCACCAGATCCATATTGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:559190/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGTCCCAGATcgcc                     <  1:134846/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:351634/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:528383/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:523098/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:505558/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:49375/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:41589/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:398461/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:387238/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:58621/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:271971/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:259456/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:230917/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:198170/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:194536/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:188251/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:183220/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1470359/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1451638/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:760661/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:833255/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:830670/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1005586/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:730215/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:726200/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:691489/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:652514/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:608772/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:595937/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1159866/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:128875/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1270008/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1256197/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1237104/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1436952/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:1162720/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1198324/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:1370325/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1353691/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1352353/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1328152/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1332481/62‑1 (MQ=255)
caAACCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1370494/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)
 aaaCCACCAGATCCATATCGTATTTCTTGATTGCATCGACCAGTACCTGGCCCAGATcgcc                     <  1:1349064/61‑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:1406732/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:880810/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1411679/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1068313/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:865688/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1418493/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:834586/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1016738/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1424550/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:237124/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:704571/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:801480/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1147689/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1148020/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1304183/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:66845/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:349845/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:356879/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1230383/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1313066/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1274803/61‑1 (MQ=255)
                    gTATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:428971/61‑1 (MQ=255)
                     tATTTCTTGATTGCATCGACCAGAACCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:220410/60‑1 (MQ=255)
                          cTTGATTGCATCGACCTGTACCTGGCCCAGATcgcc                     <  1:509627/36‑1 (MQ=25)
                                             aCCTGGCCCAGGTCGCCGCTGCCGCTCAGGGTTTCa  <  1:1256299/36‑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: