Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 861,111 G→A 35.1% L640L (CTG→TTG)  ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb861,1110GA35.1% 63.4 / 54.7 57L640L (CTG→TTG) ybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (13/24);  new base A (14/6);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

TGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTA  >  W3110S.gb/861043‑861177
                                                                    |                                                                  
tgaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGtt                                                                  <  1:2696483/71‑1 (MQ=255)
tgaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGtt                                                                  <  1:1432913/71‑1 (MQ=255)
  aCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcg                                                                >  1:4147928/1‑71 (MQ=255)
    gTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtc                                                              >  1:3523417/1‑71 (MQ=255)
     tACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtcg                                                             <  1:1282356/71‑1 (MQ=255)
      aCCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtcg                                                             <  1:2441284/70‑1 (MQ=255)
       cccGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                                            <  1:4925954/70‑1 (MQ=255)
         cgcgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGatat                                                         <  1:2144494/71‑1 (MQ=255)
           cgTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAg                                                       <  1:1374212/71‑1 (MQ=255)
             tAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGt                                                                   <  1:3599830/57‑1 (MQ=255)
             tAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                                            <  1:3861435/64‑1 (MQ=255)
             tAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGt                                                     >  1:2997631/1‑71 (MQ=255)
               aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                   >  1:4951682/1‑71 (MQ=255)
               aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCt                                                   >  1:3885783/1‑71 (MQ=255)
               aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCt                                                   >  1:2557563/1‑71 (MQ=255)
                tAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAg                                                       <  1:1831219/66‑1 (MQ=255)
                tAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAg                                                       <  1:4479774/66‑1 (MQ=255)
                  gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGt                                                     >  1:4865213/1‑66 (MQ=255)
                  gTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGa                                                 >  1:2580554/1‑70 (MQ=255)
                    tGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATa                                               <  1:865831/70‑1 (MQ=255)
                        gccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAgcgc                                          <  1:199855/71‑1 (MQ=255)
                           aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGa                                        >  1:4862934/1‑70 (MQ=255)
                           aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAg                                       >  1:1607168/1‑71 (MQ=255)
                             ccGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                   <  1:4380620/57‑1 (MQ=255)
                              cGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGAtata                                                        >  1:2966023/1‑51 (MQ=255)
                              cGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGcc                                     <  1:5145339/70‑1 (MQ=255)
                              cGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATaa                                              >  1:1348817/1‑61 (MQ=255)
                              cGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATaa                                              >  1:1051583/1‑61 (MQ=255)
                               ggACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcag                                   >  1:4038961/1‑71 (MQ=255)
                                  ccACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagcag                                <  1:2592458/71‑1 (MQ=255)
                                  ccACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCcagcag                                >  1:2838576/1‑71 (MQ=255)
                                   cACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagcag                                <  1:875159/70‑1 (MQ=255)
                                    aCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                          <  1:306357/59‑1 (MQ=255)
                                    aCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg                              >  1:2031801/1‑71 (MQ=255)
                                     cGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGAtata                                                        >  1:2932438/1‑44 (MQ=255)
                                     cGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg                              >  1:1969905/1‑70 (MQ=255)
                                     cGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGt                             <  1:2613734/71‑1 (MQ=255)
                                         cGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAg                                                       >  1:1397353/1‑41 (MQ=255)
                                          gTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCATCAGCGTATCGa                        <  1:3127468/71‑1 (MQ=255)
                                          gTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                        <  1:2414963/71‑1 (MQ=255)
                                          gTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                        <  1:4857184/71‑1 (MQ=255)
                                              cgcgGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGt                    <  1:3470682/71‑1 (MQ=255)
                                              cgcgGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGt                    <  1:2809458/71‑1 (MQ=255)
                                                  gATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                          <  1:2596604/45‑1 (MQ=255)
                                                  gATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                          <  1:2280321/45‑1 (MQ=255)
                                                  gATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                          <  1:4970165/45‑1 (MQ=255)
                                                   aTTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                   <  1:214345/35‑1 (MQ=255)
                                                   aTTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                   <  1:1811852/35‑1 (MQ=255)
                                                     tATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcg             >  1:3356060/1‑71 (MQ=255)
                                                       tGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgtc           <  1:645603/71‑1 (MQ=255)
                                                         gTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATcgtcgt         >  1:4561837/1‑71 (MQ=255)
                                                         gTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATcgtcgt         >  1:3678870/1‑71 (MQ=255)
                                                         gTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGAGAGCCAGCAGCGTATCGACAGTATCATcgtcgt         >  1:1457864/1‑71 (MQ=255)
                                                           aCTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTg       >  1:4255439/1‑71 (MQ=255)
                                                              gTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCg    >  1:4331342/1‑71 (MQ=255)
                                                              gTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCg    >  1:2636280/1‑71 (MQ=255)
                                                                ttCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTa  >  1:2102297/1‑71 (MQ=255)
                                                                    |                                                                  
TGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTA  >  W3110S.gb/861043‑861177

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: