Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,691,098 A→T 10.3% S513S (TCT→TCA)  USA300HOU_2546 ← RND superfamily resistance‑nodulation‑cell division:proton (H+) antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,691,0980AT10.3% 111.3 / 4.6 39S513S (TCT→TCA) USA300HOU_2546RND superfamily resistance‑nodulation‑cell division:proton (H+) antiporter
Reads supporting (aligned to +/- strand):  ref base A (11/24);  new base T (2/2);  total (13/26)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

AATGTTGCCATTAATGACAAGACGAATCCTAAAACTGCTTTTAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTACTCTCTGCATTAGGCCCTTGCTTAGGTATCAC  >  CP000730/2690958‑2691238
                                                                                                                                            |                                                                                                                                            
aaTGTTGCCATTAATGACAAGACGAATCCTAAAACTGCTTTTAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGCCAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTCTCGCCTCATTTAACTTTTTa                                                                                                                                              >  1:41035/1‑141 (MQ=255)
       ccATTAATGACAAGACGAATCCTAAAACTGCTTTTAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACatat                                                                                                                                       <  1:306107/141‑1 (MQ=255)
            aaTGACAAGACGAATCCTAAAACTGCTTTTAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAata                                                                                                                                  <  2:290163/141‑1 (MQ=255)
            aaTGACAAGACGAATCCTAAAACTGCTTTTAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAAGGACTGATGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAata                                                                                                                                  <  2:41804/141‑1 (MQ=255)
                   aGACGAATCCTAAAACTGCTTTTAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGa                                                                                                                           <  1:159184/141‑1 (MQ=255)
                          tCCTAAAACTGCTTTTAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTg                                                                                                                    <  2:350353/141‑1 (MQ=255)
                                aaCTGCTTTTAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGa                                                                                                              <  2:98530/141‑1 (MQ=255)
                                       tttAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTTGACATATCAATATTAATGACACTTTGTCCCGAAATCTcc                                                                                                       >  2:179426/1‑141 (MQ=255)
                                       tttAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTTGACATATCAATATTAATGACACTTTGTCCCGAAATCTcc                                                                                                       >  2:449912/1‑141 (MQ=255)
                                       tttAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTcc                                                                                                       <  1:254035/141‑1 (MQ=255)
                                               tATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGaa                                                                                               >  2:283994/1‑141 (MQ=255)
                                                    gAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCat                                                                                          <  1:50300/141‑1 (MQ=255)
                                                      aTTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCatat                                                                                        <  1:135388/141‑1 (MQ=255)
                                                       ttGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGCAATCTCCGTTTTGAACCCATAtt                                                                                       >  2:318850/1‑141 (MQ=255)
                                                               aaGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTtgt                                                                               <  1:38968/141‑1 (MQ=255)
                                                                aGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTtgtg                                                                              <  1:435524/141‑1 (MQ=255)
                                                                aGACAAGCATTAATAAGAAAAATGCTAATACAACAATGACTGAGTCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTtgtg                                                                              <  1:309513/141‑1 (MQ=255)
                                                                 gACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGc                                                                             <  1:272976/141‑1 (MQ=255)
                                                                             taaGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGCGCATCTTTATGATa                                                                 >  2:374675/1‑141 (MQ=255)
                                                                              aaGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATaa                                                                >  1:149875/1‑141 (MQ=255)
                                                                                 aaaaaTGCTAATACAACAATTACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAAtct                                                             <  2:302058/141‑1 (MQ=255)
                                                                                   acctGCTAATACAACAATGACTGAGGCAAATAAAGGTATAGCATCATCTCACTTCATAGACATCACCAACTTCATCTCAAATTGACCCGACCTCTCCGTTTTGAACCCATATTTCTCTTGTGCCTCTTTATGATAATCTCg                                                           <  1:224307/138‑1 (MQ=255)
                                                                                   aaaTGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCg                                                           <  1:278499/141‑1 (MQ=255)
                                                                                             acaacaATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGc                                                 <  1:383198/141‑1 (MQ=255)
                                                                                               aacaaTGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCAc                                               <  2:278786/141‑1 (MQ=255)
                                                                                                acaaTGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTTGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACt                                              <  1:449912/141‑1 (MQ=255)
                                                                                                acaaTGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTTGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACt                                              <  1:179426/141‑1 (MQ=255)
                                                                                                acaaTGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACt                                              <  1:310296/141‑1 (MQ=255)
                                                                                                      aCTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCACGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATAATGAACTAAATCa                                        >  2:268206/1‑141 (MQ=255)
                                                                                                          aGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTg                                    <  1:349671/141‑1 (MQ=255)
                                                                                                            gCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGTACAAATCTTTCTATTGTCTATATTTATGTTATCCACGTAAATAATGCACTCAATCATTTGTa                                  >  1:40859/1‑141 (MQ=255)
                                                                                                              aaaTAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACCTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTAct                                <  2:95527/141‑1 (MQ=255)
                                                                                                                aTAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTActct                              >  2:268617/1‑141 (MQ=255)
                                                                                                                 acctttatttgCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTActctc                             >  1:157093/11‑141 (MQ=255)
                                                                                                                   aaGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTACTCTCTg                           >  2:378884/1‑141 (MQ=255)
                                                                                                                      gTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTACTCTCTGCAt                        <  2:398883/141‑1 (MQ=255)
                                                                                                                       atttgCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTACTCTCTGCAtt                       >  1:347113/5‑141 (MQ=255)
                                                                                                                            ccTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTACTCTCTGCATTAGGcc                  <  1:450881/141‑1 (MQ=255)
                                                                                                                                     aCTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTACTCTCTGCATTAGGCCCTTGCTTAg         <  1:220697/141‑1 (MQ=255)
                                                                                                                                          ttAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTACTCTCTGCATTAGGCCCTTGCTTAGGTATc    >  1:89698/1‑141 (MQ=255)
                                                                                                                                            aGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTACTCTCTGCATTAGGCCCTTGCTTAGGTATCAc  >  1:381504/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
AATGTTGCCATTAATGACAAGACGAATCCTAAAACTGCTTTTAATGGTATTAGAATTGAACGGAAGACAATCATTAATAAGAAAAATGCTAATACAACAATGACTGAGGCAAATAAAGGTATCGCCTCATTTAACTTTTTAGACATATCAATATTAATGACACTTTGTCCCGAAATCTCCGTTTTGAACCCATATTTATCTTGTGCATCTTTATGATAATCTCGTAAATCATGCACTAAATCATTTGTACTCTCTGCATTAGGCCCTTGCTTAGGTATCAC  >  CP000730/2690958‑2691238

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