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% 111.7 / NA 29H295Q (CAT→CAGh16_A0555zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/14);  total (15/14)

GCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCATGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATT  >  AM260479/582473‑582691
                                                                                                                    |                                                                                                      
gCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATttc                                                                                         <  1:379079/132‑1 (MQ=255)
                                       tCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                 <  2:676608/133‑1 (MQ=255)
                                           tctcGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTc                                            >  2:959936/1‑134 (MQ=255)
                                                  cgGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCt                                      >  2:543253/1‑133 (MQ=255)
                                                     ggCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTAtg                                  <  1:959936/134‑1 (MQ=255)
                                                        cGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGc                               >  1:1048970/1‑134 (MQ=255)
                                                             gcgAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACt                          >  2:874211/1‑134 (MQ=255)
                                                                     gACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCt                  >  2:811485/1‑134 (MQ=255)
                                                                                 ttCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCg      >  2:265317/1‑134 (MQ=255)
                                                                                     aGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGAtt  <  1:874211/134‑1 (MQ=255)
                                                                                     aGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGAtt  <  2:308176/134‑1 (MQ=255)
                                                                                         ccTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                 >  1:186317/1‑83 (MQ=255)
                                                                                         ccTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                 <  2:186317/83‑1 (MQ=255)
                                                                                             aCCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGTAGCGCGCATCGGGGACcg                                                          >  2:867692/1‑70 (MQ=255)
                                                                                             aCCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGTAGCGCGCATCGGGGACcg                                                          <  1:867692/70‑1 (MQ=255)
                                                                                             aCCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGAcc              >  1:940833/1‑114 (MQ=255)
                                                                                             aCCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGAcc              <  2:940833/114‑1 (MQ=255)
                                                                                              ccGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                 <  1:568497/78‑1 (MQ=255)
                                                                                              ccGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                 >  2:568497/1‑78 (MQ=255)
                                                                                                gcggcgGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTAtg                                  <  2:433442/91‑1 (MQ=255)
                                                                                                gcggcgGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTAtg                                  >  1:433442/1‑91 (MQ=255)
                                                                                                 cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                 >  1:598838/1‑75 (MQ=255)
                                                                                                 cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                 <  2:598838/75‑1 (MQ=255)
                                                                                                 cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTAtgtg                                <  2:303082/92‑1 (MQ=255)
                                                                                                 cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTAtgtg                                <  2:507526/92‑1 (MQ=255)
                                                                                                 cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTAtgtg                                >  1:507526/1‑92 (MQ=255)
                                                                                                 cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTAtgtg                                >  1:303082/1‑92 (MQ=255)
                                                                                                 cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGAtt  >  1:773524/1‑122 (MQ=255)
                                                                                                 cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGAtt  <  2:773524/122‑1 (MQ=255)
                                                                                                                    |                                                                                                      
GCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCATGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATT  >  AM260479/582473‑582691

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: