Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 950,876 C→A W343C (TGG→TGT)  yddV ← predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE950,8760CA84.0% 370.7 / 80.8 169W343C (TGG→TGT) yddVpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/27);  new base A (0/142);  total (0/169)
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.80e-01

Only 100 of 169 total aligned reads displayed.

ATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTT  >  minE/950835‑950906
                                         |                              
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTcta                      <  1:771350/52‑3 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:595630/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:701940/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:683524/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:679927/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:669274/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:647011/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:629647/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:624486/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:617560/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:610276/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:71844/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:586267/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:569789/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:543154/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:528091/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:517435/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:502100/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:486951/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:480562/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:47384/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:452851/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:823384/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:996140/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:959451/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:94921/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:937243/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:919483/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:917981/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:850281/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:84967/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:83383/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:823530/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:705281/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:797348/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:760254/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:75511/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:749971/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:738255/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:728476/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:726364/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:719481/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1023594/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1195262/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:251916/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:23694/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:226179/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:208849/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:19144/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:179140/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:178561/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:163315/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:134678/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:129615/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:438467/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1178762/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1167882/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1158469/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1145154/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1140600/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:113288/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1104695/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1055227/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1047357/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:10056/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:386457/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:447848/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:426348/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:412750/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:410656/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:406092/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:376194/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:37594/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:368568/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:350727/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:3105/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:336400/52‑1 (MQ=255)
aTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:315909/52‑1 (MQ=255)
aTGATTACCTCAGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:708095/52‑1 (MQ=39)
 tGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTTaa                      <  1:660748/51‑1 (MQ=255)
 tGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:175963/51‑1 (MQ=255)
 tGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:484914/51‑1 (MQ=255)
 tGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:491281/51‑1 (MQ=255)
 tGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1150642/51‑1 (MQ=255)
   aTTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:881990/49‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACTCGGTTTCTACCTCGtttt  <  1:304498/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:316090/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:488004/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:507354/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:784533/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:927042/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:599393/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:105392/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:966967/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:974580/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCtt  <  1:348590/62‑1 (MQ=255)
          ctGCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:24702/42‑1 (MQ=255)
            gCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:1184077/40‑1 (MQ=255)
            gCATCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:940417/40‑1 (MQ=255)
               tCCTGGCGCATCTAAAGACTGGCTTTACAGAGTTCaa                      <  1:845373/37‑1 (MQ=255)
                                         |                              
ATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTT  >  minE/950835‑950906

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