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

AACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCT  >  W3110S.gb/1816506‑1816638
                                                                      |                                                               
aaCTCACCGTACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCg                                                                <  1:3122074/71‑1 (MQ=255)
aaCTCACCGGACGTGACCCGGCCTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCg                                                                >  1:335938/1‑71 (MQ=255)
aaCTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCg                                                                <  1:1889995/71‑1 (MQ=255)
 aCTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCgg                                                               <  1:2159789/71‑1 (MQ=255)
 aCTCACCGGACGTGACCCGGACTTCCACCGCAGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCg                                                                <  1:2564468/70‑1 (MQ=255)
                                    aGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTgaagaa                            <  1:2042957/71‑1 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCCggg                                                              <  1:2487647/35‑2 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCCggg                                                              <  1:670925/35‑2 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCCggg                                                              <  1:371075/35‑2 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCCggg                                                              <  1:3579953/35‑2 (MQ=255)
                                      ttGTGGGAAGCCATTGAAGCAGGCGATTTTCCggg                                                              <  1:1798565/35‑2 (MQ=255)
                                        gtgGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGAt                         >  1:1940798/1‑70 (MQ=255)
                                              aggcATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTc                   >  1:1743900/4‑70 (MQ=255)
                                              aGCCATTGAAGCAGGCGATTTTGCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTc                   >  1:2894108/1‑70 (MQ=255)
                                                 ctttGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAACAAGATGAATTc                   <  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)
                                                     gAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCg            <  1:1103787/70‑1 (MQ=255)
                                                     gAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGa           >  1:3181244/1‑71 (MQ=255)
                                                     gAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGa           >  1:1172048/1‑71 (MQ=255)
                                                     gAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGa           >  1:3808610/1‑71 (MQ=255)
                                                              gATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGAtct  >  1:1221464/1‑71 (MQ=255)
                                                               aTTTTCCGGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGAtct  <  1:3345475/71‑1 (MQ=255)
                                                                      |                                                               
AACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCT  >  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: