Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 990,149 C→G 33.0% S14C (TCC→TGC)  oppD_2 → Oligopeptide transport ATP‑binding protein OppD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE990,1490CG33.0% 36.6 / 25.2 31S14C (TCC→TGC) oppD_2Oligopeptide transport ATP‑binding protein OppD
Reads supporting (aligned to +/- strand):  ref base C (9/11);  new base G (5/5);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TTTTACTTATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTAACAACAAAAGCAATTACAAAATTA  >  USA300TCH1516_ALE/990036‑990282
                                                                                                                 |                                                                                                                                     
ttttACTTATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTGGATATTACAGGAGGGGAAGAgc                                                                                                             >  2:22887/1‑140 (MQ=255)
   tACTTATTTAGTGAAAGAATACGAGATGCAATGGAACCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTgcagg                                                                                                          <  1:11532/140‑1 (MQ=255)
    acctATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCACATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTgcaggc                                                                                                         <  2:316772/137‑1 (MQ=255)
      ttATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTgcaggcag                                                                                                       <  1:198413/140‑1 (MQ=255)
         tttAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGAAGAGGAAGTGCAGGCAGTga                                                                                                    >  2:185805/1‑140 (MQ=255)
              tGATGGATTACGTTATTCATATGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGTAAGTGCAGGCAGTGAGAGGc                                                                                               <  2:70906/140‑1 (MQ=255)
              tGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGc                                                                                               <  2:48038/140‑1 (MQ=255)
                      tACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGAttt                                                                                       <  2:248503/140‑1 (MQ=255)
                                 tttGATCTGAAAATGCGTAAATAAAAAGGGGGCATAGCAGATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACa                                                                            <  1:34620/140‑1 (MQ=255)
                                 ttggATCCGATAATGCGTAAATAGAGAGGGGGTATTGCATATGACTGAATGAATATTAGAAGTAAAGGATTTGCATGTTTCCTTTGATATTACAGCAGGGGGAGTGCAGGCAGTGAGAGGCGTAGTTTTTTATTTGAACa                                                                            <  2:170438/137‑1 (MQ=255)
                                  ttGATTCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACaa                                                                           <  1:322610/140‑1 (MQ=255)
                                      tCCTAAAATATGTAAATAAAAAGTGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTACTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAggg                                                                       <  1:138899/140‑1 (MQ=255)
                                            aaTGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGTTATTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACa                                                                 <  1:245164/140‑1 (MQ=255)
                                             aTGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACGAAGGGGACACAt                                                                >  2:222050/1‑140 (MQ=255)
                                             aTGCGTAAATAAAAAGGGGGCATAGCATATGACAGAAAGAAGATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGACGTGCAGGCAGTGAGAGGCGTAGATTTTTAGTTGAACAAAGGGGAAACAt                                                                >  2:12861/1‑140 (MQ=255)
                                              tGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACAtt                                                               >  1:298793/1‑140 (MQ=255)
                                                cGTAAAAAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTgg                                                             <  2:215178/140‑1 (MQ=255)
                                                agaaaaaaaaaAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTgg                                                             <  2:63157/137‑1 (MQ=255)
                                                  taaataaaAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCa                                                           <  1:43180/140‑1 (MQ=255)
                                                         aaaGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAAttgttg                                                    >  1:336754/1‑140 (MQ=255)
                                                                        tatGACTGAAATAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAAtcaggttcag                                     <  2:315482/140‑1 (MQ=255)
                                                                           taCTGTAAGAATATTAGAGGTAAATGATGTGCATGTTTCCTTTGATATTACAGCAGGGGAAGGGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGCATCAGGTTCAGGTa                                  <  2:265530/139‑1 (MQ=255)
                                                                              tGAAAGAATATTAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAAt                               >  1:306211/1‑140 (MQ=255)
                                                                               gAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATc                              >  2:295592/1‑140 (MQ=255)
                                                                               gAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTAAGGTAAATc                              >  2:295525/1‑140 (MQ=255)
                                                                                     atatTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTaac                        <  1:14715/140‑1 (MQ=255)
                                                                                        ttAGAAGTAAATGATTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTaacaac                     >  1:58468/1‑140 (MQ=255)
                                                                                             aGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTAAAAACAAAAg                >  1:128748/1‑140 (MQ=255)
                                                                                                 aaTGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGGAACATTGGCAATTGATGGTGAATCAGGTTCAGGTAAATCTGTAACAACAAAAGCAAt            >  1:68506/1‑140 (MQ=255)
                                                                                                     aTTTGCATGTTTGCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTAACAACAAAAGCAATTACa        >  1:180373/1‑140 (MQ=255)
                                                                                                           aTGTTTCCTTTGATATTACAGCAGCGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGGAACAACAAAAGAAATTACAAAatta  >  1:226145/1‑140 (MQ=255)
                                                                                                                 |                                                                                                                                     
TTTTACTTATTTAGTGATGGATTACGTGATGCATTTGATCCGAAAATGCGTAAATAAAAAGGGGGCATAGCATATGACTGAAAGAATATTAGAAGTAAATGATTTGCATGTTTCCTTTGATATTACAGCAGGGGAAGTGCAGGCAGTGAGAGGCGTAGATTTTTATTTGAACAAAGGGGAAACATTGGCAATTGTTGGTGAATCAGGTTCAGGTAAATCTGTAACAACAAAAGCAATTACAAAATTA  >  USA300TCH1516_ALE/990036‑990282

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