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,3880AG85.7% 10.5 / ‑2.7 7S132P (TCC→CCC) ydeNsulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (4/2);  total (4/3)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

CTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAT  >  NZ_CP009273/1576311‑1576436
                                                                             |                                                
cTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGc                                      <  1:100776/90‑1 (MQ=255)
       tCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAg                               >  1:216126/1‑90 (MQ=255)
                    aGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAAc                  >  1:253702/1‑90 (MQ=255)
                          cAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCgtgt                                        <  1:353855/62‑1 (MQ=255)
                          cAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCgtgt                                        >  2:353855/1‑62 (MQ=255)
                                   cgggAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCa   >  2:105732/1‑90 (MQ=255)
                                    gggAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGGAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAt  <  2:253702/90‑1 (MQ=255)
                                                                             |                                                
CTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAT  >  NZ_CP009273/1576311‑1576436

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