Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,393,837 C→T 59.1% R45* (CGA→TGA)  ycjZ → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,393,8370CT59.1% 8.5 / 22.4 22R45* (CGA→TGA) ycjZpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (2/7);  new base T (8/5);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 9.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTGCAGCAGCCCGCCTGAGCATGGCGCAGTCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAACGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCAACTGAAGCGGGCGAGCATCTTTT  >  W3110S.gb/1393770‑1393907
                                                                   |                                                                      
tCGTGCAGCAGCCCGCCTGAGCATGGCGCAGTCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAATGAt                                                                     >  1:753989/1‑71 (MQ=255)
       gcagcCCGCCTGAGCATGGCGCAGTCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAATGATTGGGAtt                                                              <  1:690595/71‑1 (MQ=255)
        cagcCCGCCTGAGCATGGCGCAGTCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAATGATTGGGATTg                                                             >  1:469990/1‑71 (MQ=255)
           cccGCCTGAGCATGGCGCAGTCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAACGATTGGGATTGCgg                                                          <  1:80991/71‑1 (MQ=255)
            ccGCCTGAGCATGGCGCAGTCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAATGATTGGGATTGCGGc                                                         >  1:1178924/1‑71 (MQ=255)
             cGCCTGAGCATGGCGCAGTCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAATGATTGGGATTGCGGCt                                                        <  1:79963/71‑1 (MQ=255)
             cGCCTGAGCATGGCGCAGTCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAACGATTGGGATTGCGGCt                                                        <  1:744163/71‑1 (MQ=255)
                               tCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAATGATTGGGATTGCGGCTTCTGACGCGAACCAcgcg                                      >  1:376590/1‑71 (MQ=255)
                               tCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAATGATTGGGATTGCGGCTTCTGACGCGAACCAcgcg                                      >  1:1426595/1‑71 (MQ=255)
                               tCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAACGATTGGGAtt                                                              <  1:1421885/47‑1 (MQ=255)
                                 aGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAATGAt                                                                     <  1:1057032/38‑1 (MQ=255)
                                 aGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAACGATTGGGATTGCGGCTTCTGACGCGaa                                             >  1:727789/1‑62 (MQ=255)
                                                 gTGCGTCGTATAGAAGAATGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCa                                    >  1:346817/1‑55 (MQ=255)
                                                 gTGCGTCGTATAGAAGAATGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCa                                    >  1:918942/1‑55 (MQ=255)
                                                      tcgtATAGAAGAATGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCa                         <  1:1599040/61‑1 (MQ=255)
                                                         tataGAAGAACGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCgt                                <  1:462378/51‑1 (MQ=255)
                                                             gaagaaTGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCAACTGAAGCGGGCGAGc         <  1:1521103/70‑1 (MQ=255)
                                                             gaagaaTGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCAACTGAAGCGGGCGAGCa        >  1:1333656/1‑71 (MQ=255)
                                                               agaaCGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCAACTGAAGCGGGCGAGCATc      >  1:363086/1‑71 (MQ=255)
                                                                gaaCGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCAACTGAAGCGGGCGATCATCt     <  1:1305233/71‑1 (MQ=255)
                                                                 aaCGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCAACTGAAGCGGGCGAGCATCtt    <  1:379828/71‑1 (MQ=255)
                                                                  aCGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCAACTGAAGCGGGCGAGCATCtt    <  1:592429/70‑1 (MQ=255)
                                                                   cGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCAACTGAAGCGGGCGAGCATCtttt  <  1:1745577/71‑1 (MQ=255)
                                                                   |                                                                      
TCGTGCAGCAGCCCGCCTGAGCATGGCGCAGTCAGCTTTAAGCCAGATAGTGCGTCGTATAGAAGAACGATTGGGATTGCGGCTTCTGACGCGAACCACGCGCAGCGTTGTTCCAACTGAAGCGGGCGAGCATCTTTT  >  W3110S.gb/1393770‑1393907

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: