Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,416,458 A→T I216N (ATT→AAT)  oppD_4 ← Oligopeptide transport ATP‑binding protein OppD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,416,4580AT92.9% 36.1 / NA 14I216N (ATT→AAT) oppD_4Oligopeptide transport ATP‑binding protein OppD
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (9/4);  total (10/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAA  >  USA300TCH1516_ALE/1416337‑1416579
                                                                                                                         |                                                                                                                         
gTAATCCCTCTTTACTTGTGTTCTATATTTAATTATACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTg                                                                                                                    >  1:177811/1‑129 (MQ=255)
gTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATAAGAATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTg                                                                                                                    >  1:145610/1‑129 (MQ=255)
             aaTTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAAt                                                                                                       <  2:197343/129‑1 (MQ=255)
                    ttCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTCACAATATTTCAATACATGAc                                                                                                >  2:30744/1‑129 (MQ=255)
                                                cGGATGCAAATGCTCATACTTGAAATGATGAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTTTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGt                                                                    >  2:174441/1‑129 (MQ=255)
                                                           gCTCATACTTGAACTGATTAATATTACCTCGTTCAATGAATTAAACTCCTTTTAAAAAATAATTGTACTGACAATATTTAAATACATGAATTAAGATATGTGAGATAATAAATAAAGTATGACCATGtt                                                         >  2:119728/1‑129 (MQ=255)
                                                                     gAAATGATTCCTATTACCTCGTTCACTGATTTGACCTTCTTTAAAAACAAAAAAAAACTAACAATATTACAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAat                                               >  1:4848/1‑129 (MQ=255)
                                                                     gAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAat                                               >  1:47029/1‑129 (MQ=255)
                                                                     gAAATGATGAATATTACATCGATCCACGACTTGACGTTCTTTTAAAACATAATTGTACTGACAACATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATATATGACCATGTTCTAATACAat                                               >  2:39660/1‑129 (MQ=255)
                                                                         tGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGc                                           <  1:95864/129‑1 (MQ=255)
                                                                                   tACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATATTTGACCATGTTCTAATACAATATGATGTAATAAAt                                 >  1:198909/1‑129 (MQ=255)
                                                                                                gATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGATGTAATAAATCCATCACTTGAtt                    <  1:216634/129‑1 (MQ=255)
                                                                                                              ttAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGCCTGATAATGAACAATGTTCACCCATGTTCTAATACAATATAAAGTCATAAATCCATCACTTGATATCCGCTCAAACCaa      >  2:217812/1‑128 (MQ=255)
                                                                                                                  aaCATAATTGTACTGACAATATGTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCaa  <  1:30744/129‑1 (MQ=255)
                                                                                                                         |                                                                                                                         
GTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAA  >  USA300TCH1516_ALE/1416337‑1416579

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