Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 631,602 A→G 48.2% *66* (TAA→TAG ybdD → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb631,6020AG48.2% ‑5.0 / 38.0 27*66* (TAA→TAGybdDconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/6);  new base G (6/7);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGT  >  W3110S.gb/631537‑631671
                                                                 |                                                                     
gAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAt                                                                     <  1:2865447/68‑1 (MQ=255)
gAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTg                                                                  <  1:2304092/71‑1 (MQ=255)
    agTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACg                                                              >  1:513230/1‑71 (MQ=255)
                aGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTaaa                                                  >  1:2277877/1‑71 (MQ=255)
                    gCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa                                               <  1:3725709/70‑1 (MQ=255)
                             gcgcTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTa                                                                <  1:2834630/44‑1 (MQ=255)
                                cTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTg                                  <  1:3218426/71‑1 (MQ=255)
                                  aCGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa                                               >  1:988417/1‑56 (MQ=255)
                                   cGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGaa                                <  1:96560/70‑1 (MQ=255)
                                     gTGGAAAAGGCGGGGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGc                             >  1:1193605/1‑71 (MQ=255)
                                        gAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAAc                                                 <  1:3725025/48‑1 (MQ=255)
                                        gAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCt                            >  1:2207084/1‑69 (MQ=255)
                                        gAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTg                           >  1:4008560/1‑70 (MQ=255)
                                         aaaaGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa                                               >  1:3620717/1‑49 (MQ=255)
                                           aaGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACg                       >  1:1799114/1‑71 (MQ=255)
                                           aaGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACg                       >  1:4246167/1‑71 (MQ=255)
                                           aaGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa                                               <  1:2964106/47‑1 (MQ=255)
                                             ggcggcGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCa                     <  1:4360994/71‑1 (MQ=255)
                                              gcggcgCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCaa                    >  1:5075362/1‑71 (MQ=255)
                                               cggcgCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAc                   <  1:4211185/71‑1 (MQ=255)
                                                ggcgCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACg                  <  1:3407895/71‑1 (MQ=255)
                                                 gcgcgcgcTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGt                 <  1:1677446/71‑1 (MQ=255)
                                                       gctgctgctAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCt                            >  1:682914/1‑54 (MQ=255)
                                                            tgctAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTg                                  <  1:4004061/43‑1 (MQ=255)
                                                             gctAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTa                >  1:814548/1‑60 (MQ=255)
                                                             gctAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTa                >  1:820503/1‑60 (MQ=255)
                                                              ctAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGaa                                >  1:4982321/1‑43 (MQ=255)
                                                                aaaTTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGc    >  1:327695/1‑69 (MQ=255)
                                                                 aaTTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGt  <  1:1905221/70‑1 (MQ=255)
                                                                 |                                                                     
GAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGT  >  W3110S.gb/631537‑631671

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: