Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 640,175 G→A L400L (CTG→CTA ahpF → alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb640,1750GA83.3% 139.9 / 16.9 78L400L (CTG→CTAahpFalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/13);  new base A (0/65);  total (0/78)
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.95e-01

AAATGAAAGCCGACCAGGTTCTGCAGGACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGT  >  W3110S.gb/640141‑640225
                                  |                                                  
aaaTTAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCTACattat                          <  1:1314130/61‑1 (MQ=255)
aaaTGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:660770/61‑1 (MQ=255)
aaaTGAAACCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:228058/61‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:470178/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:553550/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:545246/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:541419/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:510800/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:482910/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:474141/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:556302/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:46645/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:427877/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:406845/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:368310/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:330177/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:324222/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:29249/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1002629/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:572328/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:644796/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:704816/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:718690/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:737896/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:744231/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:786360/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:788368/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:793212/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:922140/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:928064/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:943545/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:949519/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:961783/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:971779/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1215864/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1367560/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1347708/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1327612/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:265939/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:125219/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1250169/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1250103/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:124560/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1220215/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1283212/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:118042/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:11760/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1168488/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1122715/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:109018/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1064018/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1010544/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1000327/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:260216/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1437273/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1440488/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1447529/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:249229/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:200701/58‑1 (MQ=255)
   tGAAAGCCGACCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:1459408/58‑1 (MQ=255)
            aCCAGGTTCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:315923/49‑1 (MQ=255)
                ggTTCTGCAGGACAATCTACGCAGGCTGAAAAACGTCGACattat                          <  1:588566/45‑1 (MQ=255)
                  ttCTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:27562/43‑1 (MQ=255)
                    cTGCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:301569/41‑1 (MQ=255)
                      gCAGGACAAACTACGCAGCCTGAAAAACGTCGACattat                          <  1:704559/39‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:677440/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:444293/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:172335/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:4638/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:355908/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:21988/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:587781/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:1157321/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:58051/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:1122005/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:1433536/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:285528/60‑1 (MQ=255)
                         ggACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGt  <  1:1406640/60‑1 (MQ=255)
                                  |                                                  
AAATGAAAGCCGACCAGGTTCTGCAGGACAAACTGCGCAGCCTGAAAAACGTCGACATTATTCTGAATGCGCAAACCACGGAAGT  >  W3110S.gb/640141‑640225

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: