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 11S132P (TCC→CCC) ydeNsulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/3);  total (8/3)

TACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAG  >  NZ_CP009273/1576304‑1576471
                                                                                    |                                                                                   
tACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAacac                                                                                <  2:86734/90‑1 (MQ=255)
    aTCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCgtg                                                                            >  1:337953/1‑90 (MQ=255)
                           aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                     >  2:162651/1‑90 (MQ=255)
                           aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAAAACCGTGTGCCACATAGCCGTTAGTAAAAc                                                     >  2:354985/1‑90 (MQ=255)
                                ccAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACg                                                >  2:71625/1‑90 (MQ=255)
                                                       ggTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                         <  1:43507/42‑1 (MQ=255)
                                                       ggTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGc                                                                         >  2:43507/1‑42 (MQ=255)
                                                            ttattGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCg                    >  2:137107/1‑90 (MQ=255)
                                                                  cggcgCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTg              >  1:294630/1‑90 (MQ=255)
                                                                      gcgGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt          >  2:134388/1‑90 (MQ=255)
                                                                              ggCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAg  <  1:29792/90‑1 (MQ=255)
                                                                                    |                                                                                   
TACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAG  >  NZ_CP009273/1576304‑1576471

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