Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,810,661 A→G 73.0% F96S (TTT→TCT)  yniB ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,810,6610AG73.0% 60.8 / 23.3 37F96S (TTT→TCT) yniBpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (5/5);  new base G (17/10);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGT  >  W3110S.gb/1810597‑1810727
                                                                |                                                                  
ttCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAAt                                                                >  1:2217743/1‑69 (MQ=255)
ttCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAAt                                                                >  1:481069/1‑69 (MQ=255)
ttCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAAt                                                                >  1:1303221/1‑69 (MQ=255)
  cccGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAg                                                             <  1:137259/70‑1 (MQ=255)
      cAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATg                                                         <  1:139145/70‑1 (MQ=255)
        agaagCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGaaa                                                      <  1:1582529/71‑1 (MQ=255)
          aagCGCGCCTGGCGACTCATCCGGGCGCCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAAt                                                    >  1:1010991/1‑71 (MQ=255)
          aagCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAAt                                                    >  1:143532/1‑71 (MQ=255)
                 ccTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAg                                             <  1:2316021/71‑1 (MQ=255)
                    ggCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGaa                                           >  1:2091787/1‑70 (MQ=255)
                            aTCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGtt                                  <  1:1589506/71‑1 (MQ=255)
                            aTCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGtt                                  <  1:2028409/71‑1 (MQ=255)
                             tCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGtt                                  <  1:377967/70‑1 (MQ=255)
                               cGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGaaa                                                      >  1:1923131/1‑48 (MQ=255)
                               cGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGaaa                                                      >  1:2597013/1‑48 (MQ=255)
                               cGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGaaa                                                      >  1:1062002/1‑48 (MQ=255)
                               cGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCa                                <  1:2403415/70‑1 (MQ=255)
                               cGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCa                                >  1:1884009/1‑70 (MQ=255)
                               cGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCa                                >  1:1701657/1‑70 (MQ=255)
                               cGGGCACCGGAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCAc                               >  1:1853138/1‑71 (MQ=255)
                               cGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCAc                               <  1:2294957/71‑1 (MQ=255)
                                    aCCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGccac                          <  1:2262038/71‑1 (MQ=255)
                                        gAATCCTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCAt                       >  1:2413427/1‑70 (MQ=255)
                                         aaTCCTGCAACGCCAGTCCGACAGAAATCAGGATGaa                                                       <  1:1375189/37‑1 (MQ=255)
                                             cTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATa                                       <  1:2082163/49‑1 (MQ=255)
                                             cTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCa                                >  1:1515307/1‑56 (MQ=255)
                                             cTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGccac                          >  1:2040394/1‑62 (MQ=255)
                                             cTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTaaaa                  >  1:2037924/1‑70 (MQ=255)
                                             cTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTaaaa                  >  1:2717081/1‑70 (MQ=255)
                                             cTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAAt                 >  1:6943/1‑71 (MQ=255)
                                             cTGCAACGCCAGTCCGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAAt                 >  1:233852/1‑71 (MQ=255)
                                             cTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAAt                 >  1:1931763/1‑71 (MQ=255)
                                             cTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCACAACATAACGTTCACGCCACCATTAAAAt                 >  1:114330/1‑71 (MQ=255)
                                                     ccAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAAtt                <  1:1367543/64‑1 (MQ=255)
                                                          ccGACAGAAATCAGGTTGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCggg      >  1:1669447/1‑69 (MQ=255)
                                                           cGACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGa    <  1:2220439/70‑1 (MQ=255)
                                                            gACAGAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGt  <  1:1174620/71‑1 (MQ=255)
                                                                |                                                                  
TTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGT  >  W3110S.gb/1810597‑1810727

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: