Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,071,996 +G intergenic (+60/‑123) h16_A0979 → / → h16_A0980 Hypothetical protein/Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,071,9941.G100.0% 77.2 / NA 21intergenic (+58/‑125)h16_A0979/h16_A0980Hypothetical protein/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (11/10);  total (11/10)

CGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGCCGCATCCAGCCCGACGTGGTCCCC  >  AM260479/1071872‑1072093
                                                                                                                            |                                                                                                   
cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGt                                                                                             >  2:1119632/1‑133 (MQ=255)
cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGt                                                                                             >  1:874712/1‑133 (MQ=255)
  cTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGAt                                                                                          <  1:798055/134‑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)
                                                  ctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCt                                          >  2:586563/1‑134 (MQ=255)
                                                          cGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTgg                                  <  1:784803/134‑1 (MQ=255)
                                                           gcctggACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGa                                 >  2:838676/6‑134 (MQ=255)
                                                                   gCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGcc                         <  1:863461/134‑1 (MQ=255)
                                                                        aGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGCCGCATc                    <  1:838676/134‑1 (MQ=255)
                                                                         ggggCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGCCGCATc                    >  2:46419/1‑133 (MQ=255)
                                                                                   ggCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGgcggcg                                                                                              >  2:1084479/1‑49 (MQ=38)
                                                                                   ggCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGgcggcg                                                                                              <  1:1084479/49‑1 (MQ=38)
                                                                                            gCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGCCGCATCCAGCCCGACGTGGTcccc  <  1:46419/132‑1 (MQ=255)
                                                                                                 gATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcg                                                                                  <  1:654236/47‑1 (MQ=38)
                                                                                                 gATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcg                                                                                  >  2:654236/1‑47 (MQ=38)
                                                                                                                            |                                                                                                   
CGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGCCGCATCCAGCCCGACGTGGTCCCC  >  AM260479/1071872‑1072093

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: