Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,248,092 G→C intergenic (+43/‑175) h16_A3001 → / → hemK1 transcriptional regulator, LysR‑family/methylase of polypeptide chain release factors

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,248,0920GC100.0% 88.2 / NA 24intergenic (+43/‑175)h16_A3001/hemK1transcriptional regulator, LysR‑family/methylase of polypeptide chain release factors
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (10/14);  total (10/14)

TCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCG  >  AM260479/3247961‑3248195
                                                                                                                                    |                                                                                                       
tcgtcgCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACAcc                                                                                                   >  1:224394/1‑134 (MQ=255)
     gCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCAt                                                                                              <  2:224394/134‑1 (MQ=255)
      ccGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATc                                                                                             >  1:253814/1‑134 (MQ=255)
          cGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCCGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGc                                                                                         <  2:440217/134‑1 (MQ=255)
                            gagcgcgACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCt                                                                         <  2:253814/131‑1 (MQ=255)
                                                   ggCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCa                                                    <  1:58222/134‑1 (MQ=255)
                                                          cggccggcCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc                                                                    <  2:799207/111‑1 (MQ=255)
                                                          cggccggcCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc                                                                    >  1:799207/1‑111 (MQ=255)
                                                            gccggccACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTc                                           <  2:378711/134‑1 (MQ=255)
                                                                            gcagcaACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc                           <  1:829923/134‑1 (MQ=255)
                                                                                gcaacGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTAc                        >  1:468535/5‑134 (MQ=255)
                                                                                       gACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGtt                                                         >  1:697011/1‑94 (MQ=255)
                                                                                       gACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGtt                                                         <  2:697011/94‑1 (MQ=255)
                                                                                        aCGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACAcc                                                                                                   <  2:958292/51‑1 (MQ=38)
                                                                                        aCGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACAcc                                                                                                   >  1:958292/1‑51 (MQ=38)
                                                                                          gCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAg              <  2:468535/134‑1 (MQ=255)
                                                                                                 gCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGc                                                              >  1:1027368/1‑79 (MQ=255)
                                                                                                 gCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGc                                                              <  2:1027368/79‑1 (MQ=255)
                                                                                                    gtcaagcagGGTGCCTACCTGCGCGGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc    >  2:44554/5‑134 (MQ=255)
                                                                                                    gtcaagcagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc    <  2:689008/130‑1 (MQ=255)
                                                                                                                           cAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCa                                                                       <  2:492363/44‑1 (MQ=255)
                                                                                                                           cAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCa                                                                       >  1:492363/1‑44 (MQ=255)
                                                                                                                             gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgccg  >  2:34240/1‑111 (MQ=255)
                                                                                                                             gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgccg  <  1:34240/111‑1 (MQ=255)
                                                                                                                                    |                                                                                                       
TCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCG  >  AM260479/3247961‑3248195

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: