Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 380,714 G→A 7.3% P1249P (CCG→CCA srfAA → surfactin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome380,7140GA 7.3% 171.7 / 2.6 55P1249P (CCG→CCAsrfAAsurfactin synthetase
Reads supporting (aligned to +/- strand):  ref base G (26/24);  new base A (2/2);  total (28/27)
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.35e-01

TTTATAAAGGCGAACAGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACTTTATATA  >  BS168_reference_genome/380581‑380850
                                                                                                                                     |                                                                                                                                        
tttATAAAGGCGAACAGCTTCCTGCTCTTCGTCTCCAGTATAAGGCCTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAATGAAGTCTTTTCAGGCGAGAGACCGGTGCt                                                                                                                                     >  4:67413/1‑139 (MQ=255)
       aGGCGAACAGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCtt                                                                                                                              <  3:235290/139‑1 (MQ=255)
           gAACAGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGt                                                                                                                          >  1:186618/1‑139 (MQ=255)
               aGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGa                                                                                                                      <  3:237606/139‑1 (MQ=255)
                                      tataAGGACTACGCTGTATTGCAAAGCAGACAGGCTGCTGAAGGGTACAAGACGGACCAGGCTTATTTGAAGGAAGTCTTTGCAGGCGAGCTCCCTGTGCTTCAGCTTCTGTCCTATTACCCAAGACCACCTGTTCat                                                                                                <  5:224922/138‑2 (MQ=255)
                                           ggACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAACGCttt                                                                                          >  2:126098/1‑139 (MQ=255)
                                             agtACGCTGTATGGCAAAGCAGTCAGGCTGCTGTAGGGTACATGTAGGAACTTGCTTATTGGAAGGCTGTCTTTGCAGGCGAGCTGCCGGTTCGTCAGCGGCTGTCCGATTACCGAATACCAGCTGTTCAAAGCTTTGa                                                                                        <  8:78890/137‑1 (MQ=255)
                                             aCTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCGTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCCGCGTCTGGCCGAGAACCCAAGCCCCCCTGGTCAAAGCTTTg                                                                                         >  4:96864/1‑138 (MQ=255)
                                                aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGGAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTc                                                                                                                                   <  4:163542/93‑1 (MQ=255)
                                                aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGGAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTc                                                                                                                                   >  3:163542/1‑93 (MQ=255)
                                                aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAgg                                                                                     >  6:224922/1‑139 (MQ=255)
                                                aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCAGTGCtt                                                                                                                                    >  3:238280/1‑92 (MQ=255)
                                                aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCAGTGCtt                                                                                                                                    >  3:234142/1‑92 (MQ=255)
                                                aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCAGTGCtt                                                                                                                                    <  4:234142/92‑1 (MQ=255)
                                                aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCAGTGCtt                                                                                                                                    <  4:238280/92‑1 (MQ=255)
                                                      tATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGCCCACCTGTTCAAAGCTTTGAAGGGGATCg                                                                               >  2:94252/1‑139 (MQ=255)
                                                      tATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGCAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTAAGCTTCTGTCACATTCAGCAACACCACCTGTTCAAAGATTTGAAGGGGATcc                                                                               >  7:78890/1‑138 (MQ=255)
                                                           cAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGgtgt                                                                          <  3:34964/139‑1 (MQ=255)
                                                              aGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGtcc                                                                        >  7:159049/1‑137 (MQ=255)
                                                              aGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGtcaa                                                                       >  2:221525/1‑139 (MQ=255)
                                                                cAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGtcaat                                                                      >  2:28257/1‑138 (MQ=255)
                                                                        ctgctgAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGa                                                             >  3:5472/1‑139 (MQ=255)
                                                                          gctgAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATg                                                           >  1:96116/1‑139 (MQ=255)
                                                                           ctgAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGc                                                          <  1:132907/139‑1 (MQ=255)
                                                                              aaGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCg                                                         <  1:221525/137‑1 (MQ=255)
                                                                                 ggTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCggggg                                                     >  3:42891/1‑138 (MQ=255)
                                                                                 ggTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACAACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGTATTCggggg                                                     >  5:225939/1‑138 (MQ=255)
                                                                                 ggTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACAACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCggggg                                                     >  7:282869/1‑138 (MQ=255)
                                                                                   tACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCg                                                         <  4:10071/132‑1 (MQ=255)
                                                                                   tACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCg                                                         >  3:10071/1‑132 (MQ=255)
                                                                                        gaagGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGAt                                             <  2:96116/139‑1 (MQ=255)
                                                                                           ggACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAGAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGc                                          >  2:15448/1‑139 (MQ=255)
                                                                                             aCCAGGCTTATTGGAATGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGttaa                                                                                                 <  6:237502/82‑3 (MQ=255)
                                                                                             aCCAGGCTTATTGGAATGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCa                                                                                                 >  5:237502/1‑82 (MQ=255)
                                                                                             aCCAGGCTTATTGGAAGGAAATCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCt                                        >  7:79254/1‑139 (MQ=255)
                                                                                             aCCAGGCTTATAGGAAGGAAGTCTTTGCAGCCGACCTCCCGGTGCTTCAGCTACAGTCCCAATACCCAAGACCACCTGTTAAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATCCGGGAGTAAAGTATCGcat                                        >  8:35459/1‑137 (MQ=255)
                                                                                                aGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCaa                                     >  8:259238/1‑139 (MQ=255)
                                                                                                            aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGaccac                                                                                                        <  2:134328/60‑1 (MQ=255)
                                                                                                            aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGaccac                                                                                                        >  1:134328/1‑60 (MQ=255)
                                                                                                            aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGa                                                                                        >  3:205195/1‑76 (MQ=255)
                                                                                                            aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGa                                                                                        <  4:205195/76‑1 (MQ=255)
                                                                                                            aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTgg                             >  5:243219/1‑135 (MQ=255)
                                                                                                            aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTgg                             <  6:243219/135‑1 (MQ=255)
                                                                                                                       gCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGgcg              <  1:126098/139‑1 (MQ=255)
                                                                                                                        ccggCGAGCGCCCGGGGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGgcgc             <  1:245641/137‑1 (MQ=255)
                                                                                                                        cAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTGCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGgcg              <  6:121027/138‑1 (MQ=255)
                                                                                                                         aGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTc                                       <  7:236606/112‑1 (MQ=255)
                                                                                                                         aGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTc                                       >  8:236606/1‑112 (MQ=255)
                                                                                                                            cgcgCTCCCGGGGCTTCAGCTGCTGTCCGATAACCCAAGACCACCTGTTCAACGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGGGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACt         <  6:206360/136‑1 (MQ=255)
                                                                                                                            cGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGGGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACt         <  3:80143/139‑1 (MQ=255)
                                                                                                                            cGAGCTCCCGGTGCTTCAGCTGCGGTCCGATTACCCAAGACCACCTGGTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGAGCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACt         <  6:225939/139‑1 (MQ=255)
                                                                                                                                cccccGGGGCGTCAGCGTCTGCCCGATACCCCACCACCACCTGTTCAAAGCTTTGACGGGGATCGGGTGTCAATCACGCGGGAGGCGGGGGTAACGGCTCGCCGCAATCGTTTGGCTCAAAAAAACGGCGCCACTTtat     <  4:42891/137‑1 (MQ=255)
                                                                                                                                   ccGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAACGCTTTGAAGGGGATCGGGTGTCAATCAAGCGGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACTTtata    <  4:84581/137‑1 (MQ=255)
                                                                                                                                   ccGGTGCTGCAGCTTCTGACCGATTACCCAACACCACCTGTTCAAAGCTTTGAAGGGGAGCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACTTtata    <  4:5472/137‑1 (MQ=255)
                                                                                                                                   ccGGTGCTCCAGCTTCTGTCCGATTACCCACGACCACCTGTGCAAAGCTTTGAAGGGCATCGGGGGCCAAAAAAGCTGGATGGGGGGGTAACGGAGGGCCACAATCGTTTGCCCCAACAAAACGGCGCCACTTtatata  <  8:159049/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                        
TTTATAAAGGCGAACAGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACTTTATATA  >  BS168_reference_genome/380581‑380850

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: