Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,816,575:1 +G 28.6% coding (995/2262 nt) katE → hydroperoxidase HPII(III)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,816,5751.G28.6% 12.7 / 14.1 14coding (995/2262 nt)katEhydroperoxidase HPII(III)
Reads supporting (aligned to +/- strand):  ref base . (7/3);  new base G (2/2);  total (9/5)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCT  >  W3110S.gb/1816506‑1816638
                                                                      |                                                               
aaCTCACCGTACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑g                                                                <  1:3122074/71‑1 (MQ=255)
aaCTCACCGGACGTGACCCGGCCTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑g                                                                >  1:335938/1‑71 (MQ=255)
aaCTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑g                                                                <  1:1889995/71‑1 (MQ=255)
 aCTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑gg                                                               <  1:2159789/71‑1 (MQ=255)
 aCTCACCGGACGTGACCCGGACTTCCACCGCAGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑g                                                                <  1:2564468/70‑1 (MQ=255)
                                    aGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTgaagaa                            <  1:2042957/71‑1 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑ggg                                                              <  1:2487647/35‑2 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑ggg                                                              <  1:670925/35‑2 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑ggg                                                              <  1:371075/35‑2 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑ggg                                                              <  1:3579953/35‑2 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑ggg                                                              <  1:1798565/35‑2 (MQ=255)
                                        gtgGGAAGCCATTGAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGAt                         >  1:1940798/1‑70 (MQ=255)
                                              aggcATTGAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTc                   >  1:1743900/4‑70 (MQ=255)
                                              aGCCATTGAAGCAGGCGATTTTGC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTc                   >  1:2894108/1‑70 (MQ=255)
                                                 ctttGAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAACAAGATGAATTc                   <  1:2234960/65‑1 (MQ=255)
                                                   ttGAAGCAGGCGATTTTCCGGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGa                       <  1:1085389/62‑1 (MQ=255)
                                                     gAAGCAGGCGATTTTCCGGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAAt                     >  1:3470908/1‑62 (MQ=255)
                                                     gAAGCAGGCGATTTTCCGGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAAt                     >  1:1097708/1‑62 (MQ=255)
                                                     gAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCg            <  1:1103787/70‑1 (MQ=255)
                                                     gAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGa           >  1:3181244/1‑71 (MQ=255)
                                                     gAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGa           >  1:1172048/1‑71 (MQ=255)
                                                     gAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGa           >  1:3808610/1‑71 (MQ=255)
                                                              gATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGAtct  >  1:1221464/1‑71 (MQ=255)
                                                               aTTTTCCGGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGAtct  <  1:3345475/71‑1 (MQ=255)
                                                                      |                                                               
AACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCC‑GGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCT  >  W3110S.gb/1816506‑1816638

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: ‑