Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,195,173 A→G 25.0% *118* (TAA→TGA)  rplR → 50S ribosomal subunit protein L18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,195,1730AG25.0% 30.5 / 9.7 20*118* (TAA→TGA) rplR50S ribosomal subunit protein L18
Reads supporting (aligned to +/- strand):  ref base A (4/11);  new base G (3/2);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGGTTCCAATATCATGGTCGTGTCCAGGCACTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGGTGTAAGATGGCTCACATCGAAAAACAAGCTGGCGAACTGCAGGAAAAGCTGATCGC  >  W3110S.gb/4195104‑4195238
                                                                     |                                                                 
ccGTGTTCCAATATCATGGTCGTGTCCAGGCACTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTa                                                                   <  1:1702309/70‑1 (MQ=255)
    gTTCCAATATCATGGTCGTGTCCAGGCACTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTa                                                              <  1:2509590/71‑1 (MQ=255)
           tatCATGGTCGTGTCCAGGCACTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTGAGGTAGAGgtgt                                                       <  1:751168/71‑1 (MQ=255)
            atCATGGTCGTGTCCAGGCACTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGgtgt                                                       >  1:58469/1‑70 (MQ=255)
             tCATGGTCGTGTCCAGGCACTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGGTGTaa                                                     <  1:413345/71‑1 (MQ=255)
                   tCGTGTCCAGGCACTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTagag                                                           <  1:1229128/59‑1 (MQ=255)
                   tCGTGTCCAGGCACTGGCAGATGCTGCCCGTGAAGCTGGACTTCAGTTCTAAGGTagag                                                           <  1:233695/59‑1 (MQ=255)
                          cAGGCACTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTGAGGTagag                                                           >  1:858865/1‑52 (MQ=255)
                               aCTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGGTGTAAGATGGCTCACATCGaaaa                                   <  1:831580/71‑1 (MQ=255)
                                cTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGGTGTAAGATgg                                                <  1:1664180/57‑1 (MQ=255)
                                  ggCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTGAGGTAGAGGTGTAAGATGGCTCACATCGAAAAAc                                 >  1:1625046/1‑70 (MQ=255)
                                  ggCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGGTGTAAGATGGCTCACATc                                        >  1:1952934/1‑63 (MQ=255)
                                   gCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGGTGTAAGATGGCTCACATCGAAAAACa                                <  1:1903889/70‑1 (MQ=255)
                                    cAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTGAGGTAGAGGTGTAAGATGGCTCACATCGAAAAACaa                               <  1:1300935/70‑1 (MQ=255)
                                         gctgcCCGTGAAGCTGGCCTTCAGTTCTGAGGTAGAGGTGTAAGATGGCTCACATCGAAAAACAAGCTGGc                         >  1:1889255/1‑71 (MQ=255)
                                         gctgcCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGGTGTAAGATGGCTCACATCGAAAAACAAGCTGGc                         >  1:115545/1‑71 (MQ=255)
                                           tgcCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGGTGTAAGATGGCTCACATCGa                                      <  1:278089/56‑1 (MQ=255)
                                                   aaGCTGGCCTTCAGTTCTAAGGTAGAGGTGTAAGATGGCTCACATCGAAAAACAAg                              >  1:1075164/1‑56 (MQ=255)
                                                         gCCTTCAGTTCTAAGGTAGAGGTGTAAGATGGCTCACATCGAAAAAc                                 <  1:1646764/47‑1 (MQ=255)
                                                                gTTCTAAGGTAGAGGTGTAAGATGGCTCACATCGAAAAACAAGCTGGCGAACTGCAGGAAAAGCTGATcgc  <  1:63852/71‑1 (MQ=255)
                                                                  tCTAAGGTAGAGGTGTAAGATGGCTCACATCGAAAAACAAGCTGGCGAACTGCAGGAAAAGCTg       <  1:1742106/64‑1 (MQ=255)
                                                                     |                                                                 
CCGGGTTCCAATATCATGGTCGTGTCCAGGCACTGGCAGATGCTGCCCGTGAAGCTGGCCTTCAGTTCTAAGGTAGAGGTGTAAGATGGCTCACATCGAAAAACAAGCTGGCGAACTGCAGGAAAAGCTGATCGC  >  W3110S.gb/4195104‑4195238

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: