Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,464,365 Δ1 bp coding (49/552 nt) h16_A1354 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,464,3650A.100.0% 96.0 / NA 22coding (49/552 nt)h16_A1354conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (11/11);  total (11/11)

TGGGCGCGGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTG  >  AM260479/1464245‑1464491
                                                                                                                        |                                                                                                                              
tGGGCGCGGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGa                                                                                                                  <  2:668762/134‑1 (MQ=255)
    cgcgGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCAc                                                                                                              <  2:134650/134‑1 (MQ=255)
                       caATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTg                                                                                           >  1:595460/1‑134 (MQ=255)
                            gCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGc                                                                                       >  2:477985/1‑133 (MQ=255)
                              tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTgcc                                                                                                                           >  1:1021979/1‑95 (MQ=255)
                              tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTgcc                                                                                                                           <  2:1021979/95‑1 (MQ=255)
                              tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCAt                                                                                    >  2:357321/1‑134 (MQ=255)
                                   gcgcCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGGACCGTCCCGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGc                                                                                >  2:546828/1‑133 (MQ=255)
                                        gccgcGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCgg                                                                          <  2:595460/134‑1 (MQ=255)
                                          cgcgCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTg                                                                        >  1:388757/1‑134 (MQ=255)
                                                               ggAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCca                                                                                                                       <  1:300513/66‑1 (MQ=255)
                                                               ggAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCca                                                                                                                       >  2:300513/1‑66 (MQ=255)
                                                                  gaTGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctg                                                <  2:285393‑M1/134‑6 (MQ=255)
                                                                     gCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACagccgct                                                       >  2:677657/1‑118 (MQ=255)
                                                                     gCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACagccgct                                                       <  1:677657/124‑7 (MQ=255)
                                                                      ccATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacg                                            >  1:525485‑M1/1‑125 (MQ=255)
                                                                      ccATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacg                                            >  1:512028‑M1/1‑125 (MQ=255)
                                                                                   tgACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcc                                <  1:256727‑M1/133‑22 (MQ=255)
                                                                                       cTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                           <  2:388757‑M1/134‑27 (MQ=255)
                                                                                                     ggtggGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacgg             >  1:506153‑M1/1‑94 (MQ=255)
                                                                                                             aaCCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggt     <  1:477985‑M1/134‑49 (MQ=255)
                                                                                                                cGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgag  <  1:357321‑M1/134‑52 (MQ=255)
                                                                                                                        |                                                                                                                              
TGGGCGCGGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTG  >  AM260479/1464245‑1464491

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