Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,248,064 +A intergenic (+15/‑203) 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,0631.A100.0% 90.6 / NA 25intergenic (+14/‑204)h16_A3001/hemK1transcriptional regulator, LysR‑family/methylase of polypeptide chain release factors
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (11/14);  total (11/14)

AGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAG  >  AM260479/3247928‑3248183
                                                                                                                                        |                                                                                                                        
aGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTcaag                                                                                                                       >  2:811791/1‑132 (MQ=255)
   cgtacgtGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTGGGCAGCAACGTGACGCACCTGGCCGTcaagcag                                                                                                                    <  2:373917/134‑6 (MQ=255)
           gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTa                                                                                                             >  1:378711/1‑134 (MQ=255)
           gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCCGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTa                                                                                                             >  1:440217/1‑134 (MQ=255)
                  ggCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcgc                                                                                                      <  1:811791/134‑1 (MQ=255)
                                 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)
                                                                   cgaccGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCa                                                                                                     <  1:858804/86‑1 (MQ=255)
                                                                   cgaccGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCa                                                                                                     >  2:858804/4‑89 (MQ=255)
                                                                                    ggCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCa                                        <  1:58222/134‑1 (MQ=255)
                                                                                           cggccggcCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc                                                        >  1:799207/1‑111 (MQ=255)
                                                                                           cggccggcCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc                                                        <  2:799207/111‑1 (MQ=255)
                                                                                             gccggccACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTc                               <  2:378711/134‑1 (MQ=255)
                                                                                                             gcagcaACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc               <  1:829923/134‑1 (MQ=255)
                                                                                                                 gcaacGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTAc            >  1:468535/5‑134 (MQ=255)
                                                                                                                     cGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcgc                                                                                                      <  1:508293/40‑1 (MQ=37)
                                                                                                                     cGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcgc                                                                                                      >  2:508293/1‑40 (MQ=37)
                                                                                                                        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                                                  <  2:1027368/79‑1 (MQ=255)
                                                                                                                                  gCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGc                                                  >  1:1027368/1‑79 (MQ=255)
                                                                                                                                        |                                                                                                                        
AGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAG  >  AM260479/3247928‑3248183

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