Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,331,819 A→G 17.1% F258L (TTT→CTT)  yhbE ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,331,8190AG17.1% 75.1 / 11.1 35F258L (TTT→CTT) yhbEconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (17/12);  new base G (2/4);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

TTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGC  >  W3110S.gb/3331752‑3331888
                                                                   |                                                                     
tttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAGCAg                                                                    >  1:1104592/1‑71 (MQ=255)
 ttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGt                                                                   >  1:465209/1‑71 (MQ=255)
  tAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTaaa                                                                       >  1:2626522/1‑66 (MQ=255)
  tAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTaaa                                                                       >  1:3666767/1‑66 (MQ=255)
  tAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTg                                                                  >  1:3259169/1‑71 (MQ=255)
   aaCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTg                                                                  <  1:59296/70‑1 (MQ=255)
    aCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTggg                                                                <  1:3850579/71‑1 (MQ=255)
    aCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTgg                                                                 <  1:1401576/70‑1 (MQ=255)
       tCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGt                                                              >  1:3244769/1‑70 (MQ=255)
       tCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGt                                                              >  1:2010292/1‑70 (MQ=255)
       tCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGt                                                              >  1:2236098/1‑70 (MQ=255)
         ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTaaa                                                                       >  1:2221952/1‑59 (MQ=255)
         ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAGCAGTGGGGTGAg                                                           <  1:555042/71‑1 (MQ=255)
         ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAGCAGTGGGGTGAg                                                           <  1:1050321/71‑1 (MQ=255)
               gCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAGCAGTGGGGTGAGCGtgat                                                     >  1:1446546/1‑71 (MQ=255)
               gCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgat                                                     >  1:3757518/1‑71 (MQ=255)
                    gaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGt                                                                   >  1:3880378/1‑52 (MQ=255)
                    gaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTggg                                                                >  1:3490137/1‑55 (MQ=255)
                    gaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACACCGTAAACAGTggg                                                                >  1:3280445/1‑55 (MQ=255)
                         cGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTggg                                                                >  1:2194660/1‑50 (MQ=255)
                                    agTGATAAATGATCTGTAAAAAACAGCGTAAGCAGTGGGGTGAg                                                           <  1:1084743/44‑1 (MQ=37)
                                     gTGATAAAAGATCTGAAAAAAACAGCGTAAGCAGTggg                                                                <  1:1113142/38‑1 (MQ=255)
                                     gTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTgg                                                                 <  1:2525784/37‑1 (MQ=255)
                                     gTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTgg                                                                 <  1:1722539/37‑1 (MQ=255)
                                       gATAAAAGATCTGAAAAAGACAGCGTAAACAGTGGGGTGa                                                            >  1:1862961/1‑40 (MQ=255)
                                         tAAAAGATCTGAAAAAAACAGCATAAACAGTGGGGTGAGCGtg                                                       <  1:380596/43‑1 (MQ=255)
                                          aaaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctgc                            >  1:434549/1‑69 (MQ=255)
                                           aaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                    <  1:1179275/44‑1 (MQ=255)
                                           aaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgat                                                     >  1:3107193/1‑43 (MQ=255)
                                              gATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCa                          >  1:1214292/1‑67 (MQ=255)
                                              gATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGc                        >  1:3446368/1‑69 (MQ=255)
                                                    aaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCa                  <  1:1063351/69‑1 (MQ=255)
                                                    aaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                 <  1:782002/70‑1 (MQ=255)
                                                    aaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATc                <  1:479749/71‑1 (MQ=255)
                                                    aaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATc                <  1:2658487/71‑1 (MQ=255)
                                                          aCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCt                                            >  1:2993654/1‑37 (MQ=255)
                                                          aCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCt                                            >  1:2746445/1‑37 (MQ=255)
                                                             gCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGcc       <  1:3463030/71‑1 (MQ=255)
                                                                  aaCAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGgc  >  1:3334035/1‑71 (MQ=255)
                                                                   |                                                                     
TTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGC  >  W3110S.gb/3331752‑3331888

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