Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,835,870 C→T 22.9% A14V (GCT→GTT)  ydjY → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,835,8700CT22.9% 77.2 / 19.6 40A14V (GCT→GTT) ydjYhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/10);  new base T (2/7);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 2.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CCAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGC  >  W3110S.gb/1835808‑1835930
                                                              |                                                            
ccAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCAct                                                        >  1:563441/1‑69 (MQ=255)
ccAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCAct                                                        >  1:1759821/1‑69 (MQ=255)
ccAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCAct                                                        >  1:50475/1‑69 (MQ=255)
   aaaaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTg                                                     <  1:3460044/69‑1 (MQ=255)
   aaaaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTg                                                     <  1:617609/69‑1 (MQ=255)
     aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTGtt                                                   <  1:2846243/69‑1 (MQ=255)
     aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                   <  1:2089916/69‑1 (MQ=255)
     aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                   <  1:1655772/69‑1 (MQ=255)
     aaaTGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGtt                                                   <  1:2085487/69‑1 (MQ=255)
          aaGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTg                                              >  1:1907957/1‑69 (MQ=255)
          aaGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTg                                              >  1:2446962/1‑69 (MQ=255)
                aGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTGTTTACTGGCTGtt                                        <  1:793801/69‑1 (MQ=255)
                aGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGt                                         <  1:1762221/68‑1 (MQ=255)
                 gAACGATGTTGGAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGttt                                                  >  1:2280197/1‑58 (MQ=255)
                 gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGTTTGCACTCTGt                                                    <  1:760134/56‑1 (MQ=255)
                 gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTGt                                                    >  1:891974/1‑56 (MQ=255)
                 gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTg                                       <  1:1371672/69‑1 (MQ=255)
                 gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTg                                       <  1:3716015/69‑1 (MQ=255)
                 gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTg                                       >  1:3437560/1‑69 (MQ=255)
                 gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTg                                       <  1:2289362/69‑1 (MQ=255)
                 gAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGATGCACTCTGt                                                    >  1:2494184/1‑56 (MQ=255)
                        gTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGTTGCACTCTGTTTACTGGCTGTTGCTGGGCt                                <  1:2142694/69‑1 (MQ=255)
                           gCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAg                             >  1:741927/1‑69 (MQ=255)
                            cAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTa                              <  1:1125312/67‑1 (MQ=255)
                               cATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCtt                               >  1:2958065/1‑63 (MQ=255)
                               cATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCtt                               >  1:305971/1‑63 (MQ=255)
                                 ttattCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTg                                    >  1:2289246/1‑56 (MQ=255)
                                    ttCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCg                     >  1:1006046/1‑68 (MQ=255)
                                      cAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATc                  <  1:3183810/69‑1 (MQ=255)
                                      cAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATc                  <  1:1501078/69‑1 (MQ=255)
                                           tCATGGAAAAAAGGACTGGTTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATcaacaa             <  1:2143309/69‑1 (MQ=255)
                                           tCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATcaac               >  1:867100/1‑67 (MQ=255)
                                             aTGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAag            >  1:468074/1‑68 (MQ=255)
                                             aTGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAag            >  1:152692/1‑68 (MQ=255)
                                                gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATcaacaa             >  1:195735/1‑64 (MQ=255)
                                                gAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGaa        >  1:2616100/1‑69 (MQ=255)
                                                    aaaGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAg                             >  1:1085718/1‑44 (MQ=255)
                                                    aaaGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTTCGATCAACAAGAGAATgcg    >  1:353675/1‑69 (MQ=255)
                                                    aaaGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATgcg    >  1:2800203/1‑69 (MQ=255)
                                                      aGGACTGGTTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATgcggc  >  1:1259033/1‑69 (MQ=255)
                                                              |                                                            
CCAAAAAATGAAGGATAGAACGATGTTGCAACATTATTCAGTGTCATGGAAAAAAGGACTGGCTGCACTCTGTTTACTGGCTGTTGCTGGGCTTAGCGGCTGCGATCAACAAGAGAATGCGGC  >  W3110S.gb/1835808‑1835930

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: