Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,917 C→T 50.1% intergenic (‑1000/‑168) USA300HOU_1794 ← / → USA300HOU_1796 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,9170CT50.1% ‑4.4 / 75.7 56intergenic (‑1000/‑168)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (15/13);  new base T (21/7);  total (36/20)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

TGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTA  >  CP000730/1929778‑1930053
                                                                                                                                                |                                                                                                                                        
tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCTGTGCCa                                                                                                                                         >  1:152367/1‑141 (MQ=255)
 gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                                                                        >  2:164834/1‑141 (MQ=255)
     ttgCTTGTATATCCACAGATGGTATAATAATAGAGTCGCCTATCTCTCCGGCGTCAATT‑‑GACGAAGAGAGGAGGTTCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                   <  2:103232/140‑1 (MQ=255)
        cTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttt                                                                                                                                 >  1:14277/1‑141 (MQ=255)
         ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                               <  1:54762/141‑5 (MQ=255)
             atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc                                                                                                                            <  1:285647/141‑1 (MQ=255)
              tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                                                                                                          >  1:122592/9‑141 (MQ=255)
                 tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCAAATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                       >  1:91212/6‑141 (MQ=255)
                  tGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt                                                                                                                       >  2:108555/1‑141 (MQ=255)
                  cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGCAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                      >  1:152294/5‑141 (MQ=255)
                  cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                      >  2:42977/5‑141 (MQ=255)
                   gCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGAAACTTATTTTTGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTg                                                                                                                      >  1:77367/1‑141 (MQ=255)
                    cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                    <  2:139902/139‑1 (MQ=255)
                    cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                    >  2:329878/3‑141 (MQ=255)
                     cgttgattgtataGTTGAGTCGTCTATCTCGCAGGCGTCAATTTAGACGCAGAGAGGAGGTTTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCAAATCATAGCTCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                    <  2:136383/131‑1 (MQ=255)
                     agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                   >  1:40252/2‑141 (MQ=255)
                         ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCACTT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATGAGTTGCTGTGCTGTTGCGTATTTTTCTTGTTGGCttata                                                                                                               >  1:11708/3‑137 (MQ=255)
                         ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                               >  1:310193/3‑141 (MQ=255)
                         ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TGGACGAAGAGAGGAGGTGCATT‑AC‑T‑C‑TGCTAATCATTTTCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                               <  2:94460/139‑1 (MQ=255)
                          gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                              >  2:300018/2‑141 (MQ=255)
                          atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                               >  1:289981/1‑141 (MQ=255)
                           tataATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                              <  1:321193/141‑1 (MQ=255)
                             taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AC‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                                                                            >  1:128062/1‑141 (MQ=255)
                               ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                         <  1:283091/141‑1 (MQ=255)
                               ataGTAGGGTCGCCTATCTCTCAGGCGTTAATTTATACGCAGGGAGGAGGTATA‑T‑AA‑G‑G‑GGATGCATATTTTCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTATGtgggag                                                                                                          <  2:252941/141‑6 (MQ=255)
                                  gtagAGTCTGCTATCTCTCAGGCGTCAATTTATACGCAGATAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCaa                                                                                                       <  1:31861/141‑1 (MQ=255)
                                     gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                   >  1:194162/1‑137 (MQ=255)
                                      agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTAATTGGCTTAGTAAACGCaacgaat                                                                                                  >  1:101432/1‑136 (MQ=255)
                                       gccGCCTATCTCTCAGGCGTCCATT‑‑GACGAAGTGAGGAGGTGCA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                 <  1:111180/139‑7 (MQ=255)
                                         cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaagaat                                                                                               >  2:74143/1‑133 (MQ=255)
                                          gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                                                                               <  1:236510/141‑1 (MQ=255)
                                          gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                              >  2:290875/1‑132 (MQ=255)
                                          gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                              >  1:131133/1‑132 (MQ=255)
                                          gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                              >  1:131138/1‑132 (MQ=255)
                                           ccTATCTCTCAGGCGTCAATTTAGACGCAGCGAGGCGGCGTT‑C‑AA‑T‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAGTACAAGATa                                                                                              >  2:219080/1‑141 (MQ=255)
                                           ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                             >  2:140636/1‑131 (MQ=255)
                                            cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TCAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                            <  2:169084/141‑12 (MQ=255)
                                             tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaatgattaat                                                                                           >  2:329713/1‑129 (MQ=255)
                                             tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                           >  1:119945/1‑129 (MQ=255)
                                              aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                          >  2:52411/1‑128 (MQ=255)
                                              aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                          >  1:110041/1‑128 (MQ=255)
                                              aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                          >  2:84230/1‑128 (MQ=255)
                                                   tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACat                                                                                      >  1:284260/1‑141 (MQ=255)
                                                        cGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACAtatacc                                                                                 >  1:758/1‑139 (MQ=255)
                                                                           aggaggTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcccc                                                              >  1:254420/1‑141 (MQ=255)
                                                                            ggaggTTTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                             <  2:256696/141‑1 (MQ=255)
                                                                              aggTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa                                                           <  2:242382/141‑1 (MQ=255)
                                                                                 tgtA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACtt                                                        >  1:49376/1‑140 (MQ=255)
                                                                                                     ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAATAATCCCCTCACTACCGCAAATAGTGAgg                                         <  1:211507/141‑1 (MQ=255)
                                                                                                          ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGTCAAATGATACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                    <  1:242614/141‑1 (MQ=255)
                                                                                                          ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAGTAGTGGGGGGagt                                    >  1:32641/1‑139 (MQ=255)
                                                                                                          ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                    <  1:51801/141‑1 (MQ=255)
                                                                                                                      atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAACCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTa                        >  2:206197/1‑141 (MQ=255)
                                                                                                                              aGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTa                >  2:312972/1‑141 (MQ=255)
                                                                                                                                     aGTGGCTGTGCCATTGCGTTTTTTTCATATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAACCCCCTCACTATCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATATATATTTTCGt         <  1:244606/141‑1 (MQ=255)
                                                                                                                                            gtgCCATGGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa  <  1:88590/141‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                        
TGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTA  >  CP000730/1929778‑1930053

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