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% 70.4 / NA 19*209G (TGA→GGA) h16_A0979Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (10/9);  total (10/9)

GTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCC  >  AM260479/1071840‑1072052
                                                                                               |                                                                                                                       
gTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTggacg                                                                                                                    >  1:416960/1‑96 (MQ=255)
gTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTggacg                                                                                                                    <  2:416960/100‑5 (MQ=255)
                         cTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAggg                                                         <  2:48937/134‑2 (MQ=255)
                                cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGt                                                    >  2:530721/1‑133 (MQ=255)
                                     tgGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGa                                              >  2:950297/1‑134 (MQ=255)
                                     tgGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGa                                              >  2:650381/1‑134 (MQ=255)
                                                              tGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGt                                                                         >  1:824391/1‑82 (MQ=255)
                                                              tGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGt                                                                         <  2:824391/82‑1 (MQ=255)
                                                                cAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGgccg                    <  1:650381/133‑1 (MQ=255)
                                                                         acCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcg                                         <  2:430773/103‑1 (MQ=255)
                                                                         acCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcg                                         >  1:430773/1‑103 (MQ=255)
                                                                            gTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAg                                                           >  1:570134/1‑81 (MQ=255)
                                                                            gTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAg                                                           <  2:570134/81‑1 (MQ=255)
                                                                               ctgctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGaaa     <  1:530721/133‑1 (MQ=255)
                                                                                 gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGc                          >  2:402760/1‑110 (MQ=255)
                                                                                 gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGc                          <  1:402760/110‑1 (MQ=255)
                                                                                 gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATcc  >  1:632621/1‑134 (MQ=255)
                                                                                   tgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAggg                                                         <  2:1146602/76‑2 (MQ=255)
                                                                                   tgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAggg                                                         >  1:1146602/1‑75 (MQ=255)
                                                                                         aCGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTg                                                                        >  2:254043/1‑56 (MQ=39)
                                                                                         aCGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTg                                                                        <  1:254043/56‑1 (MQ=255)
                                                                                               |                                                                                                                       
GTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCC  >  AM260479/1071840‑1072052

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: