Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,616,505 C→G A45G (GCC→GGC)  h16_A2404 → predicted esterase of the alpha‑beta hydrolase superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,616,5050CG95.0% 69.1 / ‑3.0 20coding (134/921 nt)h16_A2404predicted esterase of the alpha‑beta hydrolase superfamily
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (7/12);  minor base . (0/1);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

CGTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCG  >  AM260479/2616381‑2616630
                                                                                                                            |                                                                                                                             
cgTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgc                                                                                                                      >  2:561970/1‑134 (MQ=255)
  tCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgcgg                                                                                                                    <  1:611430/134‑1 (MQ=255)
        ccTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               <  1:79907/133‑1 (MQ=255)
        ccTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               <  2:24945/133‑1 (MQ=255)
         cTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgcgg                                                                                                             <  1:561970/134‑1 (MQ=255)
            ggTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCt                                                                                                           <  2:548133/133‑1 (MQ=255)
                    cgcgGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCata                                                                                                  >  2:795803/1‑134 (MQ=255)
                                gctgGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCCCCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               >  1:593260/1‑109 (MQ=255)
                                gctgGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               <  2:593260/109‑1 (MQ=255)
                                            ctCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGa                                                                            <  1:795803/132‑1 (MQ=255)
                                            ctCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         <  2:5979/134‑1 (MQ=255)
                                             tCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         >  1:307433/1‑134 (MQ=255)
                                               gTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         <  2:307433/132‑1 (MQ=255)
                                                      ccGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgcg                                                                                                              >  1:398097/1‑88 (MQ=255)
                                                      ccGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgcg                                                                                                              <  2:398097/88‑1 (MQ=255)
                                                             ccGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGc                                                          <  1:836543/133‑1 (MQ=255)
                                                             ccGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGc                                                          <  1:8072/133‑1 (MQ=255)
                                                                                                          aGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         >  2:680877/1‑73 (MQ=255)
                                                                                                          aGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         <  1:680877/73‑1 (MQ=255)
                                                                                                                    ggCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCg  >  1:953858/1‑134 (MQ=255)
                                                                                                                            |                                                                                                                             
CGTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCG  >  AM260479/2616381‑2616630

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