Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,576,388 A→G S132P (TCC→CCC)  ydeN ← sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,576,3880AG100.0% 56.9 / NA 17S132P (TCC→CCC) ydeNsulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/7);  total (10/7)

CCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGC  >  NZ_CP009273/1576306‑1576469
                                                                                  |                                                                                 
ccATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCg                                                                            <  1:340588/90‑1 (MQ=255)
     cTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                       <  1:761155/90‑1 (MQ=255)
     atGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                       <  2:597040/89‑1 (MQ=255)
       gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCca                                                                     >  2:192029/1‑90 (MQ=255)
            tCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAg                                                                >  2:303990/1‑90 (MQ=255)
                        tAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                       >  2:478171/1‑71 (MQ=255)
                        tAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                       <  1:478171/71‑1 (MQ=255)
                         aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGGGCCACATAGCCGTTAGTAAAAc                                                   >  1:199249/1‑90 (MQ=255)
                           acacCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGt                                                 >  2:629397/1‑90 (MQ=255)
                                     gggcgggAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCAt                                       >  1:180833/1‑90 (MQ=255)
                                         gggAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAt                                   <  1:303990/90‑1 (MQ=255)
                                                     ggTCATTATTGCGGCGCGGGAGGG‑‑CGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCg                     >  1:500520/1‑90 (MQ=255)
                                                                 ggcgCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa           >  2:345200/1‑90 (MQ=255)
                                                                  gcgcGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAt          >  1:21282/1‑90 (MQ=255)
                                                                    gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt        >  1:302700/1‑90 (MQ=255)
                                                                          ggggCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc  <  2:180833/90‑1 (MQ=255)
                                                                          ggggCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc  <  1:345200/90‑1 (MQ=255)
                                                                                  |                                                                                 
CCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGC  >  NZ_CP009273/1576306‑1576469

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