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,3630CT100.0% 39.6 / NA 12C146Y (TGT→TAT) USA300HOU_RS10125Putative multidrug export ATP‑binding/permease protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/8);  total (4/8)

CGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAG  >  USA300TCH1516_ALE/2010231‑2010483
                                                                                                                                    |                                                                                                                        
cGTGTCAATTCTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAc                                                                                                                  <  1:191321/141‑1 (MQ=255)
                     ccAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCCCATCTAAACAGCACATTATGGCTAGTGCACTAATAATTGTTATATACTCTACACCAATATTCATTAACCCGGtt                                                                                             <  1:156644/141‑1 (MQ=255)
                       aaaGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaa                                                                                           <  1:180299/141‑1 (MQ=255)
                           aaaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATg                                                                                       <  1:211402/141‑1 (MQ=255)
                                     ccGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttg                                                                             <  1:198032/141‑1 (MQ=255)
                                        ttAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgttt                                                                          >  1:294430/1‑141 (MQ=255)
                                              tgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCaa                                                                    <  1:206179/141‑1 (MQ=255)
                                                     tgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCAtt                                                             <  1:216626/141‑1 (MQ=255)
                                                                                      cacaTCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGatt                            >  1:210088/1‑141 (MQ=255)
                                                                                                     cATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAtct             >  1:335317/1‑141 (MQ=255)
                                                                                                               aGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTaaa   >  1:180090/1‑141 (MQ=255)
                                                                                                                gTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTCACCCGGTTACACTGAAATCTTTTGTTTGTTCACCATCAGTAATCCCTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAg  <  1:16183/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                        
CGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAG  >  USA300TCH1516_ALE/2010231‑2010483

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