Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 477,958 C→A 100% W634C (TGG→TGT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273477,9580CA97.1% 42.0 / ‑6.5 20W634C (TGG→TGTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (9/9);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGTTCTTTTCTTTGGTCAGATAGTAATGCGTTACCTCATTGAGCACTTTCTGTGTACGTTCCTGCGTTGCACCTGCTGGCAGCTGAACCATGGTCATAAACA  >  CP009273/477702‑478148
                                                                                                                                                                                                                                                                |                                                                                                                                                                                              
cGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTACACGATTGCGGGCATGTTAAAGGCGAAAACCATCGCCTCTTTGCTTTGCGAGCAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTTTTTTCTTCGCACGGACGATCGGCCCCATCCTTCCACGAcacg                                                                                                                                                                              >  1:144549/1‑275 (MQ=255)
  tttGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTTCAAGCAACTGGTTACGCGCCTGCGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCACAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGTCGAAAACCCTCGCATATTTGATTTGCGAGAAAGCGCGTGTGGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGCCGATGGGCCAAGTCCTTCAAGGacacata                                                                                                                                                                            >  2:410806/1‑272 (MQ=255)
              ttGGTCAAAATATCAGGGTGCCTGCCTGCGTCTGCAAGCAAATGGTTAAGGGCCTGAGGCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGGCAAAGGCGGTTGCAGTACCCAGTTACACGAGTGCGGGCAGGTTAAAGGCGAAAACGATCGCATCTTGGAGTGGCGAGAAAGCGCGTGTTGCACGAATGGTAATCGCTGCAACTTTGTGTGCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGt                                                                                                                                                                <  2:144543/274‑1 (MQ=255)
                           cAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGATGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTCTTCTGACCACg                                                                                                                                                    >  1:270330/1‑274 (MQ=255)
                                tGCTTCGCTGGTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAGCTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCga                                                                                                                                              >  1:485004/1‑274 (MQ=255)
                                                                                  gTGACCAAGGCCAGCCTGGTCAATCAGCTTAAAGGCAAAGACGGTTGCAGTACCCAGTTACACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCGTTGATGTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACGTTGTTTTCGTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGttcttttctt                                                                                            <  1:23248/275‑1 (MQ=255)
                                                                                                    gTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCa                                                                                                                                                                        >  2:237904/1‑181 (MQ=255)
                                                                                                        aTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATAc                                                                                                                                                                    <  1:237911/181‑1 (MQ=255)
                                                                                                                 aaataaaagCCGGTTGCAGTACTCCGTTCCACGATTGCAGCCAGGTTGAAGGCGAAACCCATCGCATCTTTGATTGGCGAGAAACCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGACCGGACGATCGACACAGTCCTTCAAGGAAACGAACGAAATACCGGTATTCTGACCAAGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGTTCTTTTCTTTGGTCAGATAGTAATGCGTTACATCATGGAg                                                             <  2:423545/270‑1 (MQ=255)
                                                                                                                            ggTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGCCTCCACGAGGTACTTTTCTTAGGTGCGAGGGTACTGCGTTACCGCATTGAGCACtttttttg                                                  >  2:14864/1‑270 (MQ=255)
                                                                                                                                           ttCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACCCAACGTTGTtcttttcttt                                                                                           >  2:250657/1‑219 (MQ=255)
                                                                                                                                               aCGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGTTCTTTTCTTTGGTc                                                                                       <  1:250664/219‑1 (MQ=255)
                                                                                                                                                   ttGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTa                                                                                                                              >  1:177397/1‑176 (MQ=255)
                                                                                                                                                     gCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAaacgaac                                                                                                                                                                           >  1:30219/1‑129 (MQ=255)
                                                                                                                                                     gCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGANACGAACGCAATACCGGTATTCTGACCACGGCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGTTCTTTTCTTTGGTCAGATAGTAATGCGTTACCTCATTGAGCACTTTCTGTGTACGTTCCTGCGTTGCACCTGCTgg                         <  1:410819/275‑1 (MQ=255)
                                                                                                                                                       gggCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACgg                                                                                                                          <  2:177391/176‑1 (MQ=255)
                                                                                                                                                         gCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCaa                                                                                                                                                                       <  2:30218/129‑1 (MQ=255)
                                                                                                                                                                              tCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGCACACCGACTCAACGTTGTTCTTTTCTTTGGTCAGATAGTAATGCGTTACCTCATTGAGCACTTTCTGTGTACGTTCCTGCTTTGCACCTGCTGGCAGCTGAACCATGGTCCTaaaaa  >  2:218053/1‑271 (MQ=255)
                                                                                                                                                                                                                              cTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGTTCTTTTCTTTGGt                                                                                        >  1:177149/1‑139 (MQ=255)
                                                                                                                                                                                                                                  aaCTTTGTTTTCTTCGCCCGGACGATCGGCACAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGTTCTTTTCTTTGGTCAGa                                                                                    <  2:177143/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                              
CGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATCGCATCTTTGATTTGCGAGAAAGCGCGTGTTGCACGCATGGTAATCGCTTCAACTTTGTTTTCTTCGCCCGGACGATCGGCCCAGTCCTTCAAGGAAACGAACGCAATACCGGTATTCTGACCACGTCCCGCAAAGCCGAAGCCGTTAACGGCGAACACCGACTCAACGTTGTTCTTTTCTTTGGTCAGATAGTAATGCGTTACCTCATTGAGCACTTTCTGTGTACGTTCCTGCGTTGCACCTGCTGGCAGCTGAACCATGGTCATAAACA  >  CP009273/477702‑478148

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