Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 582,759 G→T C352F (TGT→TTT)  h16_A0555 → zinc metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479582,7590GT100.0% 120.4 / NA 32C352F (TGT→TTT) h16_A0555zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/16);  total (16/16)

GGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCC  >  AM260479/582628‑582884
                                                                                                                                   |                                                                                                                             
ggggACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGtttg                                                                                                                             >  2:549377/1‑134 (MQ=255)
    aCCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGcgcg                                                                                                                         <  1:176712/134‑1 (MQ=255)
        cgGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATg                                                                                                                     >  2:248027/1‑134 (MQ=255)
           aaaaTTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGcgcg                                                                                                                         >  1:448/1‑127 (MQ=255)
           aaaaTTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGcgcg                                                                                                                         <  2:448/127‑1 (MQ=255)
               ttCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGccc                                                                                                              <  2:366978/134‑1 (MQ=255)
                                              cTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGccc                                                                                                              <  1:538972/103‑1 (MQ=255)
                                              cTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATACGGGccc                                                                                                              >  2:538972/1‑103 (MQ=255)
                                                   cTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAg                                                                           >  2:951584/1‑133 (MQ=255)
                                                   cTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCAAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAg                                                                           >  2:981176/1‑133 (MQ=255)
                                                       gAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGccc                                                                                       <  1:198251/117‑1 (MQ=255)
                                                       gAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGccc                                                                                       >  2:198251/1‑117 (MQ=255)
                                                        aGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTg                                                                     <  1:951584/134‑1 (MQ=255)
                                                        aGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTg                                                                     <  1:981176/134‑1 (MQ=255)
                                                            gATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGcc                                                                                                               >  1:665550/1‑88 (MQ=255)
                                                            gATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGcc                                                                                                               <  2:665550/88‑1 (MQ=255)
                                                             aTTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTgcgcgc                                                                <  2:361100/134‑1 (MQ=255)
                                                                 ccGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCgcggcg                                                            >  1:435458/1‑134 (MQ=255)
                                                                    aCATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTg                                                         >  1:347520/1‑134 (MQ=255)
                                                                           ggACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGcgcacg                                                   >  1:593037/1‑133 (MQ=255)
                                                                               cgcgTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGccc                                                                                       >  1:563048/1‑93 (MQ=255)
                                                                               cgcgTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGccc                                                                                       <  2:563048/93‑1 (MQ=255)
                                                                                    cAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACt                                              >  1:972039/1‑129 (MQ=255)
                                                                                    cAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACt                                              <  2:972039/129‑1 (MQ=255)
                                                                                     aGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCgcgg                                                              >  2:697640/1‑112 (MQ=255)
                                                                                     aGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCgcgg                                                              <  1:697640/112‑1 (MQ=255)
                                                                                              aGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACgg                                >  1:828656/1‑133 (MQ=255)
                                                                                                cAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCAc                             <  1:549377/134‑1 (MQ=255)
                                                                                                     gcAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTgcgc                        >  2:628772/1‑134 (MQ=255)
                                                                                                                       cTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTAc      <  1:628772/134‑1 (MQ=255)
                                                                                                                           aTACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGgcc  <  2:593037/134‑1 (MQ=255)
                                                                                                                              cGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCa                                                                                      <  2:1061237/47‑1 (MQ=255)
                                                                                                                              cGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCa                                                                                      >  1:1061237/1‑47 (MQ=255)
                                                                                                                                   |                                                                                                                             
GGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCC  >  AM260479/582628‑582884

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