Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,958,682 C→T 12.2% A69V (GCC→GTC)  yhjJ → predicted zinc‑dependent peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,958,6820CT12.2% 115.7 / 6.1 41A69V (GCC→GTC) yhjJpredicted zinc‑dependent peptidase
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base T (2/3);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACCACCCCCCAGCGTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTGCGCTAACGCAAAGCGGTGGCC  >  W3110S.gb/3958615‑3958750
                                                                    |                                                                    
cTGACCACCCCCCAGCGTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGa                                                                   <  1:2140372/71‑1 (MQ=255)
 tGACCACCCCCCAGCGTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGa                                                                   <  1:2405611/70‑1 (MQ=255)
   accaccCCCCAGCGTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGaa                                                                  <  1:2236632/69‑1 (MQ=255)
        cccccAGCGTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTacac                                                           <  1:1263939/71‑1 (MQ=255)
              gCGTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACacaaca                                                       >  1:2192850/1‑69 (MQ=255)
              gCGTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACacaaca                                                       >  1:2263816/1‑69 (MQ=255)
                gTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGc                                                   >  1:1635701/1‑71 (MQ=255)
                gTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGc                                                   >  1:140725/1‑71 (MQ=255)
                 tCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCg                                                  <  1:1470319/71‑1 (MQ=255)
                 tCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCg                                                  <  1:386049/71‑1 (MQ=255)
                       cGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTa                                              >  1:2388364/1‑69 (MQ=255)
                           cGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGTCGAAAGTACACAACAGAGCGGTTACAGTCa                                        >  1:2126181/1‑71 (MQ=255)
                             tgtTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACg                                      <  1:1611095/71‑1 (MQ=255)
                               ttGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGa                                                                   <  1:1157347/40‑1 (MQ=255)
                               ttGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGc                                     <  1:1392364/70‑1 (MQ=255)
                               ttGAAATTCGCCTGCTCGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGc                                     <  1:1348158/70‑1 (MQ=255)
                                 gAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCAt                                  >  1:2485462/1‑71 (MQ=255)
                                 gAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCAt                                  >  1:200406/1‑71 (MQ=255)
                                  aaaTTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAg                                           >  1:1661711/1‑61 (MQ=255)
                                  aaaTTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAg                                           >  1:1658513/1‑61 (MQ=255)
                                       cGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACa                                            >  1:1904859/1‑55 (MQ=255)
                                         ccTGCTGGTTAAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCg                            <  1:988793/70‑1 (MQ=255)
                                         ccTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCg                            <  1:2209205/69‑1 (MQ=255)
                                         ccTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCg                            <  1:1894278/69‑1 (MQ=255)
                                         ccTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCCGTAACAGTCACGCCATCCCTCGTa                          >  1:1721120/1‑71 (MQ=255)
                                               ggTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTgcg                     >  1:465394/1‑70 (MQ=255)
                                                gTTAATACCGGTTCGCTCGTCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATc                                 <  1:2241334/57‑1 (MQ=255)
                                                gTTAATACCGGTTCGCTCGTCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATc                                 <  1:2761439/57‑1 (MQ=255)
                                                 ttAATACCGGTTCCCTCGCCGAAAGTACACAACaga                                                     <  1:1220298/36‑1 (MQ=255)
                                                    aaTACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTGCGCTAAc                >  1:941834/1‑71 (MQ=255)
                                                     aTACCGGTTCGCTCGTCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTGCGCTAACg               >  1:294849/1‑71 (MQ=255)
                                                      tACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGtt                                               <  1:2224328/38‑1 (MQ=255)
                                                      tACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGtt                                               <  1:2227399/38‑1 (MQ=255)
                                                      tACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGc                                     <  1:887385/48‑1 (MQ=255)
                                                       aCCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTGCGCTAACg               >  1:557493/1‑69 (MQ=255)
                                                        ccGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTAt                         >  1:2505073/1‑58 (MQ=255)
                                                        ccGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTAt                         >  1:2305572/1‑58 (MQ=255)
                                                          ggTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATCGCGCTAAc                >  1:566787/1‑65 (MQ=255)
                                                           gTTCGCTCGTCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATc                                 <  1:2587256/47‑1 (MQ=255)
                                                            ttCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTGCGCTAACGCAAAGc         <  1:2721912/70‑1 (MQ=255)
                                                            ttCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTGCGCTAACGCAAAGc         <  1:1246665/70‑1 (MQ=255)
                                                             tcgctcgcCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTa                          >  1:1713211/1‑52 (MQ=255)
                                                             tcgctcgcCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTa                          >  1:176130/1‑52 (MQ=255)
                                                              cgctcgtCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATcc                                >  1:1926570/1‑45 (MQ=255)
                                                              cgctcgcCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTAtt                        >  1:1712396/1‑53 (MQ=255)
                                                                ctcgtCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTGCGCTAACGCAAAGCGGTgg    <  1:678802/71‑1 (MQ=255)
                                                                  cgcCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTGCGCTAACGCAAAGCGGTGGcc  >  1:2400685/1‑71 (MQ=255)
                                                                    |                                                                    
CTGACCACCCCCCAGCGTCCCAGCGATCGTGTTGAAATTCGCCTGCTGGTTAATACCGGTTCGCTCGCCGAAAGTACACAACAGAGCGGTTACAGTCACGCCATCCCTCGTATTGCGCTAACGCAAAGCGGTGGCC  >  W3110S.gb/3958615‑3958750

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: