| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | CP000730 | 1,929,925 | T→A | 55.9% | intergenic (‑1008/‑160) | USA300HOU_1794 ← / → USA300HOU_1796 | hypothetical protein/hypothetical protein | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 1,929,925 | 0 | T | A | 55.9% | 6.5 / 50.1 | 46 | intergenic (‑1008/‑160) | USA300HOU_1794/USA300HOU_1796 | hypothetical protein/hypothetical protein | 
| Reads supporting (aligned to +/- strand): ref base T (10/10); new base A (16/9); total (27/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.79e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01 | |||||||||||
TATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTA  >  CP000730/1929792‑1930057                                                                                                                                         |                                                                                                                                    tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                                                                                                              >  2:79705/9‑141 (MQ=255) tttcggcGATGATATATTAGTAGAGGCGCCTATCTGTCAGGCGTCAACTTAGACGCAGAGAGGATGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCGTAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                                                                                                              <  1:85934/135‑1 (MQ=255) atTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCTTTTTTTTCtt                                                                                                                              <  1:288682/141‑1 (MQ=255)   ttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAt                                                                                                                            <  2:18332/141‑1 (MQ=255)   tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                           >  1:37368/6‑141 (MQ=255)      cagacggtataATTATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCAATTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                        <  1:87477/134‑1 (MQ=255)       cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGCGAGGAGGTGTA‑TAA‑G‑GC‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                        <  1:260833/141‑1 (MQ=255)         agcatatAATAGGAGAGTCGCCTATCTCTCAGGCGTCAATTTAGCCGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa                                                                                                                      <  2:215536/138‑1 (MQ=255)          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCCTTATTGGCTTAg                                                                                                                    >  1:220004/4‑141 (MQ=255)          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                    <  1:141444/138‑4 (MQ=255)          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTACTCATTTTCGTTCACGTCATAGCACCTGTCATCAGTGGCTGTGCTGCTGCGTATTATACTTATTGGCTTAg                                                                                                                    <  2:275815/138‑1 (MQ=255)          tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTaa                                                                                                                     >  2:238635/1‑141 (MQ=255)           ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAgg                                                                                                                   >  2:291720/3‑140 (MQ=255)           ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                   >  2:319625/3‑141 (MQ=255)              ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                >  1:313259/1‑139 (MQ=255)               taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAAGATCGGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                               <  1:30533/141‑4 (MQ=255)               gaataatAGGGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGCGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTCTTGGCTTAGTaaac                                                                                                               <  2:272455/140‑4 (MQ=255)                 ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTATTAAACGc                                                                                                             >  2:159873/5‑141 (MQ=255)                 ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                              >  1:335014/1‑141 (MQ=255)                 ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCCCTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                              >  1:290362/1‑141 (MQ=255)                  taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                            <  1:297687/138‑1 (MQ=255)                  caatatagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                            <  2:196202/135‑1 (MQ=255)                    atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GAGGAAGAGAGGAGGTGCATTAC‑T‑GT‑GCTAATCATTTTCGTTCGCATCATAGCACCAGTCATCAGTCGCTGTGCTGTGGCGTCTTATACTTATTGGCTTAGTCAACGCaac                                                                                                          >  2:301096/2‑140 (MQ=255)                         gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                                                                      <  1:24325/141‑1 (MQ=255)                         gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacgaata                                                                                                     >  1:12308/1‑135 (MQ=255)                         gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                     >  2:286085/1‑135 (MQ=255)                          tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                    >  2:314126/1‑134 (MQ=255)                            gcctttccgtcAGGCTTCAATTTAAGCGCAGAGGGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTGTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                                                                                   <  2:77131/132‑1 (MQ=255)                            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAATCATCAGTGGCTTTGCTGTTGCGTATCATACTTATTGGCTTAGTAAACGCaacaaatgatt                                                                                                  >  2:102184/1‑132 (MQ=255)                            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                  >  1:21717/1‑132 (MQ=255)                            gCCTATCTCGCAGGCGTCAAGT‑‑GGCGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTGCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                  <  1:76933/141‑10 (MQ=255)                              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                <  1:312343/141‑12 (MQ=255)                              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTCAATGGCTTAGTAAACGCaacaaataattaa                                                                                                >  1:187558/1‑130 (MQ=255)                               tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTAATTGGCTTAGTAAACGCaacaaataattaaa                                                                                               >  2:105831/1‑129 (MQ=255)                               tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGCAGGTGCA‑TTA‑C‑TTGGCTAATCATTTTCACTCACAACATAGCACAAGTCATGAGTGGCTGTGCTGATGCGTATTATACTTATTGGCTTAGTAAACGCaacaactaattaat                                                                                               >  2:149117/1‑129 (MQ=255)                                 tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                             >  2:295300/1‑127 (MQ=255)                                     tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATaacatacct                                                                                          >  2:278132/1‑133 (MQ=255)                                     tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACat                                                                                          <  1:93331/141‑1 (MQ=255)                                             cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc                                                                                  >  2:352045/1‑141 (MQ=255)                                                               gaggTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAGACAAAATAGGTGCCATATAGCCGCACCAATAAAAATcccccc                                                                >  2:28321/1‑139 (MQ=255)                                                                                     cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCCCACTACCGCAAATAGTGAg                                              <  1:211986/141‑1 (MQ=255)                                                                                      ttATATTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGACTAAGTAGACGCAATAAAAAATAGGTGAAATATAGCCGCACCAATAAAAATCACCTCACTACAGAAAATAATGAgg                                             >  1:167631/1‑141 (MQ=255)                                                                                                       atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTa                            >  1:102260/1‑141 (MQ=255)                                                                                                                  catcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAACTAGTGAGGGGCTTGGTGTATAAGTAAATACTAAttt                 >  2:152766/1‑141 (MQ=255)                                                                                                                     caGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAGTGGTGTATAAGTAAACACTTATTTTCg              >  2:156389/1‑141 (MQ=255)                                                                                                                                 cagtGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATta  <  1:51878/138‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    TATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTA  >  CP000730/1929792‑1930057 | 
| 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: ‑ |