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.G96.2% 194.3 / 0.5 54intergenic (+23/+5)phrA/yjpAsecreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (25/26);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01
*BS168_reference_genome1,317,1522.T94.3% 150.5 / 3.2 54intergenic (+23/+5)phrA/yjpAsecreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base T (24/26);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑CTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATACGC  >  BS168_reference_genome/1317023‑1317273
                                                                                                                                   ||                                                                                                                          
tgCTCGTTGCGGTCGGGTTGAGCTTTTCTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTTCGGCACGCAATCAAACATGATGCAAAAAAAAAGACCCTTAGGG‑‑GTCtttttt                                                                                                                  >  1:158574/1‑139 (MQ=255)
    cGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑GTCttttttattt                                                                                                              >  8:84656/1‑135 (MQ=255)
                gTTCAGCTTTACTCAGGTGATGGGTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCa                                                                                                     <  3:126098/139‑1 (MQ=255)
                          aCTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGTGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATc                                                                                           >  4:225879/1‑139 (MQ=255)
                           cTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCTGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGATTCCATTCTTTTATCg                                                                                          >  2:308744/1‑139 (MQ=255)
                             cAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGt                                                                                         >  5:90906/1‑138 (MQ=255)
                              aGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGt                                                                                         <  7:219382/137‑1 (MQ=255)
                              aGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGt                                                                                         <  7:83679/137‑1 (MQ=255)
                                gTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGc                                                                                     >  5:259779/1‑139 (MQ=255)
                                gTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGc                                                                                     >  6:39044/1‑139 (MQ=255)
                                  gATGGATCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTc                                                                                   <  2:188117/139‑1 (MQ=255)
                                     ggTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCag                                                                                 >  4:30528/1‑138 (MQ=255)
                                     ggTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTCTTCAGCTTCCATTCTTTTATCGTCAGCTCag                                                                                 >  3:318641/1‑137 (MQ=255)
                                         cATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCagaaga                                                                             <  4:158642/138‑1 (MQ=255)
                                          aTGCAGGTGAAACAGCAAACACAGCAGGGAAAACATTACATATATCGGCACGCAATCAAACATGATCCATAAAAAAAGACCCTAAGGGGACTTTTTTATTACTTCAGCTTCCATTCTTATATCTTCATATCTGAAGATc                                                                           >  8:158418/1‑139 (MQ=255)
                                          aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGt                                                                                         >  8:275951/1‑125 (MQ=255)
                                          aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGt                                                                                         <  7:275951/125‑1 (MQ=255)
                                          aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc                                                                           >  6:235073/1‑139 (MQ=255)
                                          aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc                                                                           >  7:72008/1‑139 (MQ=255)
                                          aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc                                                                           >  3:331771/1‑139 (MQ=255)
                                          aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc                                                                           >  7:59853/1‑139 (MQ=255)
                                                  taaaCAGCAAACACAGCATGGAAAACAATTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGcc                                                                   <  8:203909/138‑1 (MQ=255)
                                                    aaCAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGccac                                                                 <  7:250252/139‑1 (MQ=255)
                                                          aaaCACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTTTTTATTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATACACTTGcccccccccg                                                           >  2:173251/1‑134 (MQ=255)
                                                           aaCACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGa                                                          <  1:173251/139‑1 (MQ=255)
                                                            acacAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATCATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGAt                                                         <  3:263583/139‑1 (MQ=255)
                                                                 gAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCa                                                    >  2:97259/1‑139 (MQ=255)
                                                                 gAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCa                                                    <  1:97259/139‑1 (MQ=255)
                                                                              ttCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATgg                                                 <  7:346422/129‑1 (MQ=255)
                                                                              ttCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATgg                                                 >  8:346422/1‑129 (MQ=255)
                                                                              ttCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCg                                       >  3:181205/1‑139 (MQ=255)
                                                                                  tatTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGc                                   <  3:70112/139‑1 (MQ=255)
                                                                                     tGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCt                                     <  4:42025/134‑1 (MQ=255)
                                                                                     tGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCt                                     >  3:42025/1‑134 (MQ=255)
                                                                                      gCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCctct                               >  1:264863/1‑139 (MQ=255)
                                                                                                aTCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAAt                     <  1:80437/139‑1 (MQ=255)
                                                                                                     aCATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCt                <  3:116023/139‑1 (MQ=255)
                                                                                                      cATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCtt               <  6:63971/139‑1 (MQ=255)
                                                                                                        tgattCATAAAAAAAGACCCCTAGGGGTCTTTTTTATTTCTTCAGCCTACATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCg             <  1:308744/139‑1 (MQ=255)
                                                                                                        tgatgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCg             <  3:59375/139‑1 (MQ=255)
                                                                                                        tgatgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCg             <  8:59853/139‑1 (MQ=255)
                                                                                                        tgatgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCg             <  8:72008/139‑1 (MQ=255)
                                                                                                         gatgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGa            >  1:232735/1‑139 (MQ=255)
                                                                                                         gatgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGa            >  5:293555/1‑139 (MQ=255)
                                                                                                         gatgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGa            >  5:236982/1‑139 (MQ=255)
                                                                                                          atgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGAt           <  2:158574/139‑1 (MQ=255)
                                                                                                           tgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGAGCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGAt           >  1:79842/1‑138 (MQ=255)
                                                                                                             cATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCg             >  7:177324/1‑134 (MQ=255)
                                                                                                             cATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCg             <  8:177324/134‑1 (MQ=255)
                                                                                                             cATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGa         >  6:126705/1‑138 (MQ=255)
                                                                                                              aTAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACa       <  2:232735/139‑1 (MQ=255)
                                                                                                               tAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATACACTTGCCACCCGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACa       >  6:332572/1‑138 (MQ=255)
                                                                                                                 aaaaaaGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATa     <  6:236982/138‑1 (MQ=255)
                                                                                                                   aaaaGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATACGc  <  1:128725/139‑1 (MQ=255)
                                                                                                                                   ||                                                                                                                          
TGCTCGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑CTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATACGC  >  BS168_reference_genome/1317023‑1317273

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: