Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,071,935 +G coding (626/627 nt) h16_A0979 → Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,071,9341.G100.0% 81.6 / NA 22*209G (TGA→GGA) h16_A0979Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (13/9);  total (13/9)

GAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGG  >  AM260479/1071829‑1072061
                                                                                                          |                                                                                                                                
gAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCt                                                                                                       <  1:558831/134‑1 (MQ=255)
     ggggTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCg                                                                                                  >  2:798055/1‑134 (MQ=255)
                 cGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAg                                                                                      >  1:237370/1‑134 (MQ=255)
                     ggCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGt                                                                                  >  1:307800/1‑134 (MQ=255)
                                           cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGt                                                             >  1:874712/1‑133 (MQ=255)
                                           cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGt                                                             >  2:1119632/1‑133 (MQ=255)
                                             cTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGAt                                                          <  1:798055/134‑1 (MQ=255)
                                                tgGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGa                                                                                                                         >  2:589921/1‑68 (MQ=255)
                                                tgGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGa                                                                                                                         <  1:589921/68‑1 (MQ=255)
                                                     gggCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcg                                                  <  2:307800/134‑1 (MQ=255)
                                                            aCCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGAc                                           >  2:175958/1‑134 (MQ=255)
                                                                               aCCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAggg                                                                  >  2:352192/1‑90 (MQ=255)
                                                                               aCCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAggg                                                                  <  1:352192/91‑2 (MQ=255)
                                                                                   cacCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCTGGc                    <  2:237370/134‑1 (MQ=255)
                                                                                   cacCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGc                    <  1:175958/134‑1 (MQ=255)
                                                                                       gTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCgg                 >  2:62068/1‑133 (MQ=255)
                                                                                       gTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCgg                 >  2:784803/1‑133 (MQ=255)
                                                                                          ctgctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAAt             >  2:246108/1‑134 (MQ=255)
                                                                                            gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAg                                                                                      <  2:214792/59‑1 (MQ=39)
                                                                                            gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAg                                                                                      >  1:214792/1‑59 (MQ=39)
                                                                                             ctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCt          >  2:586563/1‑134 (MQ=255)
                                                                                                     cGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTgg  <  1:784803/134‑1 (MQ=255)
                                                                                                          |                                                                                                                                
GAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGG  >  AM260479/1071829‑1072061

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