Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,180 C→T 32.5% intergenic (+238/‑896) 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,1800CT32.5% 78.5 / 33.2 49intergenic (+238/‑896)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (18/15);  new base T (5/11);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AATTTAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATAT  >  CP000730/76042‑76317
                                                                                                                                          |                                                                                                                                         
aaTTTAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTCCAAAGTTGAACCta                                                                                                                                         >  2:164980/1‑141 (MQ=255)
aaTTTAAAATAACTAATACTGCTGCACCTGTATGACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCta                                                                                                                                         <  2:75496/141‑1 (MQ=255)
     aaTATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGta                                                                                                                                    >  1:215712/1‑141 (MQ=255)
          aaCTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAAtt                                                                                                                               >  2:325247/1‑141 (MQ=255)
            cTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAAtttt                                                                                                                             <  2:4214/141‑1 (MQ=255)
                tatTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACa                                                                                                                         >  2:215345/1‑141 (MQ=255)
                    gctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACCAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaaa                                                                                                                     >  1:153910/1‑141 (MQ=255)
                      tgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAgg                                                                                                                   >  1:132527/1‑141 (MQ=255)
                      tgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAgg                                                                                                                   >  1:85266/1‑141 (MQ=255)
                      tgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAgg                                                                                                                   <  1:292835/141‑1 (MQ=255)
                          ccTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtc                                                                                                               <  1:269049/141‑1 (MQ=255)
                                tCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGt                                                                                                         <  2:215712/141‑1 (MQ=255)
                                                     aGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATaaa                                                                                    >  2:299408/1‑141 (MQ=255)
                                                           aGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAAttt                                                                              >  1:308947/1‑141 (MQ=255)
                                                             tAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTaa                                                                            >  1:250888/1‑141 (MQ=255)
                                                                  cAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAgg                                                                       <  1:215345/141‑1 (MQ=255)
                                                                    gggAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGAt                                                                     >  1:277265/1‑141 (MQ=255)
                                                                         tACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCCGATATAAACAATTTAACAAGGATGTTAt                                                                >  1:231888/1‑141 (MQ=255)
                                                                               cacaAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGt                                                          <  1:220342/141‑1 (MQ=255)
                                                                               cacaAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGt                                                          <  1:325247/141‑1 (MQ=255)
                                                                                   aaGGCAATTAATTTAATAATTATCGCAAAATATAGGATCTGTGCCAAAGATGAATCAAGAGGATAATTATAACAAAAAGGAGTCTTCTGTATGAACTATTTAAGATATAAACAATTTAACAAGGATGTTATCAATGTAGcc                                                      >  2:37600/1‑141 (MQ=18)
                                                                                           aaaatgaataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTAAAAATAAAAAAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGcta                                              <  2:211457/135‑1 (MQ=17)
                                                                                             atttaatAATTATAGCAAAATATAGGTTCTGTTGTAAAGTTGAAGCTATAGTATAATTTTAACAAAAATGAGTGTTTTGTTTGTTTTCTAACGGTTAAAAACACTTTAACAAGGATGTTATCACTGTTGTCGTTGGctact                                            <  2:321297/141‑1 (MQ=12)
                                                                                                            cAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCTGTTATTCACTATTTAACAAGTTTGTTACAGCTTTAGCCGTTGGCTACTATCTAAGATATGCAt                             <  2:321290/141‑1 (MQ=14)
                                                                                                            cAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCAt                             <  1:25213/141‑1 (MQ=255)
                                                                                                            agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTTTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGTTGTTGGCTACTATCTAAGATATGCAt                             >  1:203630/9‑141 (MQ=1)
                                                                                                            agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                             >  2:226687/9‑141 (MQ=1)
                                                                                                            agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                             <  1:292981/133‑1 (MQ=11)
                                                                                                             gtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAtt                            >  1:282624/8‑141 (MQ=11)
                                                                                                             agttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCAtt                            <  1:48079/134‑1 (MQ=255)
                                                                                                              taaatgaGGTTCTGGTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTg                           <  2:841/135‑1 (MQ=11)
                                                                                                              taaatgaGGGTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTg                           <  2:127806/135‑1 (MQ=11)
                                                                                                               aaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                          <  1:3634/136‑1 (MQ=9)
                                                                                                                aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACACGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAttagt                         >  1:21769/5‑138 (MQ=11)
                                                                                                                aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                         <  2:20968/137‑1 (MQ=9)
                                                                                                                 tataGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACTATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGt                        <  1:128412/141‑1 (MQ=255)
                                                                                                                 aggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgc                        <  1:124480/138‑2 (MQ=2)
                                                                                                                  ggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgct                       <  1:2280/139‑3 (MQ=9)
                                                                                                                  ataGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGtt                       >  1:225046/1‑141 (MQ=38)
                                                                                                                   taGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACACGGATGTTTTCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTa                      >  1:274754/1‑141 (MQ=255)
                                                                                                                    aGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATTTAAGATATGCATTGAGTTAt                     >  2:75770/1‑141 (MQ=0)
                                                                                                                     ggTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATc                    >  1:200037/1‑141 (MQ=34)
                                                                                                                       ttcgttttcAAAGTTGAATTTATAATATAATTTGAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGt                  <  2:327090/133‑1 (MQ=11)
                                                                                                                        tgtGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGGATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTg                 <  2:294745/139‑1 (MQ=255)
                                                                                                                            ttGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatat             <  1:33090/141‑1 (MQ=35)
                                                                                                                             tggaaaGTTGAATTTATAACACAATTTTAACAAAAAGGAGTCTTCTGTATGAACGATTTCAGATATAAACAATTTAGCAAGGATGTTAGTACGGTAGCCGTTGGCTACTATCGAAGATATGCAGTGAGATATCGTGatata            <  2:109205/138‑1 (MQ=1)
                                                                                                                                  aGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGa       <  2:197066/141‑1 (MQ=35)
                                                                                                                                     tGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGATACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAat    >  2:261204/1‑141 (MQ=21)
                                                                                                                                       aaTCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAatat  >  2:320067/1‑141 (MQ=35)
                                                                                                                                          |                                                                                                                                         
AATTTAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATAT  >  CP000730/76042‑76317

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