Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,554,822 C→T 7.7% A70T (GCT→ACT)  USA300HOU_2413 ← ABC superfamily ATP binding cassette transporter, ABC/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,554,8220CT 7.7% 146.1 / 2.5 52A70T (GCT→ACT) USA300HOU_2413ABC superfamily ATP binding cassette transporter, ABC/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (21/27);  new base T (2/2);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

CTTAGATTCCACTACGCCTAAAGGGAGACGCATTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGAAAAGACACTGATGATAGACAACAAAATTGATAA  >  CP000730/2554685‑2554960
                                                                                                                                         |                                                                                                                                          
cTTAGATTCCACTACGCCTAAAGGGAGACGCATTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGAt                                                                                                                                         >  2:388160/1‑141 (MQ=255)
cTTAGATTCCACTACGCATAAAGGGAGACACATTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTAATATGACGTAATATGTCATCGCAGCGTTGACGCTCCAAAACAAACCACTTGTGAAATATAAGCGAt                                                                                                                                         >  1:409614/1‑141 (MQ=255)
      ttCCACTACGCCTAAAGGGAGACGCATTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGc                                                                                                                                   <  2:403722/141‑1 (MQ=255)
         cACTACGCCTAAAGTGAGACGCATTATTTTTTGAGCTAATGTTTTACGAATTTCAGATAAAATTTCATATGCCGTAAAATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAAt                                                                                                                                <  2:41258/141‑1 (MQ=255)
           cTACGCCTAAAGGGAGGCGCATTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGTGATTAAAGCAATaa                                                                                                                              >  2:167148/1‑141 (MQ=255)
                       ggAGACGCATTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCa                                                                                                                  <  1:203704/141‑1 (MQ=255)
                          gaCGCATTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATaa                                                                                                               >  1:86962/1‑141 (MQ=255)
                             gCATTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCg                                                                                                            <  2:453901/141‑1 (MQ=255)
                               aTTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGaa                                                                                                          <  1:433931/141‑1 (MQ=255)
                                   aTTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTaa                                                                                                      <  1:196076/141‑1 (MQ=255)
                                    tttttGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAAt                                                                                                     >  1:154188/1‑141 (MQ=255)
                                      tttGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCg                                                                                                   <  1:319699/141‑1 (MQ=255)
                                      tttGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCg                                                                                                   <  1:244046/141‑1 (MQ=255)
                                      tttGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCg                                                                                                   <  2:431805/141‑1 (MQ=255)
                                      tttGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCg                                                                                                   >  1:426631/1‑141 (MQ=255)
                                      tttGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCg                                                                                                   >  2:117795/1‑141 (MQ=255)
                                        tGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGta                                                                                                 >  1:62094/1‑141 (MQ=255)
                                        tGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGta                                                                                                 <  2:426631/141‑1 (MQ=255)
                                              aaTTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTa                                                                                           >  2:82575/1‑141 (MQ=255)
                                               aTTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTAt                                                                                          >  2:218197/1‑141 (MQ=255)
                                                ttgtttACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAACCCATAATCGAATTAATCGTATATGTAtt                                                                                         <  2:197559/138‑1 (MQ=255)
                                                  ttttACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGt                                                                                       <  1:388160/141‑1 (MQ=255)
                                                      aCGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAat                                                                                   >  2:249215/1‑141 (MQ=255)
                                                       cGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTaacc                                                                                  <  1:304097/141‑3 (MQ=255)
                                                          aTTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAatcatc                                                                               >  1:407745/1‑141 (MQ=255)
                                                             tCAGATAAACTTTCATATGCCGTAATATGTGATAACATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGTGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGCATTGTTAATCATCAtt                                                                            <  2:94297/141‑1 (MQ=255)
                                                                 ttAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTaaaa                                                                        <  1:321888/140‑1 (MQ=255)
                                                                  tAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGTGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGCATTGTTAATCATCATTAAaat                                                                       <  1:167148/141‑1 (MQ=255)
                                                                       tttCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAAttt                                                                  <  1:400686/141‑1 (MQ=255)
                                                                       tttCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAAttt                                                                  <  2:169211/141‑1 (MQ=255)
                                                                         tCATATGCCGTAATATGTGATAACATCGATGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTa                                                                <  2:119348/141‑1 (MQ=255)
                                                                            tatGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAAt                                                             >  2:439501/1‑141 (MQ=255)
                                                                             atGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATa                                                            >  2:34352/1‑141 (MQ=255)
                                                                               gCCGTAATATGTGATAGCATCGTTGACACTCCAAAACAACACACATGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACt                                                          <  2:213062/141‑1 (MQ=255)
                                                                                ccGTATTATGTGTTAGCATTGTTGACGTTCGAAAACAACACACTTGTGAAATAAAAGCGATCAAAATAATAAAGATAAAAAACATATCCGAAATAATCGTATAGCTATTGTACATCAGCATTAATATAATGTTAAATATTg                                                         >  2:188952/1‑141 (MQ=255)
                                                                                   tAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCTTAATCGAATTAATCGTATATGTATTGGTAATCATCATTAAAATAATTTTAAATACTGctc                                                      <  2:377009/141‑3 (MQ=255)
                                                                                   tAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGccc                                                      <  2:294915/141‑1 (MQ=255)
                                                                                   tAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCAAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGccc                                                      <  2:82687/141‑1 (MQ=255)
                                                                                         gtgATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAg                                                <  1:143477/141‑1 (MQ=255)
                                                                                          tgATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAgg                                               >  2:237790/1‑141 (MQ=255)
                                                                                                 aTCGTTGACGCTCCAAAACAACACACTTGTGAAAAATAAGCGATTAAAGCAATAAAGATATAAACCCTAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAAt                                        <  1:154615/141‑1 (MQ=255)
                                                                                                        aCGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATGAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGaaa                                 <  1:203563/141‑1 (MQ=255)
                                                                                                        aCGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATGAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGaaa                                 <  1:203559/141‑1 (MQ=255)
                                                                                                         cGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGaaaa                                >  1:396906/1‑141 (MQ=255)
                                                                                                          gCTCCAAAACAACACACTTGTGAAATATAAGTGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGCATTGTTAATCATCATTAAAATAATTTTAAAGACTGCCCAATATGGAACTAATCCAGAAAAg                               >  1:127396/1‑141 (MQ=255)
                                                                                                             ccAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGAAAAGaca                            >  1:50882/1‑141 (MQ=255)
                                                                                                                     acacacTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGAAAAGACACTgatgat                    >  2:169200/1‑141 (MQ=255)
                                                                                                                      cacacTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGAAAAGACACTGATGATa                   >  2:67601/1‑141 (MQ=255)
                                                                                                                        cacTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGAAAAGACACTGATGATAGa                 >  1:54073/1‑141 (MQ=255)
                                                                                                                         acTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGAAAAGACACTGATGATAGac                <  1:231431/141‑1 (MQ=255)
                                                                                                                             gtgAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTACATACTGCCCAATAAGGAACTAATCCAGAAAAGACACTGATGATAGacaaca            <  2:155079/141‑1 (MQ=255)
                                                                                                                                       aGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGAAAAGACACTGATGATAGACAACAAAATTGATaa  >  1:362373/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CTTAGATTCCACTACGCCTAAAGGGAGACGCATTAATTTTTGAGCTAATTTTTTACGAATTTCAGATAAAATTTCATATGCCGTAATATGTGATAGCATCGTTGACGCTCCAAAACAACACACTTGTGAAATATAAGCGATTAAAGCAATAAAGATATAAACCATAATCGAATTAATCGTATATGTATTGTTAATCATCATTAAAATAATTTTAAATACTGCCCAATAAGGAACTAATCCAGAAAAGACACTGATGATAGACAACAAAATTGATAA  >  CP000730/2554685‑2554960

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