Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 582,589 T→G H295Q (CAT→CAG h16_A0555 → zinc metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479582,5890TG100.0% 41.6 / NA 12H295Q (CAT→CAGh16_A0555zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (7/5);  total (7/5)

CCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCATGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAG  >  AM260479/582474‑582703
                                                                                                                   |                                                                                                                  
ccATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATttctt                                                                                                   >  2:968454/1‑133 (MQ=255)
                                                  ggggCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTa                                                <  2:269174/134‑1 (MQ=255)
                                                    ggCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTccc                                                      <  1:420879/126‑1 (MQ=255)
                                                    ggCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTccc                                                      >  2:420879/1‑126 (MQ=255)
                                                                  ccGACCGGGTCGGCTTCCAGATCCTGACAGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCg                                >  1:319219/1‑134 (MQ=255)
                                                                              gCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAg                    >  1:1036569/1‑134 (MQ=255)
                                                                                    aGATCCTGACAGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGAtt              <  2:319219/134‑1 (MQ=255)
                                                                                               gcggcgGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTAtg                                              >  1:184394/1‑91 (MQ=255)
                                                                                               gcggcgGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTAtg                                              <  2:184394/91‑1 (MQ=255)
                                                                                               gcggcgGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGAtt              >  1:488104/1‑123 (MQ=255)
                                                                                               gcggcgGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGAtt              <  2:488104/123‑1 (MQ=255)
                                                                                                cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAg  >  1:1116456/1‑134 (MQ=255)
                                                                                                                   |                                                                                                                  
CCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCATGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAG  >  AM260479/582474‑582703

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