Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,799 C→A 44.7% intergenic (‑882/‑286) 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,7990CA44.7% 17.7 / 44.1 41intergenic (‑882/‑286)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (10/12);  new base A (8/10);  total (18/23)
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.60e-01

CATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGC  >  CP000730/1929662‑1929939
                                                                                                                                         |                                                                                                                                                
cATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgat                                                                                                                                               <  1:167988/141‑1 (MQ=255)
   aataaAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCgatgat                                                                                                                                            <  1:134826/141‑1 (MQ=255)
                     tttAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCt                                                                                                                          <  2:79189/140‑1 (MQ=255)
                           ttAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtctct                                                                                                                    >  2:331103/1‑141 (MQ=255)
                            tAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATGGTAGAGTCGCCTGactctc                                                                                                                   >  1:181922/1‑141 (MQ=255)
                              aCAATCTCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAg                                                                                                                 <  2:304180/141‑1 (MQ=255)
                               cAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAgg                                                                                                                >  1:374347/1‑141 (MQ=255)
                                aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                               <  1:323341/141‑1 (MQ=255)
                                            cctataaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTctaaattgac                                                                                                   <  2:333373/136‑10 (MQ=255)
                                                gactttGTACAATTAAGTGTAGAGGACTCAATCTATAAGCTAGAATCATATAAGTTATCAGAAGCAAGTGACTGTGCTTGTATATTTTACGATGAGCAAAAATNAGAGTCGACGATCTCTTAGTCTTCAAGTTAAACGCag                                                                                               >  2:626/4‑141 (MQ=255)
                                                gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGAGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                               >  1:372521/1‑141 (MQ=255)
                                                   tttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagaga                                                                                            <  1:82290/141‑1 (MQ=255)
                                                                             aaTCTATCAGCTAGAATCATATAACTTATCCGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTAttt                                                              <  1:45392/141‑1 (MQ=255)
                                                                                             tCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGc                                              <  1:148237/141‑1 (MQ=255)
                                                                                             tCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGc                                              <  2:119294/141‑1 (MQ=255)
                                                                                                aataACCTATCAGAAGCAACTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCAcc                                           <  2:371816/140‑1 (MQ=255)
                                                                                                         tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                                  <  1:331103/141‑1 (MQ=255)
                                                                                                         tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                                  >  1:352187/1‑141 (MQ=255)
                                                                                                            gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCCGTgg                               >  2:228990/1‑141 (MQ=255)
                                                                                                            atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                              <  1:258370/127‑1 (MQ=255)
                                                                                                            atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGGGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                              >  1:309528/15‑141 (MQ=255)
                                                                                                             tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                             <  2:159485/128‑1 (MQ=255)
                                                                                                             tcggaagaagatatgCTTGTATATTCACAGATGGTATAATAATATAGTCGCCTATCTCTCAGGCGTGAATT‑‑GATGACGAGAGGAGGGGCAGTAC‑G‑GT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTTGcg                             >  2:262988/14‑140 (MQ=255)
                                                                                                             aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGc                              >  1:286233/1‑141 (MQ=255)
                                                                                                              agccagcgtcggtgCTTGTCTATTTGCTGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGGCTATTTAGCCGCAGAGAGGAGGTGTA‑TAA‑T‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                             <  1:240800/130‑1 (MQ=255)
                                                                                                                gaagtagttatgCTTGGATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGCCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                          <  1:4203/131‑1 (MQ=255)
                                                                                                                gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                          <  1:4168/131‑1 (MQ=255)
                                                                                                                gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                          >  2:75261/11‑141 (MQ=255)
                                                                                                                gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCAATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                          <  1:161606/131‑1 (MQ=255)
                                                                                                                  agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTTTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                        <  2:320757/133‑2 (MQ=255)
                                                                                                                  agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                        <  1:36955/133‑2 (MQ=255)
                                                                                                                   gtagatatggttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                       <  2:157222/131‑3 (MQ=255)
                                                                                                                   gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGGA‑TTA‑C‑TTGGCTAATCATTTTCGTTCACATTATAGTACCAGTCATCCGTGGCTGTGctg                       >  2:325281/8‑139 (MQ=255)
                                                                                                                   gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                       >  1:356263/8‑139 (MQ=255)
                                                                                                                   gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTATCAGGCGTTAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                       >  1:254871/8‑139 (MQ=255)
                                                                                                                    tagagatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                      <  1:327937/135‑4 (MQ=255)
                                                                                                                    tGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGccc                       >  2:113043/1‑140 (MQ=255)
                                                                                                                      gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGCCCATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGctgttg                    >  2:133121/5‑136 (MQ=255)
                                                                                                                           gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTAATAAGTGGCTGTGCCATTGCGtt                >  2:170420/1‑141 (MQ=255)
                                                                                                                               gTATATTCACAGATGGTATAATAATAGCGTCGACTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGCGGTGCATTAC‑T‑TT‑GCTCATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata           <  2:19705/141‑7 (MQ=255)
                                                                                                                               gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatc           >  2:148393/1‑138 (MQ=255)
                                                                                                                                         cGATAATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGTTGTA‑TAA‑G‑GT‑GATGCTTATTTTCTTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc  <  2:307426/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                
CATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGC  >  CP000730/1929662‑1929939

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