Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,010,363 C→T C146Y (TGT→TAT)  USA300HOU_RS10125 ← Putative multidrug export ATP‑binding/permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,010,3630CT96.5% 84.8 / NA 29C146Y (TGT→TAT) USA300HOU_RS10125Putative multidrug export ATP‑binding/permease protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/18);  minor base A (1/0);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTT  >  USA300TCH1516_ALE/2010230‑2010496
                                                                                                                                     |                                                                                                                                     
aCGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAATCGGGGCTACTAAAGACATTATCACATCTAAATAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTaa                                                                                                                                >  1:143690/1‑141 (MQ=255)
  gtgtCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAcc                                                                                                                              >  2:289157/1‑141 (MQ=255)
                aaTCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGCAGATAAACAGTGCAGCTCAAGTCAATTTCACATCAAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAc                                                                                                                <  1:39660/141‑1 (MQ=255)
                   cTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCg                                                                                                             <  1:129672/141‑1 (MQ=255)
                     tCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGt                                                                                                           <  1:261395/141‑1 (MQ=255)
                           gAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAAt                                                                                                     <  1:125086/141‑1 (MQ=255)
                               aCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAACGACCCTCATCAACACAAAACTACAACTTGTTAAAAAATCTCACCAACTCTTCATTAAACAGGTTAAAAAGaaa                                                                                                 >  1:38113/1‑141 (MQ=255)
                                    cacCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCttt                                                                                            <  2:397006/141‑1 (MQ=255)
                                    cacCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCttt                                                                                            >  1:270777/1‑141 (MQ=255)
                                    cacCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCttt                                                                                            <  2:149467/141‑1 (MQ=255)
                                     acCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCtttt                                                                                           <  1:114878/141‑1 (MQ=255)
                                        gTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgtt                                                                                        >  1:367665/1‑141 (MQ=255)
                                        gTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgtt                                                                                        >  1:180329/1‑141 (MQ=255)
                                           aaaatgtaacatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTttgtttgt                                                                                     <  1:188643/141‑1 (MQ=255)
                                           aaaatgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTttgtttgt                                                                                     <  1:219075/141‑1 (MQ=255)
                                           aaaatgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTttgtttgt                                                                                     <  1:218995/141‑1 (MQ=255)
                                                         gAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAAt                                                                       <  2:258673/141‑1 (MQ=255)
                                                         gAAGATAAACAGTGCTGCTAAAGTCAAATTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCACTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAAt                                                                       <  2:143690/141‑1 (MQ=255)
                                                                            aaaGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa                                                    >  2:29036/1‑141 (MQ=255)
                                                                              aGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGAcc                                                  <  2:297020/141‑1 (MQ=255)
                                                                               gTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCt                                                 >  1:132971/1‑141 (MQ=255)
                                                                                      tCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGat                                          >  2:42759/1‑141 (MQ=255)
                                                                                      tCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGat                                          >  2:285685/1‑141 (MQ=255)
                                                                                      tCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGat                                          <  2:276386/141‑1 (MQ=255)
                                                                                                             gATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACtt                   <  1:239999/141‑1 (MQ=255)
                                                                                                                  tGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGc              <  1:310090/141‑1 (MQ=255)
                                                                                                                  tGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGc              <  1:220690/141‑1 (MQ=255)
                                                                                                                  tGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGc              <  1:1948/141‑1 (MQ=255)
                                                                                                                              tGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGtt  >  1:49972/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                     
ACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTT  >  USA300TCH1516_ALE/2010230‑2010496

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