Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,185 G→A 21.6% intergenic (+243/‑891) USA300HOU_0072 → / → USA300HOU_0075 universal stress protein/lysophospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073076,1850GA21.6% 98.5 / 25.5 51intergenic (+243/‑891)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base A (8/3);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01

AATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGG  >  CP000730/76056‑76330
                                                                                                                                 |                                                                                                                                                 
aaTATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAgg                                                                                                                                <  2:84798/149‑1 (MQ=255)
                 aTCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACt                                                                                                               >  1:90668/1‑149 (MQ=255)
                          aaCAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGat                                                                                                      <  1:295054/149‑1 (MQ=255)
                              gCTAACCAAAGTGTTAGTAATCCATTGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATaa                                                                                                  >  2:62713/1‑149 (MQ=255)
                              acTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATaa                                                                                                  >  2:272633/2‑149 (MQ=255)
                                    cAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAAtt                                                                                            >  2:167950/1‑149 (MQ=255)
                                     aaaGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAAttt                                                                                           >  1:201734/1‑149 (MQ=255)
                                     aaaGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAAttt                                                                                           >  1:20898/1‑149 (MQ=255)
                                               tAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATg                                                                                 >  2:8438/1‑149 (MQ=255)
                                               tAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTAAAAAAAAAAAAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATg                                                                                 <  2:200636/149‑1 (MQ=25)
                                                                acacAAAAGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGtt                                                                <  2:153031/149‑1 (MQ=255)
                                                                 cacaAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTg                                                               >  1:77709/1‑149 (MQ=255)
                                                                   caAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTCATTCCGTTTTTCACGATATAACAAGGATGTTATCACTGTAGCCGTTGGc                                                             <  2:13491/149‑1 (MQ=18)
                                                                        gCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctacta                                                        <  1:117267/149‑1 (MQ=255)
                                                                        gCAATTAATTGAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctacta                                                        <  2:58719/149‑1 (MQ=255)
                                                                              aatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAg                                                  <  1:106837/149‑1 (MQ=255)
                                                                                  taataaTTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGatat                                              <  2:179371/149‑1 (MQ=255)
                                                                                     taatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCa                                           >  2:72007/1‑149 (MQ=255)
                                                                                         taatacagttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGa                                       >  1:193949/14‑149 (MQ=11)
                                                                                         catatagtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATGATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                                       >  2:275308/14‑149 (MQ=11)
                                                                                         catatagtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                                       >  1:189994/14‑149 (MQ=11)
                                                                                         catatagtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGGCTTCTGTATGAACTATTTCAGATATAAAAAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                                       >  1:144028/14‑149 (MQ=11)
                                                                                           taGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTATAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGt                                     <  2:7994/149‑1 (MQ=255)
                                                                                            tacagttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgct                                    >  1:199178/11‑147 (MQ=11)
                                                                                             acagttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgcta                                   >  1:203899/10‑146 (MQ=11)
                                                                                              cAAAATCTAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATCTGCATTGAGTTAt                                  >  2:27675/1‑149 (MQ=255)
                                                                                               gtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATc                                 >  1:224009/8‑149 (MQ=11)
                                                                                               gtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATc                                 >  1:85286/8‑149 (MQ=11)
                                                                                               gtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATc                                 >  1:91290/8‑149 (MQ=255)
                                                                                               gtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGGCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATc                                 >  2:233969/8‑149 (MQ=11)
                                                                                               agttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATc                                 >  2:117019/8‑149 (MQ=11)
                                                                                               agttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATc                                 >  2:116913/8‑149 (MQ=11)
                                                                                                taaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCg                                >  2:34576/7‑149 (MQ=2)
                                                                                                taaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCg                                >  2:299398/7‑149 (MQ=2)
                                                                                                gttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCg                                >  2:82916/7‑149 (MQ=11)
                                                                                                gttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCg                                >  2:286654/7‑149 (MQ=11)
                                                                                                 aaTATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGt                               <  2:18063/149‑1 (MQ=40)
                                                                                                  taggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTg                              >  2:174450/5‑149 (MQ=11)
                                                                                                  taggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTg                              <  1:51374/145‑1 (MQ=11)
                                                                                                  taggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTg                              <  2:125815/145‑1 (MQ=11)
                                                                                                   atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGa                             <  2:248476/146‑1 (MQ=2)
                                                                                                   atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGa                             <  1:264488/146‑1 (MQ=2)
                                                                                                    ggaGGTTCTGTTGCAAAGTTGAATTTATAATAAAAAAAAAAAAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGat                            <  2:45794/147‑1 (MQ=11)
                                                                                                    ataGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGat                            <  1:93697/149‑1 (MQ=38)
                                                                                                     gaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGata                           <  1:87044/148‑1 (MQ=2)
                                                                                                         ttCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATc                       >  1:94229/1‑149 (MQ=35)
                                                                                                             gTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGaa                   <  2:85820/149‑1 (MQ=35)
                                                                                                                gCAAAGTGGAATCTCTAGTAGAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAata                <  1:244033/149‑1 (MQ=255)
                                                                                                                 cAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAatat               <  1:123223/149‑1 (MQ=35)
                                                                                                                             tataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTTTTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACgtg   <  1:239693/149‑1 (MQ=11)
                                                                                                                             tataGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtg   >  2:298579/1‑149 (MQ=32)
                                                                                                                              ataGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtgg  <  2:90702/149‑1 (MQ=17)
                                                                                                                                 |                                                                                                                                                 
AATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGG  >  CP000730/76056‑76330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: