Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 516,807 C→T 36.6% E233K (GAA→AAA)  ybbN ← predicted thioredoxin domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb516,8070CT36.6% 37.2 / 39.4 41E233K (GAA→AAA) ybbNpredicted thioredoxin domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (18/8);  new base T (11/4);  total (29/12)
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.94e-01

CGCGTCTGACCGTCTGCGGCGGTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCGGCATCTTCTGGA  >  W3110S.gb/516740‑516868
                                                                   |                                                             
cgcgTCTGACCGTCTGCGGCGGTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCAt                                                             <  1:160597/70‑1 (MQ=255)
  cgTCTGACCGTCTGCGGCGGTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTTATTgc                                                          <  1:1314109/71‑1 (MQ=255)
      tGACCGTCTGCGGCGGTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGcc                                                       >  1:2484972/1‑70 (MQ=255)
      tGACCGTCTGCGGCGGTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGcc                                                       >  1:744631/1‑70 (MQ=255)
              tGCGGCGGTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGAt                                              >  1:672406/1‑71 (MQ=255)
              tGCGGCGGTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGAt                                              >  1:3610763/1‑71 (MQ=255)
                  gcggTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAg                                          >  1:2955242/1‑71 (MQ=255)
                     gTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTg                                       >  1:3470256/1‑71 (MQ=255)
                     gTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTg                                       >  1:1884469/1‑71 (MQ=255)
                       gagaTCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGa                                               >  1:317691/1‑61 (MQ=255)
                       gagaTCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGa                                               >  1:312928/1‑61 (MQ=255)
                         gaTCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAg                                    >  1:3901314/1‑70 (MQ=255)
                               tACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTTATTGCGCCCAACCTGATGCa                                           >  1:3827722/1‑57 (MQ=255)
                               tACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAGCGcc                                <  1:3146607/68‑1 (MQ=255)
                               tACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATACAGTTGCAGCGcc                                <  1:2056667/68‑1 (MQ=255)
                                aCGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGtt                            >  1:3136956/1‑71 (MQ=255)
                                aCGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGtt                            >  1:2877455/1‑71 (MQ=255)
                                   cAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCg                         <  1:254622/71‑1 (MQ=255)
                                   cAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGAAGCGCCAGTTGCg                         <  1:1186458/71‑1 (MQ=255)
                                      aTGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGt                                         <  1:1742545/52‑1 (MQ=255)
                                      aTGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGt                                         <  1:3411565/52‑1 (MQ=255)
                                       tGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGa                                               >  1:1291824/1‑45 (MQ=255)
                                       tGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCa                                     >  1:1331413/1‑55 (MQ=255)
                                        gCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCa                                     >  1:1687816/1‑54 (MQ=255)
                                        gCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGc                     >  1:30336/1‑70 (MQ=255)
                                         cccGAACAGCAACTCCAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGcc                    <  1:1391112/70‑1 (MQ=255)
                                           cGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGt                 >  1:139949/1‑71 (MQ=255)
                                               cAGCAACTCCAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCgg             >  1:2068884/1‑71 (MQ=255)
                                               cAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCgg             >  1:382617/1‑71 (MQ=255)
                                               cAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCg              >  1:200957/1‑70 (MQ=255)
                                                aGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAGCGcc                                <  1:184282/51‑1 (MQ=255)
                                                   aaCTCCAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCGGCAtc         >  1:403502/1‑71 (MQ=255)
                                                     cTCCAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCGGCAtc         <  1:2153441/69‑1 (MQ=255)
                                                      tCCAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCGGCAtctt       >  1:3321555/1‑70 (MQ=255)
                                                      tCCAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCGGCAtctt       >  1:3506513/1‑70 (MQ=255)
                                                       ccAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTg                                       >  1:1604249/1‑37 (MQ=255)
                                                       ccAGCGCCTCTTTATTGCGCCCAACCTGATGCAGTTg                                       >  1:1941603/1‑37 (MQ=255)
                                                            gcCTCTTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCGGCATCTTCTGGa  >  1:3728830/1‑69 (MQ=255)
                                                               tctTTATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCGGCATCTTCTg    <  1:1228857/64‑1 (MQ=255)
                                                                  tCATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAg                  >  1:480422/1‑47 (MQ=255)
                                                                  tCATTGCGCCCAACCTGATGAAGTTGCAGCGCCAGTTGCGTCGCCAg                  >  1:3377125/1‑47 (MQ=255)
                                                                   |                                                             
CGCGTCTGACCGTCTGCGGCGGTGAGATCTTTACGCAGATGCCCGAACAGCAACTCCAGCGCCTCTTCATTGCGCCCAACCTGATGCAGTTGCAGCGCCAGTTGCGTCGCCAGTGCGGCATCTTCTGGA  >  W3110S.gb/516740‑516868

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: