Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,726,995 G→C 7.6% H263D (CAT→GAT)  USA300HOU_2582 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,726,9950GC 7.6% 152.5 / 2.0 49H263D (CAT→GAT) USA300HOU_2582hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/24);  new base C (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCATTGTAAAAATATCTGTGCATATTATCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTAAAACCACCAAATACCGTTGTAATACCGATATCGAATTC  >  CP000730/2726859‑2727130
                                                                                                                                        |                                                                                                                                       
tctcCTTGTAAAAATATCTGTGCATATTATCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACTTCTCACCAAATTTCAATAAATGGAcc                                                                                                                                     <  2:335539/141‑1 (MQ=255)
tctcATTGTAAAAATATTTGTGCATTTTATCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGAcc                                                                                                                                     <  2:337923/141‑1 (MQ=255)
tctcATTGTAAAAATATCTGTGCATATTATCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGAcc                                                                                                                                     <  1:6969/141‑1 (MQ=255)
 ctcATTGTAAAAATATCTGTGCATATGATCTATTTATCTCAATCTACCGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTGTCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCa                                                                                                                                    <  1:344454/141‑1 (MQ=255)
 ctcATTGTAAAAATATCTGTGCATATGATCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTCTCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCa                                                                                                                                    <  1:344766/141‑1 (MQ=255)
      tGTAAAAAGATCTGTGCATATTATCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCAGAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCaat                                                                                                                               <  1:533305/141‑1 (MQ=255)
         aaaaaTATCTGTGCATATTATCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTAAATCTCACCAAATTTCAATAAATGGACCAtaaataaa                                                                                                                            >  2:178223/1‑141 (MQ=255)
          aaaaTATCTGTGCATATTATCTATTTATCTCAAACTAGCGTCATAACCAACTTGTTGTATCATGTAAAATCTTAGCCATATAATACGCAACACTATCACGTTACATCTCACCAAATTTCAATAAATGGACCATCaataaat                                                                                                                           <  2:325545/141‑1 (MQ=255)
           aaaTATCTGTGCATATTATCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCaataaata                                                                                                                          <  2:334851/141‑1 (MQ=255)
           aaaTATCTGGGCATATTATCTATTTATCTCAATCTAACGTNTTAAACAACTTGTTGTATCATGTACAATTTTATACATAGAATGCGCAACACTACAACGGTACNTCTCTCCAAATGTCAATAACTATACGACCaataaaaa                                                                                                                          >  2:357371/1‑139 (MQ=255)
                     gCATTTTATCTATTTATCTCAATCTAACGTCTTAACTAACCTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATCGACCATCAATAAATATAAACGATaa                                                                                                                >  1:401583/1‑141 (MQ=255)
                      gatatTATCTATTTATCTCAATCTAACGTCTTAACTAATTGGTTGTATCTTGTAAAATCTTATACATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATaaa                                                                                                               <  2:47881/140‑1 (MQ=255)
                         attatCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAAtt                                                                                                            <  1:368044/141‑1 (MQ=255)
                           tatCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCa                                                                                                          <  1:138838/141‑1 (MQ=255)
                            atCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATATTCTCATGCTAAAGCCATAGATTTTTTTGCACTTGCACGGTACATATCAGCAATTTTCGATAGATGGACCAGCAATAAATATAAACGAGAAAATTCaa                                                                                                         >  1:554977/1‑141 (MQ=255)
                              ttATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGt                                                                                                       <  2:502709/140‑1 (MQ=255)
                              cTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGt                                                                                                       >  2:546928/1‑141 (MQ=255)
                                      ctcAATCTAACGTCTTAACTAAATTGTTGTGTCTTGTAACATCTTATCCATAGAATGCGCAAACCTATCACGGTACATATCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATCACATTCTATAGGa                                                                                               >  1:214882/1‑141 (MQ=255)
                                      ctcAATCTAACGTCTGAACTAAGTGGTTGTATCTTGTAAGATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACAATCAATAAATATAAACGATAAAATTAAAGTCTATAgag                                                                                               >  2:557524/1‑139 (MQ=255)
                                        ccaTCTATCGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATg                                                                                             <  1:186666/139‑1 (MQ=255)
                                        cAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAGTTTCAATAAATGGACCATCAATAAATATAAACGAAAAAATTCAAGTATATAGGaaa                                                                                             >  2:72645/1‑139 (MQ=255)
                                           tCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATTCAc                                                                                          >  1:229250/1‑141 (MQ=255)
                                              aaCGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCtt                                                                                       <  2:334280/141‑1 (MQ=255)
                                                cGTCTTAACTAACGTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCtttt                                                                                     >  1:219786/1‑141 (MQ=255)
                                                      aactaactTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGACTATAGAATGCACCTTTTGCGAGt                                                                               >  2:181913/1‑141 (MQ=255)
                                                        ctaactTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTgg                                                                             <  1:228638/141‑1 (MQ=255)
                                                            ctTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATaa                                                                         <  2:554977/141‑1 (MQ=255)
                                                            ctTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATaa                                                                         >  1:174475/1‑141 (MQ=255)
                                                                  gtATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATCATGtttt                                                                   >  2:332863/1‑140 (MQ=255)
                                                                  gtATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTtta                                                                   >  2:514850/1‑141 (MQ=255)
                                                                           aaaaTCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGca                                                          <  2:47233/141‑1 (MQ=255)
                                                                           aaaaCCGTATCCCTAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGca                                                          <  2:562245/141‑1 (MQ=255)
                                                                                 ttATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATaa                                                    >  1:465322/1‑141 (MQ=255)
                                                                                   aTCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATaaaa                                                  >  2:233937/1‑141 (MQ=255)
                                                                                   aTCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATCGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATACTGTTTATTATACGCATCATaaaa                                                  <  2:401583/141‑1 (MQ=255)
                                                                                       aTAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATATTATGCACCTTTTGCGAGTGGAAACTGGTTATTATACGCATCATAAAATTCg                                              >  1:212427/1‑141 (MQ=255)
                                                                                          gAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATACATCGACCATCAAGAAATACAAACGATAACACTACAGTCTCTAGGATGCACCGTTTGCGAGGGGCCAATGTTTATTATAAGCAACACAAAATTGGCTc                                           >  2:124764/1‑141 (MQ=255)
                                                                                           aaTGCTCATCACTAACACGTTACACCTCACCAAATTTCTCGATATGGACCATCAATAAACATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATATTGTTTATTATACGCATCATAAAATTCGCTCg                                          <  1:432793/141‑1 (MQ=255)
                                                                                           aaTGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCg                                          <  2:533912/141‑1 (MQ=255)
                                                                                            aTGCGCAACCCTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGACTGCACCTTTTGCGAGTGGATAAGGTTTATTATACGCATCATAAAATTCGccagt                                         >  2:247138/1‑137 (MQ=255)
                                                                                              gcgcAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGAACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTaa                                       <  1:84147/141‑1 (MQ=255)
                                                                                               cgcAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTaaa                                      <  1:78556/141‑1 (MQ=255)
                                                                                                      acaTCACGGTACAACTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTATATAGAATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTAAAAccacca                               <  2:262731/139‑1 (MQ=255)
                                                                                                         tCACGGTACATCTCACCAAATTTCAATACATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTAAAACCACCAAAt                            <  1:377989/141‑1 (MQ=255)
                                                                                                         tCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAACTTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTAAAACCACCAAAt                            >  1:498064/1‑141 (MQ=255)
                                                                                                                      cACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTTTTTTTGTTTATTGTTTATTATACGCATCATAAAATTCGCTCGTCAAAACACCACATACCgcggtagttc               >  2:26843/1‑132 (MQ=255)
                                                                                                                        ccAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTAAAACCACCAAATACCGTTGTAATACCg             >  1:533694/1‑141 (MQ=255)
                                                                                                                              ttCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATCGGATGCGCCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTATAACCACCACATACCGTTGTAATCCCCATATCg       >  2:530999/1‑141 (MQ=255)
                                                                                                                                   aaaaTCGACCAGCATTAAATATAAACGAGAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTAAAACCACCAAATACCGTTGTAATACCGATATCGAATtc  <  2:301375/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TCTCATTGTAAAAATATCTGTGCATATTATCTATTTATCTCAATCTAACGTCTTAACTAACTTGTTGTATCTTGTAAAATCTTATCCATAGAATGCGCAACACTATCACGGTACATCTCACCAAATTTCAATAAATGGACCATCAATAAATATAAACGATAAAATTCAAGTCTATAGGATGCACCTTTTGCGAGTGGATAATGTTTATTATACGCATCATAAAATTCGCTCGTAAAACCACCAAATACCGTTGTAATACCGATATCGAATTC  >  CP000730/2726859‑2727130

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