Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,602,642 +A intergenic (‑39/+87) h16_A3328 ← / ← h16_A3329 conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,602,6411.A100.0% 72.5 / NA 20intergenic (‑38/+88)h16_A3328/h16_A3329conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (9/11);  total (9/11)

GGTCCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCC  >  AM260479/3602510‑3602767
                                                                                                                                    |                                                                                                                              
ggTCCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaa                                                                                                                              >  1:583974/1‑133 (MQ=255)
    cATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                          <  2:1185898/134‑6 (MQ=255)
    cATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                          <  2:1153847/134‑6 (MQ=255)
               tGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCgg                                                                                                                >  1:605289/1‑134 (MQ=255)
                          cccATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTgcgc                                                                                                     <  2:583974/134‑1 (MQ=255)
                             aTCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTgcg                                                                                                      >  1:1057154/1‑130 (MQ=255)
                             aTCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTgcg                                                                                                      <  2:1057154/130‑1 (MQ=255)
                                  caccacCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                          >  1:132435/1‑99 (MQ=255)
                                  caccacCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                          <  2:132435/104‑6 (MQ=255)
                                    ccaccaTCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                          <  2:940363/102‑6 (MQ=255)
                                    ccaccaTCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                          <  2:894249/102‑6 (MQ=255)
                                    ccaccaTCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                          >  1:894249/1‑97 (MQ=255)
                                    ccaccaTCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                          >  1:940363/1‑97 (MQ=255)
                                              ggTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCgggggc                                                                                 >  1:956286/1‑133 (MQ=255)
                                                      ggTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGt                                                                        >  2:527733/1‑134 (MQ=255)
                                                        tCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGc                                                                      <  2:605289/134‑1 (MQ=255)
                                                           cTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGc                                                                      <  2:683925/131‑1 (MQ=255)
                                                           cTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGc                                                                      >  1:683925/1‑131 (MQ=255)
                                                                     gcgcGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAt                                                                                                 >  1:186589/1‑95 (MQ=255)
                                                                     gcgcGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAt                                                                                                 <  2:186589/95‑1 (MQ=255)
                                                                          aTCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCg                                                                                                                                 >  1:25979/1‑58 (MQ=255)
                                                                          aTCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCg                                                                                                                                 <  2:25979/58‑1 (MQ=255)
                                                                          aTCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGt                                                                                                         <  1:130248/82‑1 (MQ=255)
                                                                          aTCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGt                                                                                                         >  2:130248/1‑82 (MQ=255)
                                                                                    cgcCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg                                          <  2:106979/134‑1 (MQ=255)
                                                                                    cgcCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg                                          <  1:881201/134‑1 (MQ=255)
                                                                                    cgcCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg                                          <  1:527733/134‑1 (MQ=255)
                                                                                             ccccGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAt                                                                                                 <  2:458729/71‑1 (MQ=255)
                                                                                             ccccGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAt                                                                                                 >  1:458729/1‑71 (MQ=255)
                                                                                                                        aGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGcc      >  1:810743/1‑134 (MQ=255)
                                                                                                                             tAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTcc  <  2:956286/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                                              
GGTCCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCC  >  AM260479/3602510‑3602767

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