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.5 / 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)

GTTGCGCTGAGCGCGGCAGGCGGGCGGCGGCGGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCGTCAGCGCCTGTCTACATC  >  AM260480/2017982‑2018217
                                                                                                                        |                                                                                                                    
gTTGCGCTGAGCGCGGCAGGCGGGCGGCGGCGGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATgcgcccgga                                                                                                          <  1:458261/133‑6 (MQ=255)
                               gggcgCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAgcgc                                                                          >  1:651664/1‑134 (MQ=255)
                                         gCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTAACAGCGt                                                                 >  1:985887/1‑133 (MQ=255)
                                            tGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTAACAGCGTCATg                                                             <  2:985887/134‑1 (MQ=255)
                                                                    cgccCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAgcg                                                                           >  1:861659/1‑96 (MQ=255)
                                                                    cgccCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAgcg                                                                           <  2:861659/96‑1 (MQ=255)
                                                                               cGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGa                                                                                   >  1:194404/1‑77 (MQ=255)
                                                                               cGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGa                                                                                   <  2:194404/77‑1 (MQ=255)
                                                                                     cGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTc                                                                >  1:825677/1‑90 (MQ=255)
                                                                                     cGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTc                                                                <  2:825677/90‑1 (MQ=255)
                                                                                       cTTATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCGTc                  <  2:619027/134‑1 (MQ=255)
                                                                                         tATCCGGGCGCATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCGTCAg                >  1:851161/1‑134 (MQ=255)
                                                                                                  gcATGGTCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCGTCAGCGCCTGTCt       <  2:651664/134‑1 (MQ=255)
                                                                                                        tCCGGGGCATGCCCGACCATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCGTCAGCGCCTGTCTACATc  >  1:111829/1‑133 (MQ=255)
                                                                                                                        |                                                                                                                    
GTTGCGCTGAGCGCGGCAGGCGGGCGGCGGCGGGCGCTGAAGCATGGCAATGCCTTCAGCGCCCGAGCCGCCCTTGCTTCGTGCCCGCTTATCCGGGCGCATGGTCCGGGGCATGCCCGACATGCGCACCGATGGCCAGCGCCGCCCAGCTCCGATGCCAGCGCTGACAGCGTCATGGCAGCCAGTATCAGCAATCTGCGCATTGCGCCCCTCATGCGTCAGCGCCTGTCTACATC  >  AM260480/2017982‑2018217

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