Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,785,165:1 +T 12.8% coding (475/1701 nt) ydiD → short chain acyl‑CoA synthetase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,785,1651.T12.8% 152.9 / 24.7 47coding (475/1701 nt)ydiDshort chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base . (26/15);  new base T (3/3);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AATAAGTGTCAGGCAAAAATGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCA  >  W3110S.gb/1785099‑1785234
                                                                   |                                                                     
aaTAAGTGTCAGGCAAAAATGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCt                                                                   >  1:551624/1‑71 (MQ=255)
          aGGCAAAAATGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCa                                                          >  1:1173137/1‑70 (MQ=255)
          aGGCAAAAATGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCa                                                          >  1:3379237/1‑70 (MQ=255)
          aGGCAAAAATGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCa                                                          >  1:2279394/1‑70 (MQ=255)
            gCAAAAATGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATTGCTGCCGCTGCaaa                                                        <  1:1845036/71‑1 (MQ=255)
                aaaTGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAg                                                   >  1:1941211/1‑71 (MQ=255)
                  aTGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCaaaa                                                       >  1:2045179/1‑65 (MQ=255)
                  aTGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCaaaa                                                       >  1:2342070/1‑65 (MQ=255)
                  aTGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCt                                                 >  1:437320/1‑71 (MQ=255)
                          tGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCcac                                            <  1:3273478/68‑1 (MQ=255)
                           gCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCcac                                            <  1:5054265/67‑1 (MQ=255)
                           gCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAAc                                         >  1:5106757/1‑70 (MQ=255)
                           gCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAAc                                         >  1:1387152/1‑70 (MQ=255)
                                 aCGTTGTTTAAACAAACGCGTCCGGTAGATTTAATTCCTGCCGCTGCaaa                                                        <  1:3587360/50‑1 (MQ=39)
                                 aCGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTc                                               <  1:2674570/58‑1 (MQ=255)
                                 aCGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTacaac                                   >  1:754935/1‑70 (MQ=255)
                                 aCGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTacaac                                   >  1:4622074/1‑70 (MQ=255)
                                  cGTTGTTTAAACAAACTCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTacaaca                                  >  1:1517466/1‑70 (MQ=255)
                                   gttgttTAAACAAACGCGTCCGGTAGATTTAATTCCTGCCGCTGCAAAATCAGCTTCCACAACTacaaca                                  >  1:127006/1‑70 (MQ=255)
                                   gttgttTAAACAAACGCGTCCGGTAGATTTAATTCCTGCCGCTGCAAAATCAGCTTCCACAACTacaaca                                  >  1:3517031/1‑70 (MQ=255)
                                   gttgttTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACaaa                                <  1:2827606/71‑1 (MQ=255)
                                   gttgttTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACaaa                                >  1:698363/1‑71 (MQ=255)
                                     tgttTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCcac                                            <  1:1013923/57‑1 (MQ=255)
                                       tttAAACAAACGCGTCCGGTAGATTTAATTCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAAtt                              <  1:1641900/70‑1 (MQ=255)
                                              aaacGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCa                                                    >  1:1483161/1‑40 (MQ=255)
                                              aaacGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACaa                                          >  1:3357813/1‑50 (MQ=255)
                                                     tCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAAt                               <  1:844247/54‑1 (MQ=255)
                                                     tCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAAt                               <  1:3195919/54‑1 (MQ=255)
                                                     tCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAAt                               <  1:3607944/54‑1 (MQ=255)
                                                      ccGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACaaa               >  1:4708695/1‑69 (MQ=255)
                                                        ggTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTgg           <  1:690662/71‑1 (MQ=255)
                                                           aGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACaaa               >  1:85741/1‑64 (MQ=255)
                                                           aGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTc        >  1:333577/1‑71 (MQ=255)
                                                            gATTTAATTCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTc        >  1:846837/1‑71 (MQ=255)
                                                             aTTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAAtt                              <  1:4222609/47‑1 (MQ=255)
                                                             aTTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACaa                <  1:18060/61‑1 (MQ=255)
                                                             aTTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACaa                <  1:4069147/61‑1 (MQ=255)
                                                             aTTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACaa                <  1:1733217/61‑1 (MQ=255)
                                                               ttAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCg                          >  1:4236564/1‑49 (MQ=255)
                                                               ttAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTg                     <  1:5038461/54‑1 (MQ=255)
                                                               ttAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGc    >  1:586794/1‑71 (MQ=255)
                                                               ttAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGc    >  1:903505/1‑71 (MQ=255)
                                                               ttAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGc    >  1:1766332/1‑71 (MQ=255)
                                                               ttAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGc    >  1:3027286/1‑71 (MQ=255)
                                                               ttAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGAAAAACTGGCTCCCGc    >  1:4151131/1‑71 (MQ=255)
                                                               ttAACCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGc    >  1:1901279/1‑71 (MQ=255)
                                                                tAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGc    <  1:737615/70‑1 (MQ=255)
                                                                 aaTCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGt                            >  1:4313577/1‑45 (MQ=255)
                                                                 aaTCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGt                            >  1:1857382/1‑45 (MQ=255)
                                                                 aaTCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGt                            >  1:3095328/1‑45 (MQ=255)
                                                                 aaTCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCa  >  1:1740194/1‑71 (MQ=255)
                                                                 aaTCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCa  >  1:3894893/1‑71 (MQ=255)
                                                                 aaTCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCa  >  1:3628112/1‑71 (MQ=255)
                                                                 aaTCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCa  >  1:2708681/1‑71 (MQ=255)
                                                                   |                                                                     
AATAAGTGTCAGGCAAAAATGTTCTTTGCACCGACGTTGTTTAAACAAACGCGTCCGGTAGATTTAATCCTGCCGCTGCAAAATCAGCTTCCACAACTACAACAAATTGTCGGCGTGGACAAACTGGCTCCCGCCA  >  W3110S.gb/1785099‑1785234

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