Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,036,723 G→A 37.5% A184T (GCT→ACT)  ftsY → fused Signal Recognition Particle (SRP) receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,036,7230GA37.5% 27.1 / 31.0 32A184T (GCT→ACT) ftsYfused Signal Recognition Particle (SRP) receptor
Reads supporting (aligned to +/- strand):  ref base G (9/11);  new base A (5/7);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

CCGAAGCGGCAGAAGAGGCGGTGATGGTGGTTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGCCTGTTAAAAACC  >  W3110S.gb/4036655‑4036791
                                                                    |                                                                    
ccGAAGCGGCAGAAGAGGCGGTGATGGTGGTTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCGc                                                                     <  1:1683887/70‑1 (MQ=255)
  gAAGCGGCAGAAGAGGCGGTGATGGTGGTTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCACt                                                                    >  1:4021748/1‑69 (MQ=255)
                      gATGGTGGTTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCGCTCAGGAGCAGGAAAAACCGAcc                                               >  1:950380/1‑70 (MQ=255)
                      caTGGTGGTTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCGCTCAGGAGCAGGAAAAACCGAcc                                               >  1:3848188/2‑70 (MQ=255)
                       aTGGTGGTTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCACTCAGGAGCAGGAAAAACCGACCaa                                             <  1:3135866/71‑1 (MQ=255)
                        tggtggTTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCACTCAGGAGCAGGAAAAACCGACCaaa                                            <  1:1893644/71‑1 (MQ=255)
                          gtggtTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCACTCAGGAGCAGGAAAAACCGACCAaaga                                          <  1:3008818/71‑1 (MQ=255)
                           tgggTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCACTCAGGAGCAGGAAAAACCGACCAaaga                                          <  1:91909/70‑1 (MQ=255)
                             gtTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCGCTCAGGAGCAGGAAAAACCGACCAaagaa                                         >  1:437451/1‑69 (MQ=255)
                                cctccGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCACt                                                                    <  1:4401237/39‑1 (MQ=255)
                                cctccGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCACt                                                                    <  1:3956662/39‑1 (MQ=255)
                                   ccGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCGCTCAGGAGCAGGAAAAACCGACCAaag                                           >  1:2354974/1‑61 (MQ=255)
                                       cAGAAGAAGAGCAGCCGGTGGAAGAAATCGCTCAGGAGCAGGaaa                                                       <  1:2177632/45‑1 (MQ=255)
                                           agaagaGCAGCCGGTGGAAGAAATCACTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGcc                         >  1:553094/1‑71 (MQ=255)
                                              agagCAGCCGGTGGAAGAAATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGa                      <  1:487948/71‑1 (MQ=255)
                                                  cagcCGGTGGAAGAAATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGttttt                                  <  1:2593332/55‑1 (MQ=255)
                                                    gcCGGTGGAAGAAATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCa                <  1:1665310/71‑1 (MQ=255)
                                                       ggTGGAAGAAATCACTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGcc             >  1:1095138/1‑71 (MQ=255)
                                                       ggTGCAAGAAATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGcc             >  1:1704696/1‑71 (MQ=255)
                                                           gaagaaATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGCCTGt          >  1:2680586/1‑70 (MQ=255)
                                                            aagaaATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGCCTGTTa        <  1:2594359/71‑1 (MQ=255)
                                                            aagaaATCACTCAGGAGCAGGAAAAACCGACCAaagaa                                         <  1:2907506/38‑1 (MQ=255)
                                                             agaaATCACTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGCCTGTTa        >  1:3621432/1‑70 (MQ=255)
                                                              gaaATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGCCTGTTaaa      <  1:2804656/71‑1 (MQ=255)
                                                              gaaATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGCCTGTTaaa      <  1:1480804/71‑1 (MQ=255)
                                                               aaaTCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAAGGCAGCCTGTTaaa      >  1:2800461/1‑70 (MQ=255)
                                                               aaaTCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCa                <  1:5057443/60‑1 (MQ=255)
                                                                aaTCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGttt                                    <  1:4021703/39‑1 (MQ=255)
                                                                 aTCACTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGCCTGTTAAAAAc   >  1:215295/1‑71 (MQ=255)
                                                                  tCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACTCAGCCTGTTAAAAAcc  <  1:4192899/71‑1 (MQ=255)
                                                                   cGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGtttttt                                 >  1:1732623/1‑39 (MQ=255)
                                                                   cGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGtttttt                                 >  1:4068397/1‑39 (MQ=255)
                                                                    |                                                                    
CCGAAGCGGCAGAAGAGGCGGTGATGGTGGTTCCTCCGGCAGAAGAAGAGCAGCCGGTGGAAGAAATCGCTCAGGAGCAGGAAAAACCGACCAAAGAAGGTTTTTTCGCGCGCCTGAAACGCAGCCTGTTAAAAACC  >  W3110S.gb/4036655‑4036791

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: