Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,756,983 A→G 9.9% intergenic (‑85/+15) USA300HOU_2612 ← / ← nrdG hypothetical protein/anaerobic ribonucleotide reductase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,756,9830AG 9.9% 132.0 / 6.8 51intergenic (‑85/+15)USA300HOU_2612/nrdGhypothetical protein/anaerobic ribonucleotide reductase small subunit
Reads supporting (aligned to +/- strand):  ref base A (24/21);  new base G (3/2);  total (27/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 = 6.98e-01

ATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGAATAAGTGTTGTATAAATAGAC  >  CP000730/2756845‑2757119
                                                                                                                                          |                                                                                                                                        
atatGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACt                                                                                                                                        <  2:418164/141‑1 (MQ=255)
   tGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCt                                                                                                                                     >  1:492388/1‑141 (MQ=255)
    gCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTa                                                                                                                                    >  1:86845/1‑141 (MQ=255)
       aCTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATa                                                                                                                                 >  2:143026/1‑141 (MQ=255)
            tAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACatat                                                                                                                            <  1:189162/141‑1 (MQ=255)
              gTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTc                                                                                                                          >  1:258725/1‑141 (MQ=255)
               tATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCa                                                                                                                         <  1:360002/141‑1 (MQ=255)
                             tGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATAtt                                                                                                           >  1:186116/1‑141 (MQ=255)
                                                     atatatTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAATGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTg                                                                                   >  2:218535/1‑141 (MQ=255)
                                                      tatatTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGa                                                                                  <  1:50038/141‑1 (MQ=255)
                                                          tttttATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGAttgt                                                                              >  1:75631/1‑141 (MQ=255)
                                                                taAAGAATTGCTTTTATTAACTTTCAATTATGTACATAACCTTAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACa                                                                        >  1:366883/1‑141 (MQ=255)
                                                                 aatgAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACAt                                                                       <  2:434273/138‑1 (MQ=255)
                                                                 aaaGAATTTCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACAt                                                                       <  2:218736/141‑1 (MQ=255)
                                                                    gAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCta                                                                    >  1:502604/1‑141 (MQ=255)
                                                                     aaTTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATAGATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCtat                                                                   <  2:344554/141‑1 (MQ=255)
                                                                      aTTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCtata                                                                  >  1:196365/1‑141 (MQ=255)
                                                                         gCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGACTGTAATACATCTATAATg                                                               >  2:198126/1‑141 (MQ=255)
                                                                          cTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGc                                                              >  2:17792/1‑141 (MQ=255)
                                                                           ttttATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCg                                                             <  2:86845/141‑1 (MQ=255)
                                                                           ttttATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCg                                                             >  1:30880/1‑141 (MQ=255)
                                                                                ttaACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGAt                                                        <  1:435012/141‑1 (MQ=255)
                                                                                   aCTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTa                                                     <  1:220369/141‑1 (MQ=255)
                                                                                    cTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGCCTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCGTGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTaa                                                    <  2:16411/141‑1 (MQ=255)
                                                                                    cTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTaa                                                    >  1:375000/1‑141 (MQ=255)
                                                                                    cTTTCAATTATGTACCTAACCTAAAAAGAAGCAAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTaa                                                    >  2:66378/1‑141 (MQ=255)
                                                                                          aTTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAAcc                                              >  1:301890/1‑141 (MQ=255)
                                                                                           ttATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGAGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCt                                             <  2:282723/141‑1 (MQ=255)
                                                                                              tGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTta                                          >  1:248914/1‑141 (MQ=255)
                                                                                              tGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACAAATTCAACTAACTATATATTCAATCTTACGAGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTta                                          >  1:249151/1‑141 (MQ=255)
                                                                                                  cctaacctaaAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATg                                      <  2:267828/141‑1 (MQ=255)
                                                                                                   ctaacctaaAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATgg                                     <  1:327689/141‑1 (MQ=255)
                                                                                                   ctaacctaaAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATgg                                     <  1:244911/141‑1 (MQ=255)
                                                                                                      acctaaAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTaa                                  <  2:204781/141‑1 (MQ=255)
                                                                                                      acctaaAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTaa                                  <  2:417690/141‑1 (MQ=255)
                                                                                                        ctaaAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAAt                                <  1:376262/141‑1 (MQ=255)
                                                                                                                 agCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGa                       >  2:369568/1‑141 (MQ=255)
                                                                                                                 agCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGa                       <  2:315155/141‑1 (MQ=255)
                                                                                                                 agCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATATATAATGCGTTGATTTAAAGAACCTTTATCTGGTAACTCAGGTTTGa                       >  2:460364/1‑141 (MQ=255)
                                                                                                                        gCAACGAATGTTACCTTGACTTCTAATACATATTCAACTCACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGAATAAgtg                <  1:428691/141‑1 (MQ=255)
                                                                                                                        gCAACGAATGTTACCTTGACTTCTAATACATATTCAACTCACAATATTTTCAATCATACGCGAATACGAGAGTGATTGTTGTACACCTATAATGCGTTGATTTAAAGTACCTTTATGTGGTAAATCAGGTTTGAATAAgtg                <  2:480438/141‑1 (MQ=255)
                                                                                                                           acattcgttgccTTGGCTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGAATAAGtgttg             >  1:168129/11‑141 (MQ=255)
                                                                                                                           aCGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGAATAAGtgttg             >  1:348691/1‑141 (MQ=255)
                                                                                                                           aCGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGAATAAGtgttg             >  1:398906/1‑141 (MQ=255)
                                                                                                                              aaTGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGAATAAGTGTTGtat          >  1:99065/1‑141 (MQ=255)
                                                                                                                                tGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGAATAAGTGTTGTATaa        <  2:308220/141‑1 (MQ=255)
                                                                                                                                 gttgccTTGGCTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGtt                          <  1:54172/118‑1 (MQ=255)
                                                                                                                                 gttgccTTGGCTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGtt                          >  2:54109/5‑122 (MQ=255)
                                                                                                                                 gttgccTTGGCTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGtt                          >  2:54172/5‑122 (MQ=255)
                                                                                                                                 gTTGCCTTGGCTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGtt                          <  1:54109/122‑1 (MQ=255)
                                                                                                                                   tACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGAATAAGTGTTGTATAAGTa     >  2:434052/1‑141 (MQ=255)
                                                                                                                                      atttaCTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAACTCAGGTTTGAATAAGTGTTGTATAAATAGAc  <  2:129509/137‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                        
ATATGCAACTTTTAGTATTGGTATCAAAATGATTGTTAGGTCATATTCTATCAATATATTTTTATAAAGAATTGCTTTTATTAACTTTCAATTATGTACCTAACCTAAAAAGAAGCCAAGGCAACGAATGTTACCTTGACTTCTAATACATATTCAACTAACTATATATTCAATCATACGCGCATGCGAGAGTGATTGTTGTACATCTATAATGCGTTGATTTAAAGAACCTTTATATGGTAAATCAGGTTTGAATAAGTGTTGTATAAATAGAC  >  CP000730/2756845‑2757119

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