Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 861,111 G→A 46.9% 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,1110GA46.9% 2.4 / 41.4 32L640L (CTG→TTG) ybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base G (8/9);  new base A (8/7);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

TGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACT  >  W3110S.gb/861043‑861179
                                                                    |                                                                    
tgaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAAtt                                                                    >  1:482732/1‑71 (MQ=255)
   cGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgt                                                                 >  1:1646436/1‑71 (MQ=255)
              aaTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTc                                                      >  1:1086887/1‑71 (MQ=255)
               aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                     >  1:1128239/1‑71 (MQ=255)
               aTAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                     <  1:794715/71‑1 (MQ=255)
                   ttGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGAt                                                  >  1:379092/1‑70 (MQ=255)
                   ttGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGAt                                                  >  1:1541288/1‑70 (MQ=255)
                        gccAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAgg                                                        <  1:241838/59‑1 (MQ=255)
                           aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAg                                         >  1:1735579/1‑71 (MQ=255)
                           aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGa                                                              >  1:1153114/1‑50 (MQ=255)
                           aaCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAg                                         >  1:576372/1‑71 (MQ=255)
                           aaCCGGACCACGACCGTAGCGCGGATTATGCTACTGTTTCAATTCGTCGa                                                              >  1:1700995/1‑50 (MQ=38)
                            aCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGa                                                              >  1:1738525/1‑49 (MQ=255)
                            aCCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAAtt                                                                    <  1:216067/43‑1 (MQ=255)
                                gACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcag                                     >  1:1185233/1‑70 (MQ=255)
                                    aCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg                                >  1:1156623/1‑71 (MQ=255)
                                       aCCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTAt                             <  1:559479/71‑1 (MQ=255)
                                         cGTAGCGCGGATTATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCg                           <  1:943571/71‑1 (MQ=255)
                                          gTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg                                <  1:612948/65‑1 (MQ=255)
                                          gTAGCCCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCg                                <  1:1327454/65‑1 (MQ=255)
                                           tAGCGCGGATTATGGTACTGTTTCAGTTCGTCGata                                                            <  1:121567/36‑1 (MQ=255)
                                           tAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                          >  1:74334/1‑70 (MQ=255)
                                                    ttATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATa                                                 <  1:798696/38‑1 (MQ=255)
                                                     tATGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcg               >  1:1196578/1‑71 (MQ=255)
                                                      aTGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcg               <  1:686414/70‑1 (MQ=255)
                                                      aTGGTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGa                          >  1:43238/1‑59 (MQ=255)
                                                         gTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCcagcag                                  <  1:893660/48‑1 (MQ=255)
                                                         gTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgt              <  1:1242903/68‑1 (MQ=255)
                                                         gTACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgt              <  1:734737/68‑1 (MQ=255)
                                                          tACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                            <  1:20794/37‑1 (MQ=255)
                                                          tACTGTTTCAATTCGTCGATATAGGTCTGATAAGCGCGAGCca                                      <  1:1672923/43‑1 (MQ=255)
                                                                   aGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACt  <  1:358786/70‑1 (MQ=255)
                                                                    |                                                                    
TGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACT  >  W3110S.gb/861043‑861179

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: