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% 81.7 / 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)

CGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGTGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCG  >  AM260480/918897‑919133
                                                                                                             |                                                                                                                               
cggcctggcCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGcgcc                                                                                                         <  2:623537/134‑1 (MQ=255)
        ccAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGc                                                                                                 >  1:708108/1‑134 (MQ=255)
                      gTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACcg                                                                                    <  2:1129961/133‑1 (MQ=255)
                          aCCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGAc                                                                                                                         <  2:1028328/92‑1 (MQ=255)
                          aCCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGAc                                                                                                                         >  1:1028328/1‑92 (MQ=255)
                          aCCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     <  2:759269/112‑1 (MQ=255)
                          aCCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     >  1:759269/1‑112 (MQ=255)
                               ccgGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGct                                                                          >  1:925063/1‑134 (MQ=255)
                               ccgGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGct                                                                          >  2:377560/1‑134 (MQ=255)
                                              gggCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGTCCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGa                                                           <  1:391004/134‑1 (MQ=255)
                                                 cgccgcGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATg                                                        >  1:607232/1‑134 (MQ=255)
                                                  gccgcGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATgg                                                       <  1:817894/134‑1 (MQ=255)
                                                                         ccGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGAc                                                                                                                         >  1:586183/1‑45 (MQ=255)
                                                                         ccGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGAc                                                                                                                         <  2:586183/45‑1 (MQ=255)
                                                                          cGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     <  2:1117902/64‑1 (MQ=255)
                                                                          cGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     >  1:1117902/1‑64 (MQ=255)
                                                                              gTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCgtgg                           <  2:708108/134‑1 (MQ=255)
                                                                                              cGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAgg            <  1:377560/133‑1 (MQ=255)
                                                                                               gCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACcg                                                                                    <  1:712620/60‑1 (MQ=255)
                                                                                               gCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACcg                                                                                    >  2:712620/1‑60 (MQ=255)
                                                                                                ccTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGAtt         >  1:336415/1‑134 (MQ=255)
                                                                                                       cgcCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATcg  <  1:971622/134‑1 (MQ=255)
                                                                                                             |                                                                                                                               
CGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGTGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCG  >  AM260480/918897‑919133

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