Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,727,823 C→T 8.5% intergenic (‑42/‑176) USA300HOU_2582 ← / → pyrD hypothetical protein/dihydroorotate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,727,8230CT 8.5% 143.6 / 2.7 47intergenic (‑42/‑176)USA300HOU_2582/pyrDhypothetical protein/dihydroorotate oxidase
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base T (2/2);  total (25/22)
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.06e-01

GTTTCGACTCGATATGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAATATATTTCACATATATTAAGTTGAGGAT  >  CP000730/2727684‑2727962
                                                                                                                                           |                                                                                                                                           
gTTTCGACTCGATATGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAGTCTACTTACTTTATACTATAATATTGcc                                                                                                                                            >  1:17173/1‑141 (MQ=255)
 tttCGACTCGATATGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCt                                                                                                                                           <  1:196385/141‑1 (MQ=255)
  ttCGACTCGATATGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTa                                                                                                                                          >  2:374116/1‑141 (MQ=255)
    cGACTCGATATGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAAt                                                                                                                                        <  2:100601/141‑1 (MQ=255)
      aCTCGATATGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGTCTAATTa                                                                                                                                      >  2:253291/1‑141 (MQ=255)
      aCTCGATATGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGTCTAATTa                                                                                                                                      >  2:253304/1‑141 (MQ=255)
        tCGATATGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAgg                                                                                                                                    >  1:305411/1‑141 (MQ=255)
             atGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGgatgat                                                                                                                               >  1:302282/1‑141 (MQ=255)
               gCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGAttt                                                                                                                             >  1:364248/1‑141 (MQ=255)
                cTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTTCCTAATTAGGATGATTTa                                                                                                                            >  1:246731/1‑141 (MQ=255)
                    gTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCt                                                                                                                        <  1:91729/141‑1 (MQ=255)
                    gTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCt                                                                                                                        <  2:121890/141‑1 (MQ=255)
                     tttACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCATTCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCtt                                                                                                                       <  1:183849/141‑1 (MQ=255)
                      ttACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTa                                                                                                                      >  2:70568/1‑141 (MQ=255)
                              ccaccacTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGtt                                                                                                              >  2:70941/1‑141 (MQ=255)
                              ccaccacTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGtt                                                                                                              >  2:321120/1‑141 (MQ=255)
                                         tGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTg                                                                                                   <  1:289626/141‑1 (MQ=255)
                                          ggTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGc                                                                                                  >  1:66712/1‑141 (MQ=255)
                                           gTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCg                                                                                                 <  2:266887/141‑1 (MQ=255)
                                                  ctctTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTCTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAAttt                                                                                          >  2:94842/1‑141 (MQ=255)
                                                    ctTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTc                                                                                        <  1:31228/141‑1 (MQ=255)
                                                             cTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTATAAACGTTATTCAAACTTGCGTCAATTTTCAAAACtttt                                                                               <  1:293350/141‑1 (MQ=255)
                                                             cTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACtttt                                                                               <  2:307969/141‑1 (MQ=255)
                                                              ttttAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACttttt                                                                              <  1:340209/141‑1 (MQ=255)
                                                                 tAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGt                                                                           <  1:241102/141‑1 (MQ=255)
                                                                 tAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGt                                                                           >  2:144490/1‑141 (MQ=255)
                                                                   aaGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTAAAAACTTTTTAGttt                                                                         >  2:258311/1‑141 (MQ=255)
                                                                          aTGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAg                                                                  >  2:206384/1‑141 (MQ=255)
                                                                           tGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAgg                                                                 <  1:7779/141‑1 (MQ=255)
                                                                               ctAACCATTGTTCAGTCATACAAACACCTTCCAGTCTACCGACTTTACAGTATAATAGTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCCTCAATTTGCAATACTGGTTAGTTTTAACTATGagag                                                             >  2:35751/1‑138 (MQ=255)
                                                                                      ttGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATtatgtat                                                      <  1:149409/141‑1 (MQ=255)
                                                                                      ttGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATtatgtat                                                      <  1:149416/141‑1 (MQ=255)
                                                                                       tGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATtatgtatg                                                     >  2:79889/1‑141 (MQ=255)
                                                                                                      acacCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTAt                                      >  1:291753/1‑141 (MQ=255)
                                                                                                        acCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAATCGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTa                                    <  2:246731/141‑1 (MQ=255)
                                                                                                         ccTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTaa                                   >  1:183155/1‑141 (MQ=255)
                                                                                                          cTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTaaa                                  >  2:205837/1‑141 (MQ=255)
                                                                                                              cAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGaa                              <  2:18416/141‑1 (MQ=255)
                                                                                                               aaTCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAat                             <  2:302282/141‑1 (MQ=255)
                                                                                                                aTCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAata                            <  1:209157/141‑1 (MQ=255)
                                                                                                                  cTACTTACTTTATACTATAATATTGTCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAatata                          <  1:253304/141‑1 (MQ=255)
                                                                                                                  cTACTTACTTTATACTATAATATTGTCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAatata                          <  1:253291/141‑1 (MQ=255)
                                                                                                                              taCTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAATATATTTCACatatat              >  1:363696/1‑141 (MQ=255)
                                                                                                                               aCTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAATATATTTCACATATAtt             <  2:4439/141‑1 (MQ=255)
                                                                                                                                         tGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAATATATTTCACATATATTAAGTTGAGGa   >  1:347411/1‑141 (MQ=255)
                                                                                                                                          gCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAATATATTTCACATATATTAAGTTGAGGAt  >  2:164681/1‑141 (MQ=255)
                                                                                                                                          gCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAATATATTTCACATATATTAAGTTGAGGAt  >  2:32774/1‑141 (MQ=255)
                                                                                                                                          gCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAATATATTTCACATATATTAAGTTGAGGAt  >  2:165612/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                           
GTTTCGACTCGATATGCTTCGTTTACATCACCACCACTCACTGGTGAAATCTCTTTAATATCTTTTAAAGGTAAATGCTCTAACCATTGTTCATTCATACAAACACCTTCCAATCTACTTACTTTATACTATAATATTGCCTAATTAGGATGATTTAATCTTAAAAACGTTATTCAAACTTGCGTCAATTTTCAAAACTTTTTAGTTTTAACTAGGAAATTATGTATGGAAATGCTAGTCTATTAAAAGAATATATTTCACATATATTAAGTTGAGGAT  >  CP000730/2727684‑2727962

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