Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,575,730 T→C D351G (GAT→GGT)  ydeN ← sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,575,7300TC100.0% 102.4 / NA 29D351G (GAT→GGT) ydeNsulfatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/13);  total (16/13)

ATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGC  >  NZ_CP009273/1575595‑1575860
                                                                                                                                       |                                                                                                                                  
aTTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGa                                                                                                                                 <  1:329467/139‑1 (MQ=255)
        aaTCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCa                                                                                                                         <  1:447201/139‑1 (MQ=255)
                          aCCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTaaa                                                                                                        <  2:341966/138‑1 (MQ=255)
                            cGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTaaa                                                                                                        >  1:145149/1‑136 (MQ=255)
                            cGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTaaa                                                                                                        <  2:145149/136‑1 (MQ=255)
                                    aGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattg                                                                                             <  1:159961/139‑1 (MQ=255)
                                                                gggTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATCGTCATACTGTCCGTTTttcttc                                                                   <  2:481664/137‑1 (MQ=255)
                                                                gggTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATCGTCATACTGTCCGTTTttcttc                                                                   >  1:481664/1‑137 (MQ=255)
                                                                     tgAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAgttgtt                                                            >  1:380787/1‑139 (MQ=255)
                                                                            cGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAAt                                                     >  2:26309/1‑139 (MQ=255)
                                                                                 ggATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCgaga                                                       <  1:407856/132‑1 (MQ=255)
                                                                                 ggATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCgaga                                                       >  2:407856/1‑132 (MQ=255)
                                                                                   aTAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGttt                                               >  2:169701/1‑138 (MQ=255)
                                                                                       gTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTacac                                          >  2:252447/1‑139 (MQ=255)
                                                                                       gTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTacac                                          >  1:512792/1‑139 (MQ=255)
                                                                                                  ataGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCaa                                 <  2:380787/137‑1 (MQ=255)
                                                                                                          tttGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattg                                                                                             >  1:77745/1‑69 (MQ=255)
                                                                                                          tttGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAattg                                                                                             <  2:77745/69‑1 (MQ=255)
                                                                                                          tttGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAAc                       >  2:324319/1‑139 (MQ=255)
                                                                                                          tttGCGCCCCGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAAc                       >  1:234462/1‑139 (MQ=255)
                                                                                                                 cccGTTCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCg                <  2:512792/139‑1 (MQ=255)
                                                                                                                     ttCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAgaga                                                                                                    <  1:26515/51‑1 (MQ=255)
                                                                                                                     ttCAGCGGCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAgaga                                                                                                    >  2:26515/1‑51 (MQ=255)
                                                                                                                            gCAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATc     >  2:262687/1‑139 (MQ=255)
                                                                                                                             cAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCgtagta           <  2:251575/132‑1 (MQ=255)
                                                                                                                             cAGAGGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCgtagta           >  1:251575/1‑132 (MQ=255)
                                                                                                                               gagGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACgg                     >  1:75815/1‑120 (MQ=255)
                                                                                                                               gagGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACgg                     <  2:75815/120‑1 (MQ=255)
                                                                                                                               gagGACCACCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATctgc  >  1:249336/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                  
ATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCCCGTTCAGCGGCAGAGGACCATCGATAACCGCACCATTATCGGAGGTAAAGAGAATAATTGTATTGTCATACTGTCCGTTTTTCTTCAGTTGTTCGAGAATGCGTTTTACACCCTGATCAACAGAATAAACGGAAGCGTAGTAGTTATCTGC  >  NZ_CP009273/1575595‑1575860

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