Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 919,006 T→G C159G (TGC→GGC)  h16_B0812 → hydrolase or acyltransferase (alpha/beta hydrolase superfamily)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260480919,0060TG100.0% 79.5 / NA 22C159G (TGC→GGC) h16_B0812hydrolase or acyltransferase (alpha/beta hydrolase superfamily)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (10/12);  total (10/12)

GCCGGCCTGGCGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGTGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCG  >  AM260480/918887‑919133
                                                                                                                       |                                                                                                                               
gccggccTGGCGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGaa                                                                                                                   <  1:103488/134‑1 (MQ=255)
                               ggTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGcgcc                                                                                                         >  2:905581/1‑113 (MQ=255)
                               ggTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGcgcc                                                                                                         <  1:905581/113‑1 (MQ=255)
                               ggTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     <  2:9082/117‑1 (MQ=255)
                               ggTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     >  1:9082/1‑117 (MQ=255)
                               ggTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGAcc                                                                                     >  2:798619/1‑133 (MQ=255)
                                gTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACcg                                                                                    <  2:903714/133‑1 (MQ=255)
                                         ccgGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     >  1:714987/1‑107 (MQ=255)
                                         ccgGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     <  2:714987/107‑1 (MQ=255)
                                                      ccGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGcc                                                             <  1:608752/134‑1 (MQ=255)
                                                              cgcgcgTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCt                                                     <  2:210823/134‑1 (MQ=255)
                                                                               aCTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGAcc                                                                                           >  1:1038259/1‑79 (MQ=255)
                                                                               aCTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGAcc                                                                                           <  2:1038259/79‑1 (MQ=255)
                                                                                         tATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCgtggt                          <  2:127664/134‑1 (MQ=255)
                                                                                               gATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCgtgg                           <  2:655530/127‑1 (MQ=255)
                                                                                               gATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCgtgg                           >  1:655530/1‑127 (MQ=255)
                                                                                                   ccATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCttt                >  1:26728/1‑134 (MQ=255)
                                                                                                      tCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGaa                                                          >  1:969579/1‑89 (MQ=255)
                                                                                                      tCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGaa                                                          <  2:969579/89‑1 (MQ=255)
                                                                                                      tCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCGGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAg             >  2:318035/1‑134 (MQ=255)
                                                                                                                  gcCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATcg  <  2:384951/133‑1 (MQ=255)
                                                                                                                  gcCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATcg  >  1:384951/1‑133 (MQ=255)
                                                                                                                       |                                                                                                                               
GCCGGCCTGGCGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGTGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCG  >  AM260480/918887‑919133

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: