Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,180 C→T 16.7% 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,1800CT16.7% 90.8 / 6.3 35intergenic (+238/‑896)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (14/14);  new base T (3/3);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTTAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTA  >  CP000730/76042‑76319
                                                                                                                                          |                                                                                                                                           
aaTTTAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCta                                                                                                                                           >  2:149149/1‑141 (MQ=255)
     aaTATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGta                                                                                                                                      >  1:67789/1‑141 (MQ=255)
              aaTATTGCTGAACCTGTATCACATAGAACAGCTAACCAAAGTGATAGAAAACCAGAAAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAgaaaaaaaaaaa                                                                                                                             <  2:251594/141‑12 (MQ=255)
                   tgctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaa                                                                                                                        <  1:170302/141‑1 (MQ=255)
                          ccTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtc                                                                                                                 <  1:275284/141‑1 (MQ=255)
                                tCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGt                                                                                                           >  1:13582/1‑141 (MQ=255)
                                  aCTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTAt                                                                                                         <  2:238773/141‑1 (MQ=255)
                                    ttAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGa                                                                                                       <  1:202462/141‑1 (MQ=255)
                                      aGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAAc                                                                                                     <  1:303214/141‑1 (MQ=255)
                                             cgaattaaaGTGGTAGTAATCCAGGAAATATTAACACAAAGGGAATTAATTTAAGAATAATATCAAAATAGACGTTATGATGCAAAGTAAAAGCTATAGTAGAATTTTAACAAAAAGGAGTCTTCTGTTGGAACGAtaaca                                                                                              >  2:64639/7‑137 (MQ=19)
                                                             tAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGATCTGTTGCAAAGTTGAATCTATAGGATAATATGAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTaa                                                                              >  2:237731/1‑141 (MQ=255)
                                                                tCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACACTTTACCaa                                                                           >  2:267254/1‑141 (MQ=255)
                                                                     ggattACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATg                                                                      >  2:55104/5‑141 (MQ=255)
                                                                     agAATATTAAAACAAAATCAATTAATTTAATAATTATAACAAAATATAGGATATGATTAAAAGATAAATATATAATATTATTATAAAAAAAAGGAATCTAATGTATGAATTATTTGAGATATAAACAATTTataaaagaat                                                                      >  2:64029/2‑132 (MQ=255)
                                                                       aaTACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGtt                                                                    <  1:60946/141‑1 (MQ=255)
                                                                           cTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACaaagcgagtatca                                                                >  2:178301/1‑130 (MQ=255)
                                                                                     ggCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGt                                                      >  2:90283/1‑141 (MQ=255)
                                                                                       cAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTg                                                    >  2:221785/1‑141 (MQ=255)
                                                                                        aattaattTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATATAAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTgg                                                   <  2:152717/141‑1 (MQ=255)
                                                                                          ttaatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGccaacggct                                                 <  1:58559/141‑8 (MQ=255)
                                                                                          ttaatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGct                                                 <  1:240944/141‑1 (MQ=255)
                                                                                              tttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctacta                                             >  1:13762/1‑141 (MQ=255)
                                                                                                  ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTa                                         >  1:127288/1‑141 (MQ=255)
                                                                                                  ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTAAAAATAAAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTa                                         <  2:60688/141‑1 (MQ=255)
                                                                                                    aatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAg                                       <  1:176929/141‑1 (MQ=255)
                                                                                                      ttATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGat                                     <  1:70709/141‑1 (MQ=255)
                                                                                                           ataaaaTATGGTTCTGTTGAAGAGTTGAATTTATGATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATTTAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTTTATATGATATGCAt                               <  2:59234/139‑1 (MQ=1)
                                                                                                             agttaggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCAtt                              >  1:327674/8‑141 (MQ=1)
                                                                                                                 atcaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt                          <  2:286777/138‑1 (MQ=9)
                                                                                                                 aggaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgc                          <  2:331812/138‑2 (MQ=2)
                                                                                                                    aGGTTCTGTTGCAAAGTTGTATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGACGTTGGCTACTATATAAGATATGCATTGAGCTAt                       >  2:302387/1‑141 (MQ=11)
                                                                                                                    aGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTAt                       >  1:188837/1‑141 (MQ=2)
                                                                                                                               cAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATc            >  2:226681/1‑141 (MQ=35)
                                                                                                                                  aGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGa         <  2:95958/141‑1 (MQ=35)
                                                                                                                                         tCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTa  >  2:192145/1‑141 (MQ=34)
                                                                                                                                          |                                                                                                                                           
AATTTAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTA  >  CP000730/76042‑76319

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: