| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | CP000730 | 1,460,315 | T→C | 13.5% | intergenic (‑347/+53) | USA300HOU_1356 ← / ← ctpA | possible acetyltransferase/S41A family carboxy‑terminal peptidase | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 1,460,315 | 0 | T | C | 13.5% | 67.0 / 5.9 | 30 | intergenic (‑347/+53) | USA300HOU_1356/ctpA | possible acetyltransferase/S41A family carboxy‑terminal peptidase | 
| Reads supporting (aligned to +/- strand): ref base T (13/12); new base C (2/2); total (16/14) | |||||||||||
| 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.20e-01 | |||||||||||
TGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTAAATTTATTATTCGTTTCTTTA  >  CP000730/1460168‑1460461                                                                                                                                                   |                                                                                                                                                  tGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAGTTTa                                                                                                                                                   >  2:30537/1‑149 (MQ=255)tGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTAGCTTTGCAAATTTa                                                                                                                                                   >  1:290512/1‑149 (MQ=255) gCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTag                                                                                                                                                  <  1:162944/149‑1 (MQ=255)  cTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTagg                                                                                                                                                 >  1:221152/1‑149 (MQ=255)  cTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTagg                                                                                                                                                 >  2:201034/1‑149 (MQ=255)       tGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCg                                                                                                                                            >  1:9423/1‑149 (MQ=255)       tGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCg                                                                                                                                            >  1:168583/1‑149 (MQ=255)          atatTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCgtgt                                                                                                                                         >  1:147127/1‑149 (MQ=255)                    aTGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACata                                                                                                                               <  1:235409/149‑1 (MQ=255)                       aaTTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACataata                                                                                                                            >  1:6746/1‑149 (MQ=255)                        aTTGTTATCAGCGATGATAGAAATTGCGTTATTTACATAATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACataataa                                                                                                                           >  1:85055/1‑149 (MQ=255)                           gTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACataataataa                                                                                                                        >  2:48555/1‑149 (MQ=255)                                   cGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAAttt                                                                                                                >  2:66673/1‑149 (MQ=255)                                    gatgatAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAAtttt                                                                                                               <  1:56179/149‑1 (MQ=255)                                            aaaTTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAgtgt                                                                                                       <  2:9423/149‑1 (MQ=255)                                            aaaTTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTCAGGAGGCGTGTCATTAACATAATAATAATACAAATTTAGTAgtgt                                                                                                       <  2:168583/149‑1 (MQ=255)                                                        ttACATTATTTGTAAATGTCTAATTTATATCTTTTAGAGTTATACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTAtt                                                                                           <  2:288991/149‑1 (MQ=255)                                                             tcatTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTaaa                                                                                      >  1:27860/1‑149 (MQ=255)                                                                       aTGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCa                                                                            <  1:201034/149‑1 (MQ=255)                                                                                ttATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCga                                                                   <  1:32911/149‑1 (MQ=255)                                                                                    atCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAAc                                                               >  2:216511/1‑149 (MQ=255)                                                                                            gagTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAAcatcatcat                                                       >  2:239611/1‑149 (MQ=255)                                                                                            gagTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAAcatcatcat                                                       >  1:187033/1‑149 (MQ=255)                                                                                                 aGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTtta                                                  <  2:145290/149‑1 (MQ=255)                                                                                                      tCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGc                                             <  2:211313/149‑1 (MQ=255)                                                                                                          taaaaataaaCACGCCTCCTAAATTGCTATCTTTGCAAATTGAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCtttt                                         >  1:157093/1‑149 (MQ=255)                                                                                                           aaaaataaaCACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCttttt                                        <  2:2252/149‑1 (MQ=255)                                                                                                                           cctAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGt                        <  2:159953/149‑1 (MQ=255)                                                                                                                              aaaTTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTaaa                     <  2:264622/149‑1 (MQ=255)                                                                                                                              aaaTTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATAGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTaaa                     <  1:108746/149‑1 (MQ=255)                                                                                                                                                 tttAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTAAATTTATTATTCGTTTCTtta  <  2:6746/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TGCTCATTGCATATTAAACCATGAATTGTTATCAGCGATGATAGAAATTGCGTTATTTACATCATTTGTAAATGTCTAATTTATATCTTTTAGAGTTAGACGTCAATAAAAATAAACACGCCTCCTAAATTGCTATCTTTGCAAATTTAGGAGGCGTGTCATTAACATAATAATAATACAATTTTAGTAGTGTGTATCGCTTATTTTAAAATATTAATCAACTTATCGAGAACATCATCATGTTTATTAGCTTTTTCAACTAATAACTCAGTAAATTTATTATTCGTTTCTTTA  >  CP000730/1460168‑1460461 | 
| Alignment Legend | 
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG | 
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑ |