| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | CP000730 | 76,889 | G→A | 22.8% | intergenic (+947/‑187) | USA300HOU_0072 → / → USA300HOU_0075 | universal stress protein/lysophospholipase | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 76,889 | 0 | G | A | 22.8% | 73.2 / 11.5 | 35 | intergenic (+947/‑187) | USA300HOU_0072/USA300HOU_0075 | universal stress protein/lysophospholipase | 
| Reads supporting (aligned to +/- strand): ref base G (12/15); new base A (2/6); total (14/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.31e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01 | |||||||||||
TATCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTA  >  CP000730/76751‑77027                                                                                                                                          |                                                                                                                                          tATCAAAATATCAATACGGCAAAGAATACTTTAAAAGGTATTGATTGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACg                                                                                                                                          <  2:262835/141‑1 (MQ=11)       atatCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTg                                                                                                                                   >  1:230292/1‑141 (MQ=11)              aacggCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTTTAAATGGTATGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACTCTGACATGAt                                                                                                                            <  1:288281/137‑1 (MQ=1)                               tacaaGGGATTGAATGTACTTACGCTATTTATAAAAAGAAACGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGACACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCta                                                                                                           <  1:119266/138‑1 (MQ=1)                                    ggtcttgggtggaTTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATAttt                                                                                                      >  2:344883/13‑141 (MQ=1)                                           agtgcaTTTACGCTCTATATAAAAAGAACCGCACGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATAttttttactt                                                                                               <  1:20127/136‑5 (MQ=1)                                                 ttgaCGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAAc                                                                                        >  2:111397/4‑141 (MQ=1)                                                        tcgatataAAAATAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATGAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcca                                                                                 <  2:225532/138‑2 (MQ=1)                                                        tccatctaAAAAGAACCGCAGGTCTCTTCACATCTTCGGATTTTGGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAAcc                                                                                  <  2:217049/135‑1 (MQ=1)                                                        cctacagaaaaaGAACCGTAGGTCTCTTCAGATCTACGGATTTTCCCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATA‑TTTTTTACTTTGCAACAGAAcca                                                                                 <  1:218019/134‑2 (MQ=1)                                                          tatataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCAGCCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCta                                                                                >  2:162075/1‑141 (MQ=11)                                                           atataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCtat                                                                               >  1:145919/1‑141 (MQ=31)                                                             ataAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATGAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaa                                                                             <  2:21387/141‑1 (MQ=17)                                                               aaaaaGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaa                                                                           <  1:237099/141‑1 (MQ=34)                                                               aaaaaGAACCGCAGGTCTCATCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaa                                                                           >  1:238324/1‑141 (MQ=21)                                                                aaaaGAAACGCAGCTCTTTTCAGATCTACGGATGTTCGCTATGCCACGAAATTAGCACCATGATAGCAAGTTAAGCGAACACTGACATTATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATaaaaa                                                                          <  2:86919/141‑1 (MQ=11)                                                                      aCCGCAGGCCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATgg                                                                    <  2:347105/141‑1 (MQ=255)                                                                        cGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGAGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAACCATGGGc                                                                  >  1:303927/1‑141 (MQ=255)                                                                         gCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAAAACCTACAAAAAACATGGGCt                                                                 <  1:213476/141‑1 (MQ=255)                                                                         gCAGGTCTCTTCAGATCTACGGACTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCt                                                                 <  2:107346/141‑1 (MQ=255)                                                                             gTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTaaa                                                             >  2:115110/1‑141 (MQ=255)                                                                                  ttCAGATCTACGGATTTCTGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGtt                                                        <  2:18655/141‑1 (MQ=255)                                                                                  ctCAGAGCTACGGATTTTCGCCATGCCACGCAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATCTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGtt                                                        <  2:256252/140‑1 (MQ=255)                                                                                           gcgtgttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatat                                               <  2:195641/136‑1 (MQ=255)                                                                                           aCGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatat                                               >  1:122281/1‑141 (MQ=255)                                                                                            cGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTtatata                                              <  1:163463/141‑1 (MQ=255)                                                                                               aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTa                                           >  1:171136/1‑141 (MQ=255)                                                                                               aTTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTa                                           >  2:321990/1‑141 (MQ=255)                                                                                                ttttCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTAAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTAt                                          <  2:50317/141‑1 (MQ=255)                                                                                                    cGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACt                                      >  2:398822/1‑141 (MQ=255)                                                                                                             aCGAAAGTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACAtt                             <  1:96047/141‑1 (MQ=255)                                                                                                                        atcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCa                  >  1:334940/1‑141 (MQ=255)                                                                                                                         tcatGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACAGTGATTTGGGTCaa                 <  2:279530/141‑1 (MQ=255)                                                                                                                                       tAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATACAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGAGGTAtt   >  2:1362/1‑141 (MQ=255)                                                                                                                                        aaGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTa  <  2:122281/141‑1 (MQ=255)                                                                                                                                          |                                                                                                                                          TATCAAAATATCAATACAGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGCAGGTCTCTTCAGATCTACGGATTTTCGCCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAGCGAACACTGACATGATAAATTAGTGGTTAGCTATATTTTTTTACTTTGCAACAGAACCTATAAAAAACATGGGCTTAAAGGGTTAATTTATATAGTATTACTGCAAACATTGATTTGGGTCAAAAATTAGGGGTATTA  >  CP000730/76751‑77027 | 
| 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: ‑ |