Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,286,175 C→T 6.3% Q469* (CAA→TAA)  sprA → serine‑type phage integrase; phage SPbeta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,286,1750CT 6.3% 208.5 / 3.1 63Q469* (CAA→TAA) sprAserine‑type phage integrase; phage SPbeta
Reads supporting (aligned to +/- strand):  ref base C (22/37);  new base T (2/2);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GAGAACTTCATCTTCGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAGCTATTAAGGATGGTCTTAAAAGATGTTATTGT  >  BS168_reference_genome/2286040‑2286299
                                                                                                                                       |                                                                                                                            
gagaACTTCATCTTCGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAAt                                                                                                                           <  1:268290/139‑1 (MQ=255)
gagaACTTCATCTTCGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAAt                                                                                                                           <  1:412169/139‑1 (MQ=255)
gagaACTTCATCTTCGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAAt                                                                                                                           <  1:65440/139‑1 (MQ=255)
gagaACTTCATCTTCGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAAt                                                                                                                           <  7:244956/139‑1 (MQ=255)
 agaACTTCATCTTCGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAAt                                                                                                                           >  7:34154/1‑138 (MQ=255)
     cTTCATCTTCGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGa                                                                                                                       >  4:425247/1‑138 (MQ=255)
      ttCATCTTCGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGaaa                                                                                                                     >  4:299833/1‑139 (MQ=255)
              cGATATATATGAATGTGGCATTTACCGTGATGAACTTTTTTTGAAGAGTAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACa                                                                                                               <  1:426431/137‑1 (MQ=255)
              cGATAAATATGAAAGTGTCATTTACAGTGATGAACTTTTTTTTAAGAGAAAAGCTGCATTAGAAGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGa                                                                                                                 <  8:222784/135‑1 (MQ=255)
              cGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAACGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGa                                                                                                                 <  8:222788/135‑1 (MQ=255)
              agATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGa                                                                                                                 <  3:295456/134‑1 (MQ=255)
                 tAAATATGAAAGTGGCATGTACAGTTATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTa                                                                                                            <  2:361080/137‑1 (MQ=255)
                   aaTATGAAAGTGGCATTTACTGTGATTTACTTTTTTTGAAGAGAATAGCTGCATTAGATGAAGATTTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATa                                                                                                         <  2:145200/138‑1 (MQ=255)
                       tataaGTGTCATTTACAGTGATGAACTTTTTTTGAAGAGAAACGCTCCATTAGATGAAGACTATAAAGATCTCCAGAATGCAAAGAATGAGCTCAATGTATTACAAGATACACAATCAGAAATAGACATTAATACAGt                                                                                                     <  4:260178/135‑1 (MQ=255)
                             tGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATa                                                                                                                   <  7:181864/118‑1 (MQ=255)
                             tGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATa                                                                                                                   >  6:237329/1‑118 (MQ=255)
                             tGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATa                                                                                                                   >  8:181864/1‑118 (MQ=255)
                             tGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATa                                                                                                                   <  5:237329/118‑1 (MQ=255)
                               gCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGaaa                                                                                                                     <  7:11741/114‑1 (MQ=255)
                               gCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGaaa                                                                                                                     >  8:11741/1‑114 (MQ=255)
                                cATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAaacaa                                                                                            >  6:70618/1‑138 (MQ=255)
                                    tACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACata                                                                                        >  4:298223/1‑138 (MQ=255)
                                     aCAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGa                                                                                                                 >  4:144786/1‑112 (MQ=255)
                                     aCAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGa                                                                                                                 <  3:144786/112‑1 (MQ=255)
                                       aGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACataaat                                                                                     >  3:5505/1‑138 (MQ=255)
                                       aGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATAAACAATCAGAAATAGACAGTAATACAGTGAGAAACAACataaat                                                                                     >  7:167662/1‑138 (MQ=255)
                                         tgatgaACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACATAATCAGAAATAGACAGTAATACAGTGAGAAACAACataaataa                                                                                   >  5:320262/1‑138 (MQ=255)
                                         tgatgaACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACATAATCAGAAATAGACAGTAATACAGTGAGAAACAACataaataa                                                                                   >  5:314359/1‑138 (MQ=255)
                                         tgatgaACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATa                                                                                                                   >  6:23492/1‑106 (MQ=255)
                                         tgatgaACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATa                                                                                                                   <  5:23492/106‑1 (MQ=255)
                                           atgaACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGa                                                                                                                 <  8:151727/106‑1 (MQ=255)
                                           atgaACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGa                                                                                                                 >  7:151727/1‑106 (MQ=255)
                                           atgaACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATaaataaaa                                                                                 >  7:22645/1‑138 (MQ=255)
                                           atgaACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATaaataaaa                                                                                 >  3:131614/1‑138 (MQ=255)
                                           atgaACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATaaataaaa                                                                                 >  5:24198/1‑138 (MQ=255)
                                                   tttttGCAGAGAACAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         <  4:226624/138‑1 (MQ=255)
                                                   tttttGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACATAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         <  6:314359/138‑1 (MQ=255)
                                                   tttttGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACATAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         <  6:320262/138‑1 (MQ=255)
                                                   tttttGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         >  5:422842/1‑138 (MQ=255)
                                                   tttttGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         <  3:425247/138‑1 (MQ=255)
                                                   tttttGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         <  3:298223/138‑1 (MQ=255)
                                                   tttttGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         <  8:34154/138‑1 (MQ=255)
                                                   tttttGAAGAGAAAAGCTGCATTAGATGAAGAATATAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         <  1:131387/138‑1 (MQ=255)
                                                   tttttGAAGAGAAAAGCTGCATTAGATGAAGAAGTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         <  1:370389/138‑1 (MQ=255)
                                                   tatttGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAAAGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGAc                                                                         <  7:375835/136‑1 (MQ=255)
                                                                        ttAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCt                                                   <  3:299833/139‑1 (MQ=255)
                                                                                agaaTTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACAtt                                                             <  4:295716/121‑1 (MQ=255)
                                                                                agaaTTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACAtt                                                             >  3:295716/1‑121 (MQ=255)
                                                                                agaaTTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGaa                                            >  1:180700/1‑138 (MQ=255)
                                                                                   aTTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGaaa                                           >  3:21053/1‑136 (MQ=255)
                                                                                        aaGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACg                                    <  7:184457/138‑1 (MQ=255)
                                                                                        aaGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACg                                    <  6:24198/138‑1 (MQ=255)
                                                                                        aaGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACg                                    <  3:400635/138‑1 (MQ=255)
                                                                                         agagCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAg                                  <  7:226599/139‑1 (MQ=255)
                                                                                         agagCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAg                                  <  8:112050/139‑1 (MQ=255)
                                                                                         agagCTCCAGAATGCAAAGAATGAGCTCAATGGACTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAg                                  <  2:169272/139‑1 (MQ=255)
                                                                                                    atcaaaaTAATGACCTCAATGGCTTACCAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATATAATAATTGACCAGTAGCACATTGAATCAAGCTCAGAGAAAAAGAACGAGCTATTAAGGAt                       <  2:292049/135‑1 (MQ=255)
                                                                                                    aTGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAGCTATTAAGGAt                       <  5:70618/139‑1 (MQ=255)
                                                                                                                     aaTGGATTACAAGATACACAATCAGAAATAGACTGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAGCTATTAAGGATGGTCTTAAAAGATGTTa      <  4:5505/139‑1 (MQ=255)
                                                                                                                     aaTGGATTACAAGATACACAATAATAAATAGACAGTAATACAGTTAGAAACAACATAAATAAAATTATTGAACAGTATCACATTGAATCAAGATCAGAGAAAAAGAACGAGCTATTAAGGATGGTCTTAAAAGATGTTa      <  8:167662/139‑1 (MQ=255)
                                                                                                                          aTTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAGCTATTAAGGATGGTCTTAAAAGATGTTATTGt  >  4:319162/1‑138 (MQ=255)
                                                                                                                                   tACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAGCTAtt                             >  2:313966/1‑102 (MQ=255)
                                                                                                                                   tACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAGCTAtt                             <  1:313966/102‑1 (MQ=255)
                                                                                                                                       |                                                                                                                            
GAGAACTTCATCTTCGATAAATATGAAAGTGGCATTTACAGTGATGAACTTTTTTTGAAGAGAAAAGCTGCATTAGATGAAGAATTTAAAGAGCTCCAGAATGCAAAGAATGAGCTCAATGGATTACAAGATACACAATCAGAAATAGACAGTAATACAGTGAGAAACAACATAAATAAAATAATTGACCAGTATCACATTGAATCAAGCTCAGAGAAAAAGAACGAGCTATTAAGGATGGTCTTAAAAGATGTTATTGT  >  BS168_reference_genome/2286040‑2286299

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