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% 86.4 / NA 23intergenic (+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 (12/11);  total (12/11)

TCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCAAA  >  AM260479/3247961‑3248218
                                                                                                                                    |                                                                                                                              
tcgtcgCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACAcc                                                                                                                          >  1:171229/1‑134 (MQ=255)
tcgtcgCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACAcc                                                                                                                          >  1:32789/1‑134 (MQ=255)
            tCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTg                                                                                                              <  2:933379/134‑1 (MQ=255)
                                 gcgaccGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAc                                                                                            >  1:442519/5‑134 (MQ=255)
                                  cgaccGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc                                                                                           <  1:898450/131‑1 (MQ=255)
                                             ctGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGt                                                                                 <  2:171229/133‑1 (MQ=255)
                                                    gCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAg                                                                          >  1:128401/1‑134 (MQ=255)
                                                        ctcGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCa                                                                      >  2:45308/1‑134 (MQ=255)
                                                                                       gACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGc                                                   <  2:994558/123‑1 (MQ=255)
                                                                                       gACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGc                                                   >  1:994558/1‑123 (MQ=255)
                                                                                                    gtcaagcagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc                           <  2:128401/130‑1 (MQ=255)
                                                                                                        agcagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCgccgcc                       <  2:442519/134‑1 (MQ=255)
                                                                                                            gggTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAacac                                                        >  2:192547/1‑97 (MQ=255)
                                                                                                            gggTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAacac                                                        <  1:192547/97‑1 (MQ=255)
                                                                                                                     ccTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCtatata           >  2:514500/1‑133 (MQ=255)
                                                                                                                        gcgcAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgccg                         >  1:629573/1‑116 (MQ=255)
                                                                                                                        gcgcAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgccg                         <  2:629573/116‑1 (MQ=255)
                                                                                                                         cgcAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAg      <  1:514500/134‑1 (MQ=255)
                                                                                                                             gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACg                                              <  2:478143/90‑1 (MQ=255)
                                                                                                                             gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACg                                              >  1:478143/1‑90 (MQ=255)
                                                                                                                             gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc                           >  1:87357/1‑109 (MQ=255)
                                                                                                                             gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc                           <  2:87357/109‑1 (MQ=255)
                                                                                                                             gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCaaa  >  2:984620/1‑134 (MQ=255)
                                                                                                                                    |                                                                                                                              
TCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCAAA  >  AM260479/3247961‑3248218

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