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% 35.2 / NA 12S132P (TCC→CCC) ydeNsulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/5);  total (7/5)

TACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAA  >  NZ_CP009273/1576304‑1576477
                                                                                    |                                                                                         
tACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                                      <  2:48322/90‑1 (MQ=255)
     tCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCgtgt                                                                                 <  2:571504/90‑1 (MQ=255)
         gAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCca                                                                             >  1:140612/1‑90 (MQ=255)
                   aTTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGtt                                                                   <  2:592994/90‑1 (MQ=255)
                        aaTAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCCCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTaa                                                              <  2:366634/90‑1 (MQ=255)
                           aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                           >  1:311733/1‑90 (MQ=255)
                           aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                           >  2:652571/1‑90 (MQ=255)
                                                               ttGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCg                       <  1:69130/90‑1 (MQ=255)
                                                                    gcgcGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAt                  >  2:203097/1‑90 (MQ=255)
                                                                      gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt                >  1:171310/1‑90 (MQ=255)
                                                                      gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt                >  1:369578/1‑90 (MQ=255)
                                                                      gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt                >  2:612111/1‑90 (MQ=255)
                                                                                    aaaCACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCaa  <  1:98191/90‑1 (MQ=255)
                                                                                    |                                                                                         
TACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAA  >  NZ_CP009273/1576304‑1576477

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: