Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,602,688 (G)6→5 intergenic (‑85/+41) h16_A3328 ← / ← h16_A3329 conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,602,6830G.100.0% 100.4 / NA 23intergenic (‑80/+46)h16_A3328/h16_A3329conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (13/10);  total (13/10)

GTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTG  >  AM260479/3602557‑3602808
                                                                                                                               |                                                                                                                             
gTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGGGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCgggggc                                                                                                                          >  2:33316/1‑132 (MQ=255)
   ggCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATg                                                                                                                     <  2:627765/134‑1 (MQ=255)
     cAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGGGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGc                                                                                                                   <  1:33316/134‑1 (MQ=255)
       ggTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGt                                                                                                                 >  1:83352/1‑134 (MQ=255)
                                     cgcCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg                                                                                   <  2:83352/134‑1 (MQ=255)
                                                            cgcAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCa                                                                                               >  2:68624/1‑99 (MQ=255)
                                                            cgcAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCa                                                                                               <  1:68624/99‑1 (MQ=255)
                                                                    ccGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGcgccgc                                                    >  1:546790/1‑134 (MQ=255)
                                                                     cGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCcgcg                                                   >  1:379633/1‑134 (MQ=255)
                                                                                 tgcgaaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGctc                                       >  1:66733/6‑134 (MQ=255)
                                                                                    gaaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                                    <  2:379633/132‑1 (MQ=255)
                                                                                     aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCgg                                                                                                             >  2:7338/2‑60 (MQ=255)
                                                                                     aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCgg                                                                                                             <  1:7338/59‑1 (MQ=255)
                                                                                     aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCtctc                                     >  2:12340/2‑132 (MQ=255)
                                                                                      aTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGc                                                                                                               <  2:440197/57‑1 (MQ=255)
                                                                                      aTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGc                                                                                                               >  1:440197/1‑57 (MQ=255)
                                                                                      aTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAg                                  >  2:962455/1‑134 (MQ=255)
                                                                                          aGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCgggggc                                                                                                                          >  1:740727/1‑42 (MQ=255)
                                                                                          aGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCgggggc                                                                                                                          <  2:740727/43‑2 (MQ=255)
                                                                                                       gATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                                    <  1:817126/115‑1 (MQ=255)
                                                                                                       gATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                                    >  2:817126/1‑115 (MQ=255)
                                                                                                          cGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGtgtg                                                                                                  >  1:801850/1‑50 (MQ=255)
                                                                                                          cGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGtgtg                                                                                                  <  2:801850/50‑1 (MQ=255)
                                                                                                               gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg                                                                                   <  2:452228/60‑1 (MQ=255)
                                                                                                               gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg                                                                                   >  1:452228/1‑60 (MQ=255)
                                                                                                                      gcatgcaGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTg  >  1:505775/1‑134 (MQ=255)
                                                                                                                               |                                                                                                                             
GTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTG  >  AM260479/3602557‑3602808

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: