Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,566,298 G→A L384L (CTG→TTG)  dos ← cAMP phosphodiesterase, heme‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,566,2980GA83.1% 184.6 / 28.0 83L384L (CTG→TTG) doscAMP phosphodiesterase, heme‑regulated
Reads supporting (aligned to +/- strand):  ref base G (14/0);  new base A (0/69);  total (14/69)
Fisher's exact test for biased strand distribution p-value = 3.78e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTG  >  W3110S.gb/1566243‑1566355
                                                       |                                                         
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2394894/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2749859/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2737476/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2679268/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:267044/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2549910/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2531762/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2500324/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2468186/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:279435/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2354002/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:234605/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2169332/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2118354/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1018218/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2035960/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1911079/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1855009/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:722052/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:973851/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:92790/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:88904/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:858608/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:839105/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:753507/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:746383/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:728161/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2054658/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:685097/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:642236/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:52550/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:51613/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:45450/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:414942/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2824520/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:2798407/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1208783/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1587256/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1585329/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1850352/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1514280/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1512883/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1454699/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:133339/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1248643/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1567032/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:115058/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1143757/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1140652/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1109419/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1108104/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1103135/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1050359/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1034099/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1751410/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1736185/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:160168/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1625729/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1631341/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1646085/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1677728/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1673169/62‑1 (MQ=255)
aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1650990/62‑1 (MQ=255)
aCGGCTTTGTAGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1810219/62‑1 (MQ=255)
 cGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:519564/61‑1 (MQ=255)
 cGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1194026/61‑1 (MQ=255)
 cGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1600432/61‑1 (MQ=255)
 cGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1522553/61‑1 (MQ=255)
                      cGAGGTAATTGTGCAGGTTATTGCGATTTGGCAAACCGGt                                                     <  1:1703270/40‑1 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCt   >  1:2556552/1‑61 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:1298717/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:2463160/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:524498/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:401317/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:334804/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:1149852/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:2719070/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:2098934/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:776904/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:2692049/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:1759432/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:1972580/1‑62 (MQ=255)
                                                   ggCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg  >  1:1965326/1‑62 (MQ=255)
                                                       |                                                         
ACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTG  >  W3110S.gb/1566243‑1566355

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