Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 1,518,085 C→T A85T (GCC→ACC)  h16_B1356 ← conserved hypothetical membrane spanning protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604801,518,0850CT100.0% 89.1 / NA 25A85T (GCC→ACC) h16_B1356conserved hypothetical membrane spanning protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/13);  total (12/13)

GATGGCGGGGGCCGCCGCGGCGGCCATGCTGGCAAGCCCGAGCAGGGCCGCGCACAGGTAACGGTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGCCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGTCATCACACCCGCGCCGAGCGGTGTGCCGAGCGTGGCCGGCGAGAA  >  AM260480/1517955‑1518202
                                                                                                                                  |                                                                                                                     
gATGGCGGGGGCCGCCGGGGCGGCCATGCTGGCAAGCCCGAGCAGGGCCGCGCACAGGTAACGGTTGGGGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGtcc                                                                                                                     <  2:1086895/133‑1 (MQ=255)
    gCGGGGGCCGCCGCGGCGGCCATGCTGGCAAGCCCGAGCAGGGCCGCGCACAGGTAACGGTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTac                                                                                                                <  1:341997/134‑1 (MQ=255)
         ggCCGCCGCGGCGGCCATGCTGGCAAGCCCGAGCAGGGCCGCGCACAGGTAACGGTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACg                                                                                                             >  2:418832/1‑132 (MQ=255)
           ccgccgCGGCGGCCATGCTGGCAAGCCCGAGCAGGGCCGCGCACAGGTAACGGTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACg                                                                                                             >  1:681863/1‑130 (MQ=255)
                   gcggcCATGCTGGCAAGCCCGAGCAGGGCCGCGCACAGGTAACGGTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGt                                                                                                  >  2:945009/1‑133 (MQ=255)
                    cggcCATGCTGGCAAGCCCGAGCAGGGCCGCGCACAGGTAACGGTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTca                                                                                                <  2:264673/134‑1 (MQ=255)
                            cTGGCAAGCCCGAGCTGGGCCGCGCACAGGTAACGGTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGCGCGGGTCCAGGTCGCGGTAGGTCCGGTACACGTAGATCCCGGTCAGCAGCAc                                                                                         <  2:72871/133‑1 (MQ=255)
                                   gCCCGAGCAGGGCCGCGCACAGGTAACGGTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCgaga                                                                                 <  2:681863/134‑1 (MQ=255)
                                                               gTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGcc                                                     <  1:945009/134‑1 (MQ=255)
                                                                     gggTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTca                                                                                                <  2:537981/85‑1 (MQ=255)
                                                                     gggTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTca                                                                                                >  1:537981/1‑85 (MQ=255)
                                                                       gTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGtca                                             >  1:760587/1‑134 (MQ=255)
                                                                       gTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGtca                                             >  2:539832/1‑134 (MQ=255)
                                                                       gTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGtca                                             >  2:1144329/1‑134 (MQ=255)
                                                                         ttCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGtcat                                            >  2:1026500/1‑133 (MQ=255)
                                                                             tGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGTCATCacac                                       >  2:330233/1‑134 (MQ=255)
                                                                                  ggCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGTCATCACACCcgcg                                  <  1:418832/134‑1 (MQ=255)
                                                                                     aCTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTcagc                                                                                              <  2:757001/71‑1 (MQ=255)
                                                                                     aCTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTcagc                                                                                              >  1:757001/1‑71 (MQ=255)
                                                                                     aCTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGTCATCACACCCGCGCCg                               <  1:539832/134‑1 (MQ=255)
                                                                                               cACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGTCATCACACCCGCGCCGAGCGGTGTGc                     <  1:1144329/134‑1 (MQ=255)
                                                                                                       cgttGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGa                                                               >  2:859682/1‑84 (MQ=255)
                                                                                                       cgttGAGCGGATCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGa                                                               <  1:859682/84‑1 (MQ=255)
                                                                                                                  tCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGcc                                                     <  2:1145151/83‑1 (MQ=255)
                                                                                                                  tCCAGGTCGCGGTTGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGcc                                                     >  1:1145151/1‑83 (MQ=255)
                                                                                                                               tGGTCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGTCATCACACCCGCGCCGAGCGGTGTGCCGAGCGTGGCCGGCGAGaa  <  1:1029425/121‑1 (MQ=255)
                                                                                                                                  |                                                                                                                     
GATGGCGGGGGCCGCCGCGGCGGCCATGCTGGCAAGCCCGAGCAGGGCCGCGCACAGGTAACGGTTGGCGGGTTTCATGCCCGGCACTCCTTGCGCACGCGTTCGTTGAGCGGATCCAGGTCGCGGTTGGCCCGGTACACGTAGATCCCGGTCAGCAGCACGGCCGAGACGATCACCGCCACGCCGACGGGAATGCCCACCGTCATCACACCCGCGCCGAGCGGTGTGCCGAGCGTGGCCGGCGAGAA  >  AM260480/1517955‑1518202

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