Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,416,062 A→G 25.8% N27D (AAT→GAT)  yjfK → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,416,0620AG25.8% 50.2 / 17.9 31N27D (AAT→GAT) yjfKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (12/11);  new base G (5/3);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCG  >  W3110S.gb/4416005‑4416130
                                                         |                                                                    
cTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATa                                                                   >  1:1948603/1‑61 (MQ=255)
 tGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTa                                                          >  1:695072/1‑69 (MQ=255)
 tGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTa                                                          >  1:625253/1‑69 (MQ=255)
   tttGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTa                                                          >  1:2063590/1‑67 (MQ=255)
   tttGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTa                                                          <  1:524708/67‑1 (MQ=255)
   tttGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTAc                                                         <  1:2313264/68‑1 (MQ=255)
       gCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATa                                                                   <  1:1033880/54‑1 (MQ=255)
         aaGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCaa                                                                     <  1:1531189/50‑1 (MQ=255)
                 taaGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgt                                        >  1:1759041/1‑71 (MQ=255)
                     ccATCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCt                                                             <  1:2458457/46‑1 (MQ=255)
                      cAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTc                                                     <  1:1848464/53‑1 (MQ=255)
                        gCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgtttcgttt                                 <  1:482342/71‑1 (MQ=255)
                          tATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGc                               >  1:696672/1‑71 (MQ=255)
                              gCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTgg                            >  1:2110407/1‑70 (MQ=255)
                               cTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTc                                                     >  1:1694785/1‑44 (MQ=255)
                               cTCGTGGTCCGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTg                             >  1:2187814/1‑68 (MQ=255)
                                        cGCTGGGACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAAtt                  <  1:1978512/70‑1 (MQ=255)
                                          cTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGaa                    >  1:1174170/1‑66 (MQ=255)
                                          cTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGaa                    >  1:1100478/1‑66 (MQ=255)
                                          cTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCt               <  1:1323653/71‑1 (MQ=255)
                                             ggACTTCATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTg              <  1:1152141/69‑1 (MQ=255)
                                               aCTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTgg                            >  1:42314/1‑53 (MQ=255)
                                               aCTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTg          <  1:2334385/71‑1 (MQ=255)
                                                 ttCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAAt                   <  1:555666/60‑1 (MQ=255)
                                                   cATCTCGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGaa                          >  1:958999/1‑51 (MQ=255)
                                                   cATCTCAATAGTGGCTTTACGCTCGACACGTTAGcgtttcg                                    <  1:1009104/41‑1 (MQ=255)
                                                   cATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTg                                >  1:1829299/1‑45 (MQ=255)
                                                   cATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTg                                >  1:1571989/1‑45 (MQ=255)
                                                     tctcGATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGc                               <  1:1081264/44‑1 (MQ=255)
                                                     tctcAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCt               <  1:1197752/60‑1 (MQ=255)
                                                      ctcAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGcc   >  1:1709798/1‑71 (MQ=255)
                                                       tcAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCg  >  1:411876/1‑71 (MQ=255)
                                                         aaTAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATg                      >  1:2108220/1‑49 (MQ=255)
                                                         |                                                                    
CTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCG  >  W3110S.gb/4416005‑4416130

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