Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,294,438 T→C 38.6% L395S (TTG→TCG)  yraM → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,294,4380TC38.6% 26.9 / 49.0 44L395S (TTG→TCG) yraMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/12);  new base C (7/10);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TTAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAAT  >  W3110S.gb/3294368‑3294499
                                                                      |                                                             
ttAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGtt                                                               >  1:3258321/1‑71 (MQ=255)
      cAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTg                                                           >  1:2124068/1‑69 (MQ=255)
       aGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTg                                                           >  1:2223577/1‑68 (MQ=255)
        ggTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGaa                                                         >  1:864634/1‑69 (MQ=255)
           tCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGc                                                     >  1:900557/1‑70 (MQ=255)
             agcagGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTg                                                           <  1:832109/62‑1 (MQ=255)
              gcagGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCaa                                                   <  1:3824781/69‑1 (MQ=255)
                   aTGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGaa                                                         >  1:2659820/1‑58 (MQ=255)
                   aTGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGaa                                                         >  1:3169260/1‑58 (MQ=255)
                      gcgcGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCt                                          >  1:1951248/1‑70 (MQ=255)
                          gAGTATTGTGGTCGGTCCGTTGCTGAAGAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACg                                     >  1:4304512/1‑71 (MQ=255)
                          gAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGc                                                     <  1:2973677/55‑1 (MQ=255)
                           aGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGc                                                     <  1:1254888/54‑1 (MQ=255)
                            gTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTa                                   >  1:4524909/1‑71 (MQ=255)
                             tATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCaa                                                   >  1:1855133/1‑54 (MQ=255)
                               ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCaa                                                   >  1:3296714/1‑52 (MQ=255)
                               ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTg                                >  1:2660539/1‑71 (MQ=255)
                               ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTg                                >  1:4185070/1‑71 (MQ=255)
                               ttGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTg                                >  1:248797/1‑71 (MQ=255)
                                tgtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGc                                                     >  1:3015518/1‑49 (MQ=255)
                                tgtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGt                                    <  1:1950125/66‑1 (MQ=255)
                                tgtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGt                                    <  1:1162707/66‑1 (MQ=255)
                                 gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACg                                     <  1:2164713/64‑1 (MQ=255)
                                 gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACg                                     <  1:1985731/64‑1 (MQ=255)
                                 gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACg                                     <  1:5034734/64‑1 (MQ=255)
                                 gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTgg                               <  1:137500/70‑1 (MQ=255)
                                 gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTgg                               <  1:428341/70‑1 (MQ=255)
                                 gtgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTgg                               <  1:4180912/70‑1 (MQ=255)
                                  tgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTc                                              >  1:3854255/1‑54 (MQ=255)
                                  tgGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCa                             >  1:1334837/1‑71 (MQ=255)
                                    gtcggtcCGTTGCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACt                                               >  1:1104613/1‑51 (MQ=255)
                                            gTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGt                                    <  1:315751/54‑1 (MQ=255)
                                            gTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGt                                    <  1:171338/54‑1 (MQ=255)
                                            gTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGt                                    <  1:1505848/54‑1 (MQ=255)
                                            gTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGt                                    <  1:833824/54‑1 (MQ=255)
                                               gCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCg                 >  1:2531871/1‑70 (MQ=255)
                                               gCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCg                 <  1:4660769/70‑1 (MQ=255)
                                               gCTGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAAcc                      >  1:3639278/1‑65 (MQ=255)
                                                 tGAAAAATAACGTTGAAGAGTCGCTGAAGAGCAACACTCCGCTg                                         <  1:4275863/44‑1 (MQ=255)
                                                   aaaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgag              <  1:4483456/69‑1 (MQ=255)
                                                   aaaaaTAACGTTGAAGAGTCTCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgaga             <  1:973355/70‑1 (MQ=255)
                                                    aaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGgaga             >  1:4616745/1‑69 (MQ=255)
                                                    aaaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGaa            >  1:926759/1‑70 (MQ=255)
                                                     aaaTAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAata          <  1:3475215/71‑1 (MQ=255)
                                                             ttGAAGAGTCGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAAt  <  1:1106274/71‑1 (MQ=255)
                                                                      |                                                             
TTAAGCCAGGTTCAGCAGGATGGCGCGAGTATTGTGGTCGGTCCGTTGCTGAAAAATAACGTTGAAGAGTTGCTGAAGAGCAACACTCCGCTGAACGTACTGGCACTGAACCAGCCGGAGAATATCGAAAAT  >  W3110S.gb/3294368‑3294499

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: