Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,845,036 T→G 8.6% K93Q (AAA→CAA)  USA300HOU_2686 ← ABC superfamily ATP‑binding cassette transporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,845,0360TG 8.6% 125.8 / 3.8 46K93Q (AAA→CAA) USA300HOU_2686ABC superfamily ATP‑binding cassette transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base T (20/22);  new base G (2/2);  total (22/24)
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.02e-01

AAATGTTCTATCATGCCAACTTTTTCAGCCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAATTAATCCGTTTATGCAATTTGCAAAAG  >  CP000730/2844897‑2845170
                                                                                                                                           |                                                                                                                                      
aaaTGTTCTATCATGCCAACTTTTTCAGCCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATTTCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCtttt                                                                                                                                       <  2:90357/141‑1 (MQ=255)
 actGTTCTATCCTCCCAAATTTTTCAGTCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTa                                                                                                                                      <  2:166/139‑1 (MQ=255)
 aaTGTTCTATCATGCCAACTTTTTCAGCCCAATAACTTACATTTTTCTTTATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTGTa                                                                                                                                      >  2:191355/1‑141 (MQ=255)
  aTGTTCTATCATGCCAACTTTTTCAGCCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTaa                                                                                                                                     <  2:45060/141‑1 (MQ=255)
   tGTTCTATCATGCCAACTTTTTCAGCCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTaaa                                                                                                                                    >  1:201441/1‑141 (MQ=255)
   tGTTCTATCATGCCAACTTTTTCAGCCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGCCCATCTGTATCTTTTaaa                                                                                                                                    >  1:221851/1‑141 (MQ=255)
        tATCATGCCAACTTTTTCAGCCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGt                                                                                                                               >  2:50870/1‑141 (MQ=255)
             tGCCAACTTTTTCAGCCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACcaa                                                                                                                          <  2:49425/141‑1 (MQ=255)
                      tttCAGCCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAg                                                                                                                 <  2:205596/141‑1 (MQ=255)
                                aTAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATg                                                                                                       >  1:58823/1‑141 (MQ=255)
                                    cTTACATTTGTCTTCATATCAGCTTGTTCAACACAAGTAGTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTaaa                                                                                                   <  1:97042/141‑1 (MQ=255)
                                     ttACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAaac                                                                                                  >  1:40666/1‑141 (MQ=255)
                                     ttACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAaac                                                                                                  >  1:222431/1‑141 (MQ=255)
                                              tcttcATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACgg                                                                                         >  1:245301/1‑141 (MQ=255)
                                                    tttcatcTTGATCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCAt                                                                                   >  2:216112/3‑141 (MQ=255)
                                                    tatCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCAt                                                                                   >  2:216125/1‑141 (MQ=255)
                                                        atcTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTg                                                                               <  2:306218/141‑1 (MQ=255)
                                                             gTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTTTCCAGCTGGATCGTTTAAGACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGattat                                                                          >  1:317086/1‑141 (MQ=255)
                                                              ttCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGattatt                                                                         >  1:183850/1‑141 (MQ=255)
                                                                 aaCACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTGATCCTCTGAGCGGTCGGACGGGAACACGTATGACTTTACGTACAACATTTGGTCTATCTTGGAAGCAACTAGCGGGTGCATCTTGAttatttat                                                                      >  1:197787/1‑141 (MQ=255)
                                                                  acacAATTATTTTCTAATAGAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTtata                                                                     >  2:229580/1‑141 (MQ=255)
                                                                  acacAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTtata                                                                     >  1:267312/1‑141 (MQ=255)
                                                                    acaATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCAATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTtatata                                                                   >  1:313130/1‑141 (MQ=255)
                                                                       attattTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAg                                                                >  2:239750/1‑141 (MQ=255)
                                                                        ttattTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGt                                                               <  2:136936/141‑1 (MQ=255)
                                                                          attTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTc                                                             >  2:195851/1‑141 (MQ=255)
                                                                          attTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTGTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTc                                                             <  1:191020/141‑1 (MQ=255)
                                                                          attTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTGTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTc                                                             <  1:191355/141‑1 (MQ=255)
                                                                          attTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCATATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTc                                                             <  2:90049/141‑1 (MQ=255)
                                                                             ttCTAATAAAATGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTTTCTTTTAAAGCTGTACCCTCAACATTATATCTATCATTTAAATATCTTACGGTTGCCTCTTGATTATTTATATATAGTTcccc                                                          <  2:24331/141‑1 (MQ=255)
                                                                              tCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCa                                                         >  1:109591/1‑141 (MQ=255)
                                                                                 aataaAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGtt                                                      >  1:241527/1‑141 (MQ=255)
                                                                                     aaaaGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAAATAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGt                                                  <  2:213538/141‑1 (MQ=255)
                                                                                                 cTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGttt                                      <  1:69047/141‑1 (MQ=255)
                                                                                                  ttCAGCTGCTGTTAAGCCTATGAACTTTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGtttt                                     <  2:274931/141‑1 (MQ=255)
                                                                                                  ttCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGtttt                                     <  2:275011/141‑1 (MQ=255)
                                                                                                        tgctgTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATTTTACCTTTAGTTTTAAATgg                               >  1:303044/1‑141 (MQ=255)
                                                                                                           tgTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAAt                            <  1:229580/141‑1 (MQ=255)
                                                                                                            gTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACTTCTAGCGGTTGCTTCTTGATTATTTATATATAGTTCCCCTGTTATGTTACCTTTTGTTTTAAATGGAAtt                           <  2:58823/141‑1 (MQ=255)
                                                                                                            gTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAAtt                           <  1:27006/141‑1 (MQ=255)
                                                                                                              tAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAATTaa                         <  1:23193/141‑1 (MQ=255)
                                                                                                              tAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAATTaa                         <  2:181467/141‑1 (MQ=255)
                                                                                                                          cTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAATTAATCCGTTTATGCa             <  1:216112/141‑1 (MQ=255)
                                                                                                                          cTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAATTAATCCGTTTATGCa             <  1:216125/141‑1 (MQ=255)
                                                                                                                                aTCTGTATCTTTTAAAACTGTACCACCAACATTAGCTCTCTCATGTAAACACATCACGATTGCATCTTTATTCTTTACATATAGTTCCCCAGTTATGTTACCTTTAGTTTTACATGCAATTACTTCGTTTCTGAAATTTGc       <  2:67137/141‑1 (MQ=255)
                                                                                                                                     tATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAATTAATCCGTTTATGCAATTTGCAAAAg  <  2:172983/141‑1 (MQ=255)
                                                                                                                                     tATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAATTAATCCGTTTATGCAATTTGCAAAAg  >  1:152672/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                      
AAATGTTCTATCATGCCAACTTTTTCAGCCCAATAACTTACATTTTTCTTCATATCATCTTGTTCAACACAATTATTTTCTAATAAAAAGGCCATATCTTCAGCTGCTGTTAAGCCTATGAACTGTCCATCTGTATCTTTTAAAACTGTACCAACAACATTAGATCTATCATGTAAACAACTAACGGTTGCATCTTGATTATTTATATATAGTTCCCCAGTTATGTTACCTTTAGTTTTAAATGGAATTAATCCGTTTATGCAATTTGCAAAAG  >  CP000730/2844897‑2845170

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: ‑