Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 2,295,879 T→A Q603L (CAG→CTG)  h16_B2042 ← Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604802,295,8790TA100.0% 50.6 / NA 15Q603L (CAG→CTG) h16_B2042Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (6/9);  total (6/9)

GTTGCTGGTCGCGGGCCTGGTCCTGCCACTGGCGTTGCTGGGCGTCGCGTTGCCGGGCCTCCTGCTGTTGCCGCTGCAGATCAGGCGCCTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCTGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGCCGGTCCGCCCGGCTGGCCCGGCTGGCCAACAAAGCCGCGGGCTTCGTCGGGCCGGC  >  AM260480/2295748‑2296008
                                                                                                                                   |                                                                                                                                 
gTTGCTGGTCGCGGGCCTGGTCCTGCCACTGGCGTTGCTGGGCGTCGCGTTGCCGGGCCTCCTGCTGTTGCCGCTGCAGATCAGGCGCCTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCgctgcagc                                                                                                                                 <  2:145565/134‑1 (MQ=255)
                      cTGCCACTGGCGTTGCTGGGCGTCGCGTTGCCGGGCCTCCTGCTGTTGCCGCTGCAGATCAGGCGCCTGCCGCTGCTGCACGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCaaaa                                                                                                           <  2:928652/134‑1 (MQ=255)
                                 gTTGCTGGGCGTCGCGTTGCCGGGCCTCCTGCTGTTGCCGCTGCAGATCAGGCGCCTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGt                                                                                                <  1:890172/134‑1 (MQ=255)
                                  ttGCTGGGCGTCGCGTTGCCGGGCCTCCTGCTGTTGCCGCTGCAGATCAGGCGCCTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGtt                                                                                               >  2:64535/1‑134 (MQ=255)
                                                         cctcctGCTGTTGCCGCTGCAGATCAGGCGCCTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGcgcccgc                                                                        <  1:948844/134‑1 (MQ=255)
                                                         cctcctGCTGTTGCCGCTGAAGATCAGGCGCCTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGcgcccgc                                                                        <  1:897534/134‑1 (MQ=255)
                                                                            cAGATCAGGCGCCTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGCCgg                                                       <  1:315613/132‑1 (MQ=255)
                                                                            cAGATCAGGCGCCTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGCCGGTc                                                     >  1:910299/1‑134 (MQ=255)
                                                                                        cTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGc                                                          <  1:229217/117‑1 (MQ=255)
                                                                                        cTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGc                                                          >  2:229217/1‑117 (MQ=255)
                                                                                        cTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGCCGGTCCgcccggctggc                                         >  1:1043755/1‑134 (MQ=255)
                                                                                              ctgctgctCGCGCAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTGGCGGGTCGGCCGGTCCgcccggctggcccgcaa                                   <  2:860795/134‑4 (MQ=255)
                                                                                               tgctgctCGCGCAGCGCGCGCTGCTTTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGCCGGTCCgcccggctggcccggctg                                  >  2:929896/1‑134 (MQ=255)
                                                                                                        cgcAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGt                                                                                                <  1:930486/63‑1 (MQ=255)
                                                                                                        cgcAGCGCGCGCTGCTGTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGt                                                                                                >  2:930486/1‑63 (MQ=255)
                                                                                                                      ctgTTCCCGCTGCAGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGCCGGTCCGCCCGGCTGGCCCGGCTGGCCAACAAAGCCGCGGGCTtcg            <  2:390235/133‑1 (MQ=255)
                                                                                                                               ctgcagctCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGCCGGTCCGCCCGGCTGGCCCGGCTGGCCAACAAAGCCGCGGGCTTCGTCGGGCCGgc  <  1:929896/134‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                 
GTTGCTGGTCGCGGGCCTGGTCCTGCCACTGGCGTTGCTGGGCGTCGCGTTGCCGGGCCTCCTGCTGTTGCCGCTGCAGATCAGGCGCCTGCCGCTGCTGCTCGCGCAGCGCGCGCTGCTGTTCCCGCTGCTGCTCCATCTGCTGGCGCTGCAAAACCTGCTGCTGTTCGCGCTGCATGCCCTGCGCCCGCCGTTGCGGGTCGGCCGGTCCGCCCGGCTGGCCCGGCTGGCCAACAAAGCCGCGGGCTTCGTCGGGCCGGC  >  AM260480/2295748‑2296008

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