Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 978,518 A→G S117G (AGC→GGC)  h16_B0862 → cytosine deaminase or related metal‑dependent hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260480978,5180AG100.0% 146.0 / NA 38S117G (AGC→GGC) h16_B0862cytosine deaminase or related metal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/19);  total (19/19)

GCTGCCGCGCATGCTGCTGATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAAGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGGCACGCGCAGCGACGCCGGGC  >  AM260480/978391‑978629
                                                                                                                               |                                                                                                               
gctgccgcGCATGCTGCTGATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGcc                                                                                                           >  2:1063907/1‑134 (MQ=255)
      gcgcATGCTGCTGATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAgg                                                                                                     <  1:1063907/134‑1 (MQ=255)
      gcgcATGCTGCTGATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAgg                                                                                                     <  1:362025/134‑1 (MQ=255)
        gcATGCTGCTGATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCg                                                                                                   <  1:435186/134‑1 (MQ=255)
                   aTGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCacctgta                                                                                        <  1:337139/134‑1 (MQ=255)
                              gATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTcc                                                                              >  1:842799/1‑133 (MQ=255)
                                             ccaccTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCacctgta                                                                                        >  1:233133/1‑108 (MQ=255)
                                             ccaccTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCacctgta                                                                                        <  2:233133/108‑1 (MQ=255)
                                                    aCGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCt                                                       <  1:657049/134‑1 (MQ=255)
                                                     cGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTg                                                      >  2:728831/1‑134 (MQ=255)
                                                       cAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGc                                                                          <  2:783742/112‑1 (MQ=255)
                                                       cAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGc                                                                          >  1:783742/1‑112 (MQ=255)
                                                       cAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTgcg                                                    >  2:455596/1‑134 (MQ=255)
                                                         gTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGcc                                                                                                           >  2:443959/1‑77 (MQ=255)
                                                         gTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGcc                                                                                                           <  1:443959/77‑1 (MQ=255)
                                                          tCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCgcgg                                                 <  2:842799/134‑1 (MQ=255)
                                                             ttCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGcggcgg                                              >  2:1037528/1‑134 (MQ=255)
                                                             ttCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGcggcgg                                              >  1:489479/1‑134 (MQ=255)
                                                              tCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTgcg                                                    >  2:36674/1‑127 (MQ=255)
                                                              tCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTgcg                                                    <  1:36674/127‑1 (MQ=255)
                                                              tCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCggcggc                                             >  2:866695/1‑134 (MQ=255)
                                                               cGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAg                                                                                                      <  1:164473/76‑1 (MQ=255)
                                                               cGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAg                                                                                                      >  2:164473/1‑76 (MQ=255)
                                                                              gtTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAg                                                                                                      >  1:1071158/1‑61 (MQ=255)
                                                                              gtTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAg                                                                                                      <  2:1071158/61‑1 (MQ=255)
                                                                               ttCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGa                            <  2:489479/134‑1 (MQ=255)
                                                                                         cATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGGCAcg                  <  2:922873/134‑1 (MQ=255)
                                                                                                  tCTTCCGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGGCACGCGCAGCgca         >  1:207227/1‑132 (MQ=255)
                                                                                                      ccGGCGCGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGGCACGCGCAGCGACGc       <  1:1037528/132‑1 (MQ=255)
                                                                                                         gcgcGTGTGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGGCACGCGCAGCGACGCCGGGc  >  1:640433/1‑134 (MQ=255)
                                                                                                             gtgtGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGGCACGCGCAGCg           <  2:465641/121‑1 (MQ=255)
                                                                                                             gtgtGGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGGCACGCGCAGCg           >  1:465641/1‑121 (MQ=255)
                                                                                                               gtgGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCg                                   <  2:305424/95‑1 (MQ=255)
                                                                                                               gtgGGTGCCGCTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCg                                   >  1:305424/1‑95 (MQ=255)
                                                                                                                         cTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGgcggcg                                                                  >  1:759450/1‑54 (MQ=255)
                                                                                                                         cTGGAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGgcggcg                                                                  <  2:759450/54‑1 (MQ=255)
                                                                                                                           ggAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTgcg                                                    >  2:302536/1‑66 (MQ=255)
                                                                                                                           ggAAGGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTgcg                                                    <  1:302536/66‑1 (MQ=255)
                                                                                                                               |                                                                                                               
GCTGCCGCGCATGCTGCTGATGCCCGGCATGATCGACACGCACCACCACCTGACGCAGTCGTTCGGCAAGTCGCTGGTGTTCGGCGAACCATCGGAGATCTTCCGGCGCGTGTGGGTGCCGCTGGAAAGCAGCCTGAAGGCGGAGCACCTGTACCTGTCGTCCAAGCTGGCGGCGCTGGAAGCGCTGCGCGGCGGCTTCACCACCGTGGTCGATGCCGGCACGCGCAGCGACGCCGGGC  >  AM260480/978391‑978629

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: