Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,111,013 +G coding (553/1089 nt) h16_A1942 → Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,111,0121.G100.0% 53.2 / NA 15D184E (GAT→GAGh16_A1942Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (8/7);  total (8/7)

GCGTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGC  >  AM260479/2110880‑2111127
                                                                                                                                     |                                                                                                                    
gcgTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATg                                                                                                                     <  1:1068206/134‑1 (MQ=255)
gcgTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGAGGATg                                                                                                                     <  1:28226/134‑1 (MQ=255)
                   gCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGAt                                                                                                   >  2:421723/1‑134 (MQ=255)
                    ccAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATc                                                                                                  >  2:947486/1‑134 (MQ=255)
                        cATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGAcg                                                                                              <  1:384048/134‑1 (MQ=255)
                                            cgcAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGcc                                                                          >  1:72436/1‑134 (MQ=255)
                                                             aTCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTccgccg                                                         >  1:399854/1‑134 (MQ=255)
                                                                    aTGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGccctc                                                  <  2:72436/134‑1 (MQ=255)
                                                                                         ccGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGcgtgccgcg                             <  2:134515/134‑1 (MQ=255)
                                                                                            aCGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGcgccc                          >  2:495080/1‑134 (MQ=255)
                                                                                                                     cgacgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTgcgc  <  1:16050/133‑1 (MQ=255)
                                                                                                                      gacgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTg                                                                             <  1:486215/57‑1 (MQ=37)
                                                                                                                      gacgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTg                                                                             >  2:486215/1‑57 (MQ=37)
                                                                                                                       acgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTg     <  1:478119/128‑1 (MQ=255)
                                                                                                                       acgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTg     >  2:478119/1‑128 (MQ=255)
                                                                                                                        cgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTcgcg             >  2:317145/1‑119 (MQ=255)
                                                                                                                        cgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTcgcg             <  1:317145/119‑1 (MQ=255)
                                                                                                                                     |                                                                                                                    
GCGTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGC  >  AM260479/2110880‑2111127

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