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

ACGCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCC  >  AM260479/1071808‑1072052
                                                                                                                               |                                                                                                                       
aCGCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGcc                                                                                                                   >  1:924191/1‑134 (MQ=255)
 cgcgCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGccc                                                                                                                  >  1:41515/1‑134 (MQ=255)
              ccGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCa                                                                                                        <  2:41515/131‑1 (MQ=255)
                          ggggTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAgg                                                                                                              <  2:597144/113‑1 (MQ=255)
                          ggggTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAgg                                                                                                              >  1:597144/1‑113 (MQ=255)
                          ggggTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAAc                                                                                                      >  1:361916/1‑121 (MQ=255)
                          ggggTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAAc                                                                                                      <  2:361916/121‑1 (MQ=255)
                          ggggTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCg                                                                                         >  1:523926/1‑134 (MQ=255)
                                   cGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAg                                                                                 <  2:924191/133‑1 (MQ=255)
                                       gCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGc                                                                            >  1:223991/1‑134 (MQ=255)
                                                         cTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAggg                                                         <  2:223991/134‑2 (MQ=255)
                                                                cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAgggg                                                                                                            <  1:690366/77‑1 (MQ=255)
                                                                cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAgggg                                                                                                            >  2:690366/1‑77 (MQ=255)
                                                                cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAggg                                                                                                             >  2:433714/1‑76 (MQ=255)
                                                                cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAggg                                                                                                             <  1:433714/76‑1 (MQ=255)
                                                                                                    aCCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCt                                                                                              <  1:648695/55‑1 (MQ=39)
                                                                                                    aCCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCt                                                                                              >  2:648695/1‑55 (MQ=39)
                                                                                                                 gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCg                                                                           >  2:232130/1‑61 (MQ=255)
                                                                                                                 gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCg                                                                           <  1:232130/61‑1 (MQ=255)
                                                                                                                 gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATcc  >  2:762153/1‑134 (MQ=255)
                                                                                                                               |                                                                                                                       
ACGCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCC  >  AM260479/1071808‑1072052

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: