Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,010,363 C→T 100% 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,3630CT100.0% 121.5 / NA 39C146Y (TGT→TAT) USA300HOU_RS10125Putative multidrug export ATP‑binding/permease protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/21);  total (18/21)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTG  >  USA300TCH1516_ALE/2010227‑2010497
                                                                                                                                        |                                                                                                                                      
ttCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATc                                                                                                                                    >  1:124247/1‑141 (MQ=255)
   aCGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTaa                                                                                                                                 >  1:129550/1‑141 (MQ=255)
            ttttCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATAtt                                                                                                                        >  1:211552/1‑141 (MQ=255)
              ttCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCa                                                                                                                      <  2:171141/141‑1 (MQ=255)
                   aaTCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAc                                                                                                                 >  1:103679/1‑141 (MQ=255)
                      ctaccAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCg                                                                                                              <  2:207563/138‑1 (MQ=255)
                          cAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCGATTTCACATCTAAAAAGAACATTAGGGATAGTGCAATAATAATTGCTATCTCTTGTAACCAAATATTCATTAACCCGGTTa                                                                                                          >  1:371733/1‑141 (MQ=255)
                           aaaGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaa                                                                                                         <  2:104093/141‑1 (MQ=255)
                            aaGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaaa                                                                                                        >  1:247136/1‑141 (MQ=255)
                                aaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGTCTCTAATGTACAGTGCCCTTTGCGCGGTGATATAGTCTAACCAAATATTCATTAACCCGGTTAAAATGa                                                                                                    <  1:291885/141‑1 (MQ=255)
                                 aatgTACACCGTTAAAATGTAAAATGGGAAGATAAACTGGTCTTCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaa                                                                                                   <  2:310663/138‑1 (MQ=255)
                                 aaCGTGGACGGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTAACATCTAAAAATAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCTAATATTCATTAACCCGGTTAAAATGaa                                                                                                   <  2:109955/141‑1 (MQ=255)
                                  aCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaaa                                                                                                  <  1:385470/141‑1 (MQ=255)
                                               aaatgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTtgtttgtt                                                                                     >  2:112260/1‑141 (MQ=255)
                                                aatgtaaagtgGGGAGATAAACAGTGCTGCTAAAGTCAATTTCACAGCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTc                                                                                    <  1:396602/141‑1 (MQ=255)
                                                aatgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTc                                                                                    <  2:385619/141‑1 (MQ=255)
                                                    taaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAca                                                                                <  2:110295/141‑1 (MQ=255)
                                                     aaaatgGGAAGATAAAAAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACAGAATGAGAATTTGCCGGGTGAAAATGAAATCTTTTGTTTGTTCCAcat                                                                               >  1:204623/1‑141 (MQ=255)
                                                      aaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAcatc                                                                              <  2:72578/141‑1 (MQ=255)
                                                         tgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTCTTGTTTGTTGAACATCAtt                                                                           >  1:26375/1‑141 (MQ=255)
                                                             aaGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCATATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATc                                                                       <  2:161599/141‑1 (MQ=255)
                                                                  aaaCAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActct                                                                  >  2:245606/1‑141 (MQ=255)
                                                                   aaCAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTa                                                                 <  1:157342/141‑1 (MQ=255)
                                                                   aaCAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGCTTGTTCAACATCACTAATCACTCTa                                                                 >  2:279785/1‑141 (MQ=255)
                                                                         gctgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTTTTCAACATAATTAATCACTCTAGATAtt                                                           >  1:226160/1‑141 (MQ=255)
                                                                          ctgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTa                                                          <  2:247136/141‑1 (MQ=255)
                                                                                                    aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATaa                                >  1:382334/1‑141 (MQ=255)
                                                                                                      gAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATaaaa                              >  2:195401/1‑141 (MQ=255)
                                                                                                              tGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTCCTTGACCTACGTGATTATTAGCATAAAATCTCGCAc                      <  1:112260/141‑1 (MQ=255)
                                                                                                                    gTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCAGCATCTTCTTAACAATACCTACTTCCTGATTAGCATAAAATCTCGCACTTAAAg                >  1:288334/1‑141 (MQ=255)
                                                                                                                        aataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAACATAAAATATAGCACTTAAAGCttt            >  2:345135/1‑140 (MQ=255)
                                                                                                                         ataataatTGTTATATAATCTAACCAAATATTCATTAGCCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGt           <  1:161435/141‑1 (MQ=255)
                                                                                                                         ataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGt           >  1:389804/1‑141 (MQ=255)
                                                                                                                          taataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTa          <  1:330856/141‑1 (MQ=255)
                                                                                                                            ataatTTTTATATAATCTAACCAAATATTCATTAACTTGCATAAGATTAAATCTTTTGTTTGTTCAACATCATTAATCATTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTaaa        <  1:270341/141‑1 (MQ=255)
                                                                                                                                ttGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAACTCTCTCACTTAAAGCTTGTAAAtact    >  2:158004/1‑138 (MQ=255)
                                                                                                                                 tGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGtt   <  1:165038/141‑1 (MQ=255)
                                                                                                                                 tGTTATATAATCTAACCAAATATTCATTAACCCCGATACACTCATATTTTGTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGtt   <  1:373802/141‑1 (MQ=255)
                                                                                                                                  gTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTg  <  1:61745/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTG  >  USA300TCH1516_ALE/2010227‑2010497

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