Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 281,512 C→T 35.1% S11F (TCC→TTC)  yagE → predicted lyase/synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb281,5120CT35.1% 62.3 / 54.9 57S11F (TCC→TTC) yagEpredicted lyase/synthase
Reads supporting (aligned to +/- strand):  ref base C (20/17);  new base T (9/11);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAA  >  W3110S.gb/281444‑281578
                                                                    |                                                                  
ttAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcg                                                                  <  1:1113284/71‑1 (MQ=255)
 taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTTcgc                                                                 >  1:2613515/1‑71 (MQ=255)
 taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTTcgc                                                                 >  1:1564438/1‑71 (MQ=255)
 taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcgc                                                                 >  1:4561368/1‑71 (MQ=255)
 taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcg                                                                  <  1:2130411/70‑1 (MQ=255)
   aataaaCGAACGCTGTTCTATAATGTACAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt                                                               >  1:2937744/1‑71 (MQ=255)
       aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTTCGCgttgtt                                                           <  1:440020/71‑1 (MQ=255)
       aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTTCGCgttgtt                                                           <  1:4134836/71‑1 (MQ=255)
       aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgttgtt                                                           <  1:1321301/71‑1 (MQ=255)
               cTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGaa                                                   >  1:2537110/1‑71 (MQ=255)
                  ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTTCGCGTTGTTCACGGGAATc                                                 <  1:3936171/70‑1 (MQ=255)
                  ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTTCGCGTTGTTCACGGGAATc                                                 <  1:4599203/70‑1 (MQ=255)
                  ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTTCGCGTTGTTCACGGGAATc                                                 <  1:1060852/70‑1 (MQ=255)
                  ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAAt                                                  <  1:3051892/69‑1 (MQ=255)
                  ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAAt                                                  <  1:2500371/69‑1 (MQ=255)
                   tCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTTCGCGTTGTTCACGGGAATCa                                                >  1:1722652/1‑70 (MQ=255)
                   tCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCAt                                               >  1:4330783/1‑71 (MQ=255)
                       taATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTc                                                          >  1:2684504/1‑56 (MQ=255)
                         aTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCa                                                <  1:3098020/64‑1 (MQ=255)
                         aTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCa                                                <  1:3939169/64‑1 (MQ=255)
                          tGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCt                                        <  1:506249/71‑1 (MQ=255)
                             aGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGtc                                     >  1:156529/1‑71 (MQ=255)
                              gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCa                                                         >  1:4959120/1‑50 (MQ=255)
                              gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCa                                                >  1:3861728/1‑59 (MQ=255)
                                 cAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCa                                                         >  1:3855969/1‑47 (MQ=255)
                                         ttCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCAttttt                         <  1:4498779/71‑1 (MQ=255)
                                         ttCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCAttttt                         <  1:478625/71‑1 (MQ=255)
                                          tCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACg                                                       >  1:1803/1‑40 (MQ=255)
                                             gCAAGGAGATCTCATGCCGCAGTTCGCGTTGTTCAc                                                        >  1:633848/1‑36 (MQ=255)
                                             gCAAGGAGATCTCATGCCGCAGTTCGCGTTGTTCAc                                                        >  1:323982/1‑36 (MQ=255)
                                             gCAAGGAGATCTCATGCCGCAGTTCGCGTTGTTCAc                                                        >  1:185258/1‑36 (MQ=255)
                                              cAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTccacca                              >  1:138058/1‑61 (MQ=255)
                                              cAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTccacca                              >  1:2861817/1‑61 (MQ=255)
                                              cAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTccacca                              >  1:2002334/1‑61 (MQ=255)
                                                aGGAGATCTCATGCCGCAGTTCGCGTTGTTCACGGGaa                                                   >  1:322524/1‑38 (MQ=255)
                                                aGGAGATCTCATGCCGCAGTTCGCGTTGTTCACGGGAATCATTCCCCCt                                        <  1:1430523/49‑1 (MQ=255)
                                                aGGAGATCTCATGCCGCAGTTCGCGTTGTTCACGGGAATCATTCCCCCTGtctc                                   >  1:1450580/1‑54 (MQ=255)
                                                 ggAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTccacca                              <  1:2376831/58‑1 (MQ=255)
                                                 ggAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                   <  1:3624436/69‑1 (MQ=255)
                                                 ggAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                   <  1:4778421/69‑1 (MQ=255)
                                                  gagaTCTCATGCCGCAGTTCGCGTTGTTCACGGGAATCATTccc                                           >  1:606002/1‑44 (MQ=255)
                                                   agaTCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACg               >  1:4815429/1‑71 (MQ=255)
                                                     aTCTCATGCCGCAGTTCGCGTTGTTCACGGGAATCATTCCCCCt                                        <  1:607206/44‑1 (MQ=255)
                                                     aTCTCATGCCGCAGTTCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCAttttt                         <  1:3432388/59‑1 (MQ=255)
                                                     aTCTCATGCCGCAGTTCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCAttttt                         <  1:4721943/59‑1 (MQ=255)
                                                     aTCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTccacca                              <  1:3519482/54‑1 (MQ=255)
                                                       ctcATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCAtttt                          >  1:2089980/1‑56 (MQ=255)
                                                       ctcATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCAttt                           >  1:3672794/1‑55 (MQ=255)
                                                       ctcATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGccg           >  1:3607966/1‑70 (MQ=255)
                                                        tcATGCCGCAGTCCGCGTTGTTCACGGGAATCAtt                                              >  1:2145923/1‑35 (MQ=255)
                                                        tcATGCCGCAGTCCGCGTTGTTCACGGGAATCAtt                                              >  1:3441225/1‑35 (MQ=255)
                                                         cATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAg          >  1:3964814/1‑70 (MQ=255)
                                                          aTGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGcc            <  1:240697/67‑1 (MQ=255)
                                                           tGCCGCAGTTCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTc       <  1:2319499/71‑1 (MQ=255)
                                                           tGCCGCAGTTCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTc       <  1:403982/71‑1 (MQ=255)
                                                             ccGCAGTCCGCGTTGTTCAGGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGa     <  1:1755981/71‑1 (MQ=255)
                                                                cAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATaa  <  1:2723815/71‑1 (MQ=255)
                                                                    |                                                                  
TTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAA  >  W3110S.gb/281444‑281578

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: