Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,971,531 G→T 22.2% L330I (CTA→ATA)  ygeD ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,971,5310GT22.2% 64.5 / 6.6 27L330I (CTA→ATA) ygeDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (6/15);  new base T (2/4);  total (8/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCT  >  W3110S.gb/2971467‑2971598
                                                                |                                                                   
tCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACg                                                                <  1:3409901/70‑1 (MQ=255)
  gCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCa                                                              <  1:3574314/70‑1 (MQ=255)
  gCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCa                                                              <  1:627546/70‑1 (MQ=255)
  gCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCa                                                              <  1:1654915/70‑1 (MQ=255)
     aaGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGa                                                          <  1:373285/71‑1 (MQ=255)
          tttGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCAt                                                             >  1:1654258/1‑63 (MQ=255)
                cTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTATCAACg                                                                <  1:437875/54‑1 (MQ=255)
                  gCAATCGCGTTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCaa                                                                  >  1:1218043/1‑50 (MQ=255)
                  gCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCaa                                                                  >  1:2052929/1‑50 (MQ=255)
                   cAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTATCAACGCATTGAGCGGAACGACaaaa                                            <  1:3225231/71‑1 (MQ=255)
                   cAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACaaa                                             <  1:1032154/70‑1 (MQ=255)
                   cAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACaaa                                             <  1:2253973/70‑1 (MQ=255)
                   cAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACaaa                                             <  1:1556481/70‑1 (MQ=255)
                    aaTCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTATCAACGCATTGAGCGGAACGACaaaa                                            <  1:1295337/70‑1 (MQ=255)
                      tCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGc                                                               <  1:1369050/49‑1 (MQ=255)
                        gCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTATCAACGCATTGAGCGGAACGACaaaa                                            >  1:2943817/1‑66 (MQ=255)
                             ccccgccccgACGCTTTTTTTACCCCGCTCCTGTAGCAACGCAtt                                                            >  1:62771/1‑45 (MQ=255)
                              cccgccccgACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCAt                                 <  1:1488914/71‑1 (MQ=255)
                                       aCGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGcc                           <  1:3089846/68‑1 (MQ=255)
                                        cGCTTTTTTTACCCCGCTCCTGTATCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATc                       <  1:2661652/71‑1 (MQ=255)
                                         gCTTTTTTTACCCCGCTCCTGTATCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAAt                        >  1:2079699/1‑69 (MQ=255)
                                                  aCCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATcaac              <  1:2153621/70‑1 (MQ=255)
                                                   ccccGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATcaac              <  1:1130279/69‑1 (MQ=255)
                                                     ccGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATCAACAAgg          <  1:1058094/71‑1 (MQ=255)
                                                       gCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATCAACAAGGCa        <  1:1256204/71‑1 (MQ=255)
                                                             gTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCt  >  1:3563030/1‑71 (MQ=255)
                                                             gTAGCAACGCATTGAGCGGAAAGACAAAAAAGCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCt  >  1:879904/1‑71 (MQ=255)
                                                                |                                                                   
TCGCCAAGGTTTTGTACTGCAATCGCATTCCCCGCCCCGACGCTTTTTTTACCCCGCTCCTGTAGCAACGCATTGAGCGGAACGACAAAAAAGCCCCCCATCACGCCAATCAGCATCAACAAGGCATAGGCT  >  W3110S.gb/2971467‑2971598

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: