| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | CP000730 | 35,876 | C→T | 46.3% | intergenic (+107/+1) | USA300HOU_0025 → / ← USA300HOU_0026 | hypothetical protein/IS431mec transposase | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 35,876 | 0 | C | T | 46.3% | 32.0 / 49.9 | 45 | intergenic (+107/+1) | USA300HOU_0025/USA300HOU_0026 | hypothetical protein/IS431mec transposase | 
| Reads supporting (aligned to +/- strand): ref base C (11/13); new base T (9/12); total (20/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01 | |||||||||||
AGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAC  >  CP000730/35736‑36015                                                                                                                                            |                                                                                                                                           aGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGc                                                                                                                                             >  1:216508/1‑141 (MQ=255)                 tttttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTCATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatga                                                                                                                            >  2:276640/1‑141 (MQ=255)                    ttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGc                                                                                                                         <  1:23039/141‑1 (MQ=255)                     tGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCt                                                                                                                        <  1:252774/141‑1 (MQ=255)                      tgataaaaaTGATTACTTGAACAATAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTa                                                                                                                       <  2:65442/137‑1 (MQ=255)                           aaaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAACAATAGAAGGTTCTGTTGCAAAGTAAAATAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTc                                                                                                                  >  1:329679/1‑141 (MQ=255)                             aaTGATTACTTGAACAAGGGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTGAAAAAATATAGCTACCCACTAATTTATCATGTCAGTGTTAGCTTAGCTTGCTAGCATGAGGCTAATTTCGt                                                                                                                >  2:38533/1‑141 (MQ=21)                                      ttGAAAAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGTAATGGct                                                                                                       <  2:289270/141‑2 (MQ=255)                                           cAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAt                                                                                                  <  1:58720/141‑1 (MQ=255)                                           cAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAt                                                                                                  <  1:323583/141‑1 (MQ=255)                                               agCTAAAATAAAATTGTGATCTAATAGAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGACACTCCGt                                                                                              >  2:306491/1‑141 (MQ=255)                                                    aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATc                                                                                         >  2:1121/1‑141 (MQ=255)                                                    aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATc                                                                                         <  1:305847/141‑1 (MQ=255)                                                    aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATc                                                                                         <  1:18041/141‑1 (MQ=255)                                                     aataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACAACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGAGGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCt                                                                                        >  2:76765/1‑141 (MQ=255)                                                      acaaaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                                                     >  1:138666/15‑141 (MQ=11)                                                      acaaaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                                                     >  1:133051/15‑141 (MQ=11)                                                         aaatTGTGATCTAATAAGAATAGAAGGTGCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCCTGTCAGTGTTCCTTTAACTTGCTAGCATGATGCTAATGTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                                                    <  2:113847/141‑1 (MQ=14)                                                           atagtattttaataaAAAT‑GATGGTTCCGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                <  1:247329/132‑1 (MQ=11)                                                           agagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                <  1:56848/132‑1 (MQ=11)                                                           agagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                <  1:247202/132‑1 (MQ=11)                                                            gagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                                               <  1:276739/133‑1 (MQ=11)                                                             tgtgATCTAATAAAAATCGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                <  2:106749/141‑1 (MQ=255)                                                               tgAGCTAATAAAAATAAAAGGTTCTTTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagaccaca                                                                              <  2:252766/141‑6 (MQ=255)                                                                 ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATGATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCGTGGCGAAAATCAGTAGATCTGAAGAGGCCTGCgg                                                                          >  1:133254/4‑141 (MQ=11)                                                                    taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACGTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                                         >  2:167419/1‑141 (MQ=255)                                                                     aataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTGGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                                                      <  2:53817/141‑1 (MQ=11)                                                                           aaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                                                                  <  1:133753/141‑1 (MQ=34)                                                                                         ttacAAAGTAAAAAAATATAGATAACCACTAATTTATCATGTCAGTGTTCGCTTAACTGGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTaaa                                                    <  2:170172/138‑1 (MQ=0)                                                                                                                    accaatatctCATGGTAGTGTTCGTTGAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCTAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAGAGTATTCt                         <  2:333292/132‑1 (MQ=1)                                                                                                                      caagataccATGTTAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCttt                       <  1:249583/133‑1 (MQ=1)                                                                                                                                    aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt         >  2:58900/1‑141 (MQ=11)                                                                                                                                    aGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt         >  2:2305/1‑141 (MQ=14)                                                                                                                                     gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTTGTGGTATGGTGAAAATCCGTATATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt        >  2:279205/1‑141 (MQ=11)                                                                                                                                     gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt        <  2:242828/141‑1 (MQ=255)                                                                                                                                     gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt        <  2:165168/141‑1 (MQ=12)                                                                                                                                     gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGCACAGACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt        <  1:175654/141‑1 (MQ=255)                                                                                                                                     gggtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt        <  1:15005/139‑1 (MQ=255)                                                                                                                                      tgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACttt       >  1:115136/1‑141 (MQ=11)                                                                                                                                      tgtTCGTCTAACTTGCTAGCATGATGCTAATTTCGGGGCATGGGTAAACTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTGTTTGCCGTATTGATACttt       >  2:235536/1‑141 (MQ=11)                                                                                                                                       gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTg      >  1:91551/1‑141 (MQ=17)                                                                                                                                       gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTg      >  1:218223/1‑141 (MQ=17)                                                                                                                                       gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACATTTTAAAGTATTCTTTGCTGTATTGATACtttt      >  1:277969/1‑140 (MQ=17)                                                                                                                                         tCGTTTAACGTGCTAGCATGATGCTAATTGCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt    <  1:254666/141‑1 (MQ=255)                                                                                                                                         tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt    >  1:34949/1‑141 (MQ=32)                                                                                                                                           gCTTGACTGGCTAACATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCGGGCGAGACGAGGGGTTCTTTTTAGATATGGCGTAAATACACTCCATACCTTATAAAGTATTCCTTGCTATCTTGATAcgtggagac  >  2:18553/1‑133 (MQ=2)                                                                                                                                            |                                                                                                                                           AGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAC  >  CP000730/35736‑36015 | 
| 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: ‑ |