Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 500,156 A→G intergenic (‑4/‑214) fsr ← / → ushA fosmidomycin MFS transporter/bifunctional UDP‑sugar hydrolase/5'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273500,1560AG87.5% 39.5 / ‑0.7 16intergenic (‑4/‑214)fsr/ushAfosmidomycin MFS transporter/bifunctional UDP‑sugar hydrolase/5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (8/6);  total (9/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 = 6.40e-01

ACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAA  >  NZ_CP009273/500069‑500241
                                                                                       |                                                                                     
aCCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGgaa                                                                                     >  2:324816/1‑90 (MQ=255)
           cAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTc                                                                          >  1:535215/1‑90 (MQ=255)
                     gcGGGCTTTGGTCGTCGCCGCGGCCGCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGa                                                                <  2:535215/89‑1 (MQ=255)
                           cTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTaaa                                                          >  2:170904/1‑90 (MQ=255)
                           cTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTaaa                                                          >  2:133992/1‑90 (MQ=255)
                              ggTCGTTGCCGCCCCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGaa                                                       <  1:932079/90‑1 (MQ=255)
                                        gCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAtt                                             >  1:149642/1‑90 (MQ=255)
                                                        caGGCTGGGGTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGa                             <  1:543910/90‑1 (MQ=255)
                                                        caGGCTGGGGTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGa                             <  1:483996/90‑1 (MQ=255)
                                                          ggCTGGGGTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACt                           <  1:393434/90‑1 (MQ=255)
                                                               gggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACAt                      <  1:873486/90‑1 (MQ=255)
                                                                ggTTGTTCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATa                     <  1:133992/90‑1 (MQ=255)
                                                                   tgttCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCt                  >  1:441567/1‑90 (MQ=255)
                                                                    gttCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTg                 >  1:620371/1‑90 (MQ=255)
                                                                    gttCACTCATTGCCATAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTg                 >  1:781437/1‑90 (MQ=255)
                                                                    gttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTg                 >  1:342563/1‑90 (MQ=255)
                                                                                   tAGGGAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTgaa  >  2:44411/1‑90 (MQ=255)
                                                                                       |                                                                                     
ACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAA  >  NZ_CP009273/500069‑500241

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