Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,185 G→A 9.7% 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,1850GA 9.7% 114.5 / 4.1 54intergenic (+243/‑891)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base A (2/3);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

TGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGA  >  CP000730/76061‑76324
                                                                                                                            |                                                                                                                                           
tgctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaa                                                                                                                             >  2:232084/1‑141 (MQ=255)
tgctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaa                                                                                                                             >  2:210103/1‑141 (MQ=255)
      aCCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGt                                                                                                                       <  1:142585/141‑1 (MQ=255)
       ccTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtc                                                                                                                      >  1:104872/1‑141 (MQ=255)
       ccTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtc                                                                                                                      >  1:104167/1‑141 (MQ=255)
         tGTATCAATTAGAACAGAAAAAAAAAGAGAAAGTAATCCATGTAATACTAACACAAAGGAAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAATGTTGTATCTATAGTAAAAAAAAAAAAAAAAGGAGtctt                                                                                                                    <  2:109531/141‑1 (MQ=255)
                   aGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAAc                                                                                                          <  1:266315/141‑1 (MQ=255)
                     aaCAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTa                                                                                                        <  2:294881/141‑1 (MQ=255)
                         gCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTc                                                                                                    <  1:264321/141‑1 (MQ=255)
                                aaaCTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGAtata                                                                                             <  2:181844/141‑1 (MQ=255)
                                     gtTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACaa                                                                                        <  1:177789/141‑1 (MQ=255)
                                         gTAATCCAGGGAATACGAACACAAAGGCAATTTATTTAATAATTAGAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATGATAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTa                                                                                    <  2:275417/141‑1 (MQ=21)
                                                   gAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGt                                                                          >  1:272632/1‑141 (MQ=255)
                                                    aaTACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGtt                                                                         >  2:197693/1‑141 (MQ=255)
                                                    aaTACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGtt                                                                         <  2:248620/141‑1 (MQ=255)
                                                          aaCACAAACGCAATTAATTTAATTATTAACGCAATATATAGGTTATGTTGAAAAGTATAATCTAAAGTATAATTTTAACTAAAAGGAGTCTTCCGTATGCACTAGTTCAGATATAAACAATTTAACAAGAATGTAAtacct                                                                   >  2:150105/1‑137 (MQ=17)
                                                           acacAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTAAACAAAAAGGAGTCTTCTGTATGAACTATTTCCGATATAAACAATTTAACAAGGATGTTATCAcgg                                                                  <  2:217671/141‑3 (MQ=255)
                                                             acaAAGGCAATTAATATAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTa                                                                <  1:233074/141‑1 (MQ=255)
                                                              caAAGGCAATTAATTTAATAATTATAGCAAAATATAGATTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAg                                                               <  1:288107/141‑1 (MQ=255)
                                                                    cAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTg                                                         >  2:103729/1‑141 (MQ=255)
                                                                     aatgaagtTAACAACGATAGAAAAAAAAAGGTTCTTTTGCAAAGTTGAATCTATAGTAAAAAAAAAAAAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTgg                                                        <  2:196424/141‑1 (MQ=11)
                                                                           tttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCCATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctacta                                                  >  1:116791/1‑141 (MQ=255)
                                                                              aataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCt                                               >  2:335496/1‑141 (MQ=255)
                                                                               ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTa                                              >  1:83984/1‑141 (MQ=255)
                                                                                          gtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAtt                                   <  2:237635/134‑1 (MQ=255)
                                                                                          gtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGCCTACCATCTAAGATATGCAtt                                   <  1:191025/134‑1 (MQ=11)
                                                                                          gtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATAGAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAtt                                   <  1:51004/134‑1 (MQ=11)
                                                                                           gttagaaGGTTCTTTTGCAAAGTTGAATTTATAATATAATTAAAAACAAAAGTAATCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTg                                  <  2:145906/135‑1 (MQ=11)
                                                                                           aaaatgaGGTTCTGTTGCAAAGTTGAATGTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTg                                  <  2:165349/135‑1 (MQ=255)
                                                                                           aaaTATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTg                                  <  1:330115/141‑1 (MQ=255)
                                                                                            aaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                                 >  1:305998/6‑141 (MQ=9)
                                                                                            aaTAGAGGTTCTGGTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGCGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGa                                 <  2:262436/141‑1 (MQ=255)
                                                                                             aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                                <  2:171997/137‑1 (MQ=2)
                                                                                             aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                                <  1:333272/137‑1 (MQ=2)
                                                                                              aggaGGTTCTTTTGCAAAGGTGAATTTATAATATAATTATAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgc                               <  2:365222/138‑2 (MQ=11)
                                                                                               tgaGGTTCTGTTGCAAAGTTGCATTTATAGTATAATTTTAAACAATAGAAGAGTTCCGTACAAACTATTTCAGATATAAACATTTTAACAATGATTTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGtt                              <  2:359269/139‑1 (MQ=11)
                                                                                               ataGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAAAAAGGAGGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTgcggt                              >  1:134857/1‑137 (MQ=255)
                                                                                                gaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGCGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGttt                             >  2:345130/2‑140 (MQ=2)
                                                                                                gaGGTTCTGTTGCAAAGTTGAATTTATAATATAATCTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgctc                             >  2:66066/2‑138 (MQ=255)
                                                                                                 aGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATTTTGGGCATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTAt                            >  1:271527/1‑141 (MQ=1)
                                                                                                 aGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTAt                            >  1:325511/1‑141 (MQ=0)
                                                                                                         ttGCAAAGTTGAACGTATATTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTCGCCGTTGGCTACTATCTAAGAGATGCATAGCGCTATCGTGataa                    >  2:322340/1‑140 (MQ=11)
                                                                                                          tGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatata                   >  2:340717/1‑141 (MQ=35)
                                                                                                          tGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatata                   >  2:300244/1‑141 (MQ=35)
                                                                                                            gaaaGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTTTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATc                 <  2:112078/140‑1 (MQ=255)
                                                                                                            cAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATc                 >  1:128590/1‑141 (MQ=35)
                                                                                                               aGTTGAATGTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGACGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGa              <  1:51454/141‑1 (MQ=12)
                                                                                                                agtGAATCTATAGTAAAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTGGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGaa             <  1:290331/139‑1 (MQ=255)
                                                                                                                   tgttCTGTAGTATAATTTTAACAAAAAGGAGGCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTGTCACTGTAGCCGTTGGCTACTATCTAAGATATGCGTTGAGTTATCGTGATATATCTGAAata          >  2:3970/4‑141 (MQ=11)
                                                                                                                     atttataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTCACACATGAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATAtt        >  1:367714/4‑141 (MQ=11)
                                                                                                                     atttataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATAtt        >  2:27225/4‑141 (MQ=0)
                                                                                                                     atttataCTCTAACTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATAtt        >  2:266406/4‑141 (MQ=1)
                                                                                                                      tttataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTa       <  2:21677/139‑1 (MQ=9)
                                                                                                                      tttataGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTa       >  1:207323/3‑141 (MQ=14)
                                                                                                                      tttataGTAAGATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGCGTTATCGTGATATATCTGAAatataa       >  1:273653/3‑139 (MQ=11)
                                                                                                                       ttataGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTaa      <  1:171819/140‑1 (MQ=17)
                                                                                                                        tataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCATATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAg     <  1:122693/141‑1 (MQ=12)
                                                                                                                        tataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAg     <  1:25807/141‑1 (MQ=12)
                                                                                                                           aGTATAATTTTAACAAAAAGGAGTTTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGa  >  1:237521/1‑141 (MQ=17)
                                                                                                                            |                                                                                                                                           
TGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGA  >  CP000730/76061‑76324

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