Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,475,701 T→C S352P (TCG→CCG)  yheS → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,475,7010TC94.1% 47.0 / ‑3.4 17S352P (TCG→CCG) yheSABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (9/7);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

ATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAA  >  NZ_CP009273/3475618‑3475774
                                                                                   |                                                                         
aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                     >  1:174605/1‑90 (MQ=255)
aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGGAAACCGACAt                                                                     >  2:164334/1‑90 (MQ=255)
aTCGCATTATTCTCGACTCGATTAAACGGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                     >  2:405221/1‑90 (MQ=255)
       tattCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                              <  1:345262/90‑1 (MQ=255)
               actcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAg                                                      <  1:20795/90‑1 (MQ=255)
                            gAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTgc                                         >  1:130316/1‑90 (MQ=255)
                                  ggTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGt                                   >  2:229476/1‑90 (MQ=255)
                                   gTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTg                                                <  1:286061/76‑1 (MQ=255)
                                   gTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTg                                                >  2:286061/1‑76 (MQ=255)
                                                   ttGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTc                  <  2:318969/90‑1 (MQ=255)
                                                    tGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAtt                                                                    >  1:220246/1‑39 (MQ=255)
                                                    tGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAtt                                                                    <  2:220246/39‑1 (MQ=38)
                                                           ttAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaaa          <  1:162815/90‑1 (MQ=255)
                                                                  gCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCa   >  2:384255/1‑90 (MQ=255)
                                                                  gCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCa   >  2:270968/1‑90 (MQ=255)
                                                                            gggTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCaa  >  1:33985/1‑81 (MQ=255)
                                                                            gggTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCaa  <  2:33985/81‑1 (MQ=255)
                                                                                   |                                                                         
ATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAA  >  NZ_CP009273/3475618‑3475774

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