Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 2,018,102 +C intergenic (+111/+98) h16_B1778 → / ← h16_B1779 putative peptidase/putative peptidase, S33 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604802,018,1011.C100.0% 50.0 / NA 14intergenic (+110/+99)h16_B1778/h16_B1779putative peptidase/putative peptidase, S33 family
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (7/7);  total (7/7)

CAGCGTTCAGGTCGTTGCGCTGAGCGCGGCAGGCGGGCGGCGGCGGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCGTC  >  AM260480/2017969‑2018201
                                                                                                                                     |                                                                                                    
cAGCGTTCAGGTCGTTGCGCTGAGCGCGGCAGGCGGGCGGCGGCGGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGAc                                                                                                     <  1:246712/134‑1 (MQ=255)
     ttCAGGTCGTTGCGCTGAGCGCGGCAGGCGGGCGGCGGCGGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGAccatgc                                                                                                >  1:495193/1‑129 (MQ=255)
                 cgcTGAGCGCGGCAGGCGGGCGGCGGCGGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGc                                                                                      >  2:879609/1‑133 (MQ=255)
                   cTGAGCGCGGCAGGCGGGCGGCGGCGGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCACTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAg                                                                                   >  1:699238/1‑134 (MQ=255)
                                    gcggcggcgGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGAt                                                                  >  1:238228/1‑134 (MQ=255)
                                           cgggcgCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGcc                                                                                     <  1:398526/108‑1 (MQ=255)
                                           cgggcgCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGcc                                                                                     >  2:398526/1‑108 (MQ=255)
                                                           gCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGc                                           <  2:495193/134‑1 (MQ=255)
                                                             aaTGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAg                                         >  2:465063/1‑134 (MQ=255)
                                                              aTGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGc                                        >  2:404999/1‑134 (MQ=255)
                                                                     cAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTAt                                  <  1:419019/133‑1 (MQ=255)
                                                                                 cgccCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTgcgc                     <  1:879609/134‑1 (MQ=255)
                                                                                                  cGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGAc                                                      <  1:541518/84‑1 (MQ=255)
                                                                                                  cGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGAc                                                      >  2:541518/1‑84 (MQ=255)
                                                                                                  cGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCg    <  2:150601/134‑1 (MQ=255)
                                                                                                    cTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCGTc  <  1:404999/134‑1 (MQ=255)
                                                                                                                                     |                                                                                                    
CAGCGTTCAGGTCGTTGCGCTGAGCGCGGCAGGCGGGCGGCGGCGGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCGTC  >  AM260480/2017969‑2018201

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: