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% 61.4 / NA 17D184E (GAT→GAGh16_A1942Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (8/9);  total (8/9)

CGCGTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCG  >  AM260479/2110879‑2111131
                                                                                                                                      |                                                                                                                        
cgcgTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGAt                                                                                                                           >  1:856903/1‑134 (MQ=255)
         tggGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTcgac                                                                                                                  <  2:1070667/134‑1 (MQ=255)
               cGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACgaggat                                                                                                           <  2:856903/134‑4 (MQ=255)
               cGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACgaggat                                                                                                           >  2:341584/1‑131 (MQ=255)
                     ccAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATc                                                                                                      >  1:1066786/1‑134 (MQ=255)
                                      ccACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCcgccgcgcc                                                                                        <  1:195457/131‑1 (MQ=255)
                                            gcgcAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATggact                                                                                                                     <  1:893366/95‑5 (MQ=255)
                                            gcgcAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATggact                                                                                                                     >  2:893366/1‑91 (MQ=255)
                                                      tCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATg                                                                                                                         <  1:446097/81‑1 (MQ=255)
                                                      tCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATg                                                                                                                         >  2:446097/1‑81 (MQ=255)
                                                        gggCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGc                                                                    <  1:341584/133‑1 (MQ=255)
                                                          gcgcATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTcc                                                                 <  2:766004/134‑1 (MQ=255)
                                                                              ggCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTc                                             >  2:405394/1‑134 (MQ=255)
                                                                                        cGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGAcgccgc                                                                                              >  1:86402/1‑75 (MQ=255)
                                                                                        cGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGAcgccgc                                                                                              <  2:86402/75‑1 (MQ=255)
                                                                                                   ccGACGAGTATGAGACCTTCGACGAACTGGTGGATg                                                                                                                         <  2:227507/36‑1 (MQ=255)
                                                                                                   ccGACGAGTATGAGACCTTCGACGAACTGGTGGATg                                                                                                                         >  1:227507/1‑36 (MQ=255)
                                                                                                                    ttCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTgcg       >  2:202270/1‑134 (MQ=255)
                                                                                                                      cgacgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTgcg       <  1:405394/132‑1 (MQ=255)
                                                                                                                         cgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTg                                                                                    >  1:826875/1‑52 (MQ=37)
                                                                                                                         cgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTg                                                                                    <  2:826875/52‑1 (MQ=37)
                                                                                                                         cgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGcc                                                                              >  1:744247/1‑58 (MQ=37)
                                                                                                                         cgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGcc                                                                              <  2:744247/58‑1 (MQ=37)
                                                                                                                         cgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCgcggcg  >  1:256926/1‑134 (MQ=255)
                                                                                                                                      |                                                                                                                        
CGCGTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCG  >  AM260479/2110879‑2111131

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