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,4580AT91.4% 109.8 / 4.4 35I216N (ATT→AAT) oppD_4Oligopeptide transport ATP‑binding protein OppD
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base T (16/16);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

CACATGTTTTAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAAC  >  USA300TCH1516_ALE/1416315‑1416588
                                                                                                                                               |                                                                                                                                  
cacaTGTTTTAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTGCCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTa                                                                                                                                <  2:278254/148‑1 (MQ=255)
      ttttAACTCAATCATAGTACTCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACa                                                                                                                          <  1:173862/148‑1 (MQ=255)
      ttttAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACa                                                                                                                          <  1:329115/148‑1 (MQ=255)
          aaCTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAatat                                                                                                                      <  2:86633/148‑1 (MQ=255)
          aaCTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAatat                                                                                                                      <  1:123875/148‑1 (MQ=255)
           aCTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATAtt                                                                                                                     <  2:341191/148‑1 (MQ=255)
                       tAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGAc                                                                                                         >  2:272638/1‑148 (MQ=255)
                         aTCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACtt                                                                                                       <  1:41747/148‑1 (MQ=255)
                                  tAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTAtgt                                                                                              >  1:99886/1‑148 (MQ=255)
                                                  ttAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTAAATGATTTGAACTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTAAAGTTATGTGTGATAATAAATAATg                                                                              >  2:237152/1‑148 (MQ=255)
                                                              tCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGtt                                                                  >  1:274885/1‑148 (MQ=255)
                                                                    taCGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATa                                                            <  1:102352/148‑1 (MQ=255)
                                                                    taCGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTAAATAAAAAACCTCAGTTATGTGTGACAATAAATAATGTTTGACCCTGTTCTacca                                                            <  2:103536/148‑4 (MQ=255)
                                                                       ggATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACaa                                                         >  2:218084/1‑148 (MQ=255)
                                                                         aTGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAata                                                       <  2:254307/148‑1 (MQ=255)
                                                                                      tACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTTAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAAt                                          >  1:289061/1‑148 (MQ=255)
                                                                                      tACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAAt                                          >  1:321799/1‑148 (MQ=255)
                                                                                       aCTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATc                                         <  1:67745/148‑1 (MQ=255)
                                                                                       aCTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATc                                         <  1:310592/148‑1 (MQ=255)
                                                                                       aCTTGAAATGAATGAAATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATc                                         <  1:127109/148‑1 (MQ=255)
                                                                                          tGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCAt                                      >  2:74914/1‑148 (MQ=255)
                                                                                           gAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATc                                     >  2:138193/1‑148 (MQ=255)
                                                                                           gAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATc                                     >  2:231574/1‑148 (MQ=255)
                                                                                           gAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATc                                     >  1:237232/1‑148 (MQ=255)
                                                                                             aaTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCAc                                   <  1:250051/148‑1 (MQ=255)
                                                                                             aaTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCAc                                   <  1:146006/148‑1 (MQ=255)
                                                                                                gATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTg                                >  2:59172/1‑148 (MQ=255)
                                                                                                   tAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAATATGATTGGTTTGACAATATTTCAATACATGACTTAAGTTATGGGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATACATCACTCGAtt                             >  2:226607/1‑148 (MQ=255)
                                                                                                           ccTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCa                     <  1:10952/148‑1 (MQ=255)
                                                                                                             tCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCaaa                   >  2:374287/1‑148 (MQ=255)
                                                                                                             tCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCaaa                   >  1:246696/1‑148 (MQ=255)
                                                                                                                ttCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCa                >  2:141354/1‑148 (MQ=255)
                                                                                                                          tGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAatgatg      >  1:242752/1‑148 (MQ=255)
                                                                                                                              cttcttTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAAc  >  2:48931/1‑148 (MQ=255)
                                                                                                                              cttcttTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAAc  >  2:76407/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                  
CACATGTTTTAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAAC  >  USA300TCH1516_ALE/1416315‑1416588

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