Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,907,305 2 bp→AT coding (181‑182/1221 nt) h16_A2693 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,907,3050GA100.0% 29.1 / NA 9D61I (GAC→ATC) h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/4);  total (5/4)
*AM2604792,907,3060AT100.0% 28.4 / NA 9D61I (GAC→ATC) h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (5/4);  total (5/4)

CGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGG  >  AM260479/2907205‑2907410
                                                                                                    ||                                                                                                        
gcatgatGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAgg                                                                       >  2:1063896/7‑134 (MQ=255)
                                            ccGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGcc                           <  1:1063896/134‑1 (MQ=255)
                                                         cAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTg               >  1:857340/1‑133 (MQ=255)
                                                                           gTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGcc                                                                     >  1:897847/1‑61 (MQ=255)
                                                                           gTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGcc                                                                     <  2:897847/61‑1 (MQ=255)
                                                                                agctgagctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCt                                                                <  1:147822/61‑1 (MQ=255)
                                                                                agctgagctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCt                                                                >  2:147822/1‑61 (MQ=255)
                                                                                     agctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcgg  <  1:387886/118‑1 (MQ=255)
                                                                                     agctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcgg  >  2:387886/1‑118 (MQ=255)
                                                                                                    ||                                                                                                        
CGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGG  >  AM260479/2907205‑2907410

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