Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 1,317,152 +GT intergenic (+23/+5) phrA → / ← yjpA secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,317,1521.G94.8% 208.0 / 4.2 59intergenic (+23/+5)phrA/yjpAsecreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (35/20);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*BS168_reference_genome1,317,1522.T93.2% 183.0 / 5.2 59intergenic (+23/+5)phrA/yjpAsecreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (35/20);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CTCGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑C‑‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATAC  >  BS168_reference_genome/1317025‑1317271
                                                                                                                                 ||                                                                                                                          
cTCGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑G‑‑‑tcttttt                                                                                                                 <  3:162857/136‑3 (MQ=255)
        ggTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑G‑‑‑TCttttttatttat                                                                                                          >  7:335171/1‑129 (MQ=255)
                           cAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATc                                                                                         >  1:372854/1‑136 (MQ=255)
                            aGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATATCATATTGCGGCACGCAATCAAACATGATGCATAAAAAACGACCCTTAGCGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTAt                                                                                          <  3:184104/134‑1 (MQ=255)
                            aGGTGATGGGTCATGCAGGTGAAACAGCAAACACAGAATGGAATACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACACTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTAt                                                                                          <  4:313620/134‑1 (MQ=255)
                                 aTGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAg                                                                                    <  8:13265/135‑1 (MQ=255)
                                 aTGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAg                                                                                    <  5:26779/135‑1 (MQ=255)
                                 aTGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGc                                                                                   >  3:154331/1‑136 (MQ=255)
                                 aTGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGc                                                                                   >  7:312804/1‑136 (MQ=255)
                                    gTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTc                                                                                 <  4:154331/135‑1 (MQ=255)
                                    gTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCa                                                                                >  8:383597/1‑136 (MQ=255)
                                      tCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCag                                                                               >  5:249535/1‑135 (MQ=255)
                                           cAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc                                                                         >  4:219856/1‑136 (MQ=255)
                                           cAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc                                                                         >  4:202317/1‑136 (MQ=255)
                                            aGGTGAAACAGCAAACACAGAAGGGAAAACATATCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATcc                                                                        <  4:298810/136‑1 (MQ=255)
                                                 aaaCAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAg                                                                                    <  1:69033/119‑1 (MQ=255)
                                                 aaaCAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAg                                                                                    >  2:69033/1‑119 (MQ=255)
                                                 aaaCAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTg                                                                   >  6:154636/1‑136 (MQ=255)
                                                     aGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGccac                                                               <  7:65532/136‑1 (MQ=255)
                                                      gCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGt                                                                                       >  3:246919/1‑111 (MQ=255)
                                                      gCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGt                                                                                       <  4:246919/111‑1 (MQ=255)
                                                      gCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGccac                                                               >  4:332491/1‑135 (MQ=255)
                                                      gCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGccac                                                               >  8:338127/1‑135 (MQ=255)
                                                           cacaGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCgg                                                         >  4:245181/1‑136 (MQ=255)
                                                           cacaGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCgg                                                         >  2:357922/1‑136 (MQ=255)
                                                           cacaGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCg                                                          >  8:316363/1‑135 (MQ=255)
                                                           cacaGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCg                                                          >  1:114156/1‑135 (MQ=255)
                                                             caGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCATTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCTCCAGCGGAt                                                       <  4:150154/136‑1 (MQ=255)
                                                                  gggAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCa                                                  >  5:324981/1‑136 (MQ=255)
                                                                           tttCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGAtt                                         >  3:263970/1‑136 (MQ=255)
                                                                             tCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTc                                       >  6:6094/1‑136 (MQ=255)
                                                                             tCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTc                                       >  7:103617/1‑136 (MQ=255)
                                                                               atatTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTcc                                      <  7:383597/135‑1 (MQ=255)
                                                                                tatTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTTTTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGAGCCGCATGGCCGAgttaggc                                    >  7:359404/1‑132 (MQ=255)
                                                                                tatTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGc                                    >  7:14092/1‑136 (MQ=255)
                                                                                     cGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGcc                                <  5:380441/135‑1 (MQ=255)
                                                                                     cGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGcc                                <  5:6094/135‑1 (MQ=255)
                                                                                      ggCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCtt                                                                                                      <  7:412023/64‑1 (MQ=255)
                                                                                      ggCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCtt                                                                                                      >  8:412023/1‑64 (MQ=255)
                                                                                         aCGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTc                           <  3:245181/136‑1 (MQ=255)
                                                                                             aaTCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCCCTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGGCTCTTCCAGt                       >  5:115268/1‑136 (MQ=255)
                                                                                                      tgatgCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCt              <  6:280615/136‑1 (MQ=255)
                                                                                                      tgatgCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCt              <  4:263970/136‑1 (MQ=255)
                                                                                                      tgatgCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCt              <  7:338127/136‑1 (MQ=255)
                                                                                                        atgCATAAAAAAAGACCCGTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCt              <  2:325953/134‑1 (MQ=255)
                                                                                                          gCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGCTCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCGTCGa          >  4:45372/1‑136 (MQ=255)
                                                                                                          gCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCCTGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTc            >  4:45373/1‑134 (MQ=255)
                                                                                                          gCATAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGa          >  7:169682/1‑136 (MQ=255)
                                                                                                            ctacaaacaGACCCTTAGGG‑‑CTCTTTTTTATTTCTTCAGCGTCCATTCTTGTATCGACCGCACAGAAGATCCACTTGCCACCAGCGGAGCCGCAGGGCCGATTTCCGCGGCCGCTTCCAGTGAATCTGCTTCGATg        <  7:331404/132‑1 (MQ=255)
                                                                                                            aTAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGCAGATCCCCTTGCCCCCAGCGGATCCGCATGGCCGATTTCCGCGGCCTCTTCCAGTGAATCGGCGTCGCTg        >  6:53386/1‑136 (MQ=255)
                                                                                                            aTAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGGTTTCCGCGGCCGGGTGCAGGGTATCTGCGTCGATg        >  3:232980/1‑136 (MQ=255)
                                                                                                            aTAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGAt         >  8:71893/1‑135 (MQ=255)
                                                                                                            aTAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGAt         >  5:278/1‑135 (MQ=255)
                                                                                                            aTAAAAAAAGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGAt         >  1:82127/1‑135 (MQ=255)
                                                                                                                aaaaaGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGcc                                >  4:17297/1‑108 (MQ=255)
                                                                                                                aaaaaGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGcc                                <  3:17297/108‑1 (MQ=255)
                                                                                                                aaaaaGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACa     >  5:85814/1‑135 (MQ=255)
                                                                                                                aaaaaGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACa     >  5:134946/1‑135 (MQ=255)
                                                                                                                  aaaGACCCTTAGGGGTCT‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATAc  <  6:324981/136‑1 (MQ=255)
                                                                                                                                 ||                                                                                                                          
CTCGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑C‑‑‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATAC  >  BS168_reference_genome/1317025‑1317271

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