Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 581,821 +C coding (117/1734 nt) h16_A0555 → zinc metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479581,8201.C100.0% 53.9 / NA 15S39T (AGC→ACC) h16_A0555zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (7/8);  total (7/8)

GCCGGCCATTGCCGCGCTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTCCGA  >  AM260479/581697‑581943
                                                                                                                            |                                                                                                                           
gcCGGCCATTGCCGCGCTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGt                                                                                                                    >  2:113200/1‑134 (MQ=255)
                                        ggcggcCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGccc                                                                            >  2:1165277/1‑134 (MQ=255)
                                                            ctgGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAg                                                        <  1:113200/134‑1 (MQ=255)
                                                                      cGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTc                                              <  1:251347/134‑1 (MQ=255)
                                                                       gTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCt                                             <  1:1165277/134‑1 (MQ=255)
                                                                               cgcCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGc                                                                                                                            >  1:733516/1‑46 (MQ=255)
                                                                               cgcCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGc                                                                                                                            <  2:733516/46‑1 (MQ=255)
                                                                                     ccTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAAc                               >  2:710339/1‑134 (MQ=255)
                                                                                             cAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTAt                                           >  1:250848/1‑114 (MQ=255)
                                                                                             cAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTAt                                           <  2:250848/114‑1 (MQ=255)
                                                                                                  tACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGc                                                       >  2:1098172/1‑97 (MQ=255)
                                                                                                  tACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGc                                                       <  1:1098172/97‑1 (MQ=255)
                                                                                                      ccgcGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCggg              <  1:1042254/134‑1 (MQ=255)
                                                                                                        gcgcACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGc                                                       <  1:381213/91‑1 (MQ=255)
                                                                                                        gcgcACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGc                                                       >  2:381213/1‑91 (MQ=255)
                                                                                                              ggCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGt      >  2:722995/1‑134 (MQ=255)
                                                                                                                  ccGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTCCGa  <  1:710339/134‑1 (MQ=255)
                                                                                                                            |                                                                                                                           
GCCGGCCATTGCCGCGCTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTCCGA  >  AM260479/581697‑581943

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