Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 162,467 A→G 12.1% L63S (TTA→TCA)  USA300HOU_0153 ← ABC superfamily ATP binding cassette transporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730162,4670AG12.1% 82.1 / 4.3 33L63S (TTA→TCA) USA300HOU_0153ABC superfamily ATP binding cassette transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base A (13/16);  new base G (2/2);  total (15/18)
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.28e-01

CCACAAAGATTGCTGCTAATAACAAATCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATAGATCACCTATTTGAGGTATTCCTATCATTAAA  >  CP000730/162327‑162595
                                                                                                                                            |                                                                                                                                
ccACAAAGATTGCTGCTAATAACAAATCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGGTGCTACCTa                                                                                                                                  >  1:9223/1‑141 (MQ=255)
ccACAAAGATTGCTGCTAATAACAAATCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTa                                                                                                                                  >  2:260994/1‑141 (MQ=255)
 cacaAAGATTGCTGCTAATAACAAATCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGAGGCTACCTaa                                                                                                                                 >  2:231141/1‑141 (MQ=255)
               ctAATAACAAATCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAGCGGATTCCACCTAATACTGTACTTACGa                                                                                                                   <  1:267071/141‑1 (MQ=255)
                   aaaCAAATCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGc                                                                                                               <  1:32171/140‑1 (MQ=255)
                      cAAATCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATCAGAACGCAATTGAAATTGAATCGATGCTACCTAATACTGTACTTACGATAGCCAt                                                                                                            <  1:63574/141‑1 (MQ=255)
                      cAAATCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCAt                                                                                                            <  1:106191/141‑1 (MQ=255)
                        aaTCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTc                                                                                                          <  1:256253/141‑1 (MQ=255)
                         aTCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCg                                                                                                         >  2:231732/1‑141 (MQ=255)
                             ggAATAGTACGAATTATATTTAAAATAAAGCGCGAGGCTGTTGGAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATg                                                                                                     <  1:116208/141‑1 (MQ=255)
                             ggAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATg                                                                                                     >  1:243428/1‑141 (MQ=255)
                              gAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATgg                                                                                                    >  1:111601/1‑141 (MQ=255)
                                aTAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGAAAGTGGTACGGTAGCGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTa                                                                                                  >  1:7732/1‑141 (MQ=255)
                                aTAGTACGAACTATATTTAAAATAAAGCGCGAGGGGATTGAAATCCACTTTTGATAACCAATCTTCCTAGCACAAAATCACGCCACTGGTATTGAAACGACGCTACCAAATACTGTCCTTACGCTGACCAGCGGAATGGTa                                                                                                  >  2:64941/1‑141 (MQ=255)
                                                 tAAAATAAAGCGCGAGTGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGt                                                                                 <  1:34384/141‑1 (MQ=255)
                                                   aaaTAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTaa                                                                               <  2:257305/141‑1 (MQ=255)
                                                            gcgAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTa                                                                      <  2:249651/141‑1 (MQ=255)
                                                                   tATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTc                                                               <  2:113584/141‑1 (MQ=255)
                                                                            ccACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAgg                                                      <  1:231732/141‑1 (MQ=255)
                                                                             cACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATCGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTTTTGTAAATACTCGAAATCAggg                                                     >  2:28687/1‑140 (MQ=255)
                                                                                  ttGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATttgttg                                                                        <  1:225005/117‑1 (MQ=255)
                                                                                  ttGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATttgttg                                                                        >  2:225005/1‑117 (MQ=255)
                                                                                   tgatgaACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGGGGAAta                                               >  2:113027/1‑140 (MQ=255)
                                                                                         aCGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTGATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCAtttgt                                         >  2:265575/1‑141 (MQ=255)
                                                                                         aCGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTGATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCAtttgt                                         >  1:52371/1‑141 (MQ=255)
                                                                                             tatTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGTTTCCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGa                                     <  1:6538/141‑1 (MQ=255)
                                                                                              atTGCTAGCAAATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGaa                                    <  2:234507/141‑1 (MQ=255)
                                                                                                        caTAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACCATGGTGTTGTAATTGTTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATAGATCACCt                          <  2:111601/141‑1 (MQ=255)
                                                                                                                     ttGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATGGATGACCTATTTGAGGTATTc             >  1:4941/1‑141 (MQ=255)
                                                                                                                      tGGTATTGAAACGATGCTACCTGATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATAGATCACCTATTTGAGGTATTcc            <  1:265575/141‑1 (MQ=255)
                                                                                                                      tGGTATTGAAACGATGCTACCTGATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATAGATCACCTATTTGAGGTATTcc            <  2:52371/141‑1 (MQ=255)
                                                                                                                      tGGTATTGAAACGATGCTACCTAATACTGTACTGACGTAATCCTTTCGAAGGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATAGATCACCTATTTGAGGTATTcc            <  1:183318/141‑1 (MQ=255)
                                                                                                                        gTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATGGATCACCTATTTGAGGTATTCCTa          >  2:45819/1‑141 (MQ=255)
                                                                                                                                aCGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATAGATCACCTATTTGAGGTATTCCTATCATTaaa  >  2:119353/1‑141 (MQ=255)
                                                                                                                                aCGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATAGATCACCTATTCGAGGTATTCCGATCATTaaa  >  2:154138/1‑141 (MQ=255)
                                                                                                                                aCGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATAGATCACCAATTTGAGGTATTCCTATCATTaaa  >  2:119361/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                
CCACAAAGATTGCTGCTAATAACAAATCTGGAATAGTACGAACTATATTTAAAATAAAGCGCGAGGGTATTGAAATCCACTTTTGATGAACGATATTGCTAGCACATAATAACGCAATTGGTATTGAAACGATGCTACCTAATACTGTACTTACGATAGCCATTCGAATGGTATCTAACATTGGCGTTGTAATTTGTTGTAAATACTCGAAATCAGGTGGAATCATTTGTTTGAATAGATCACCTATTTGAGGTATTCCTATCATTAAA  >  CP000730/162327‑162595

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