Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 675,819 C→T 11.4% L13L (TTG→TTA ydjJ ← putative membrane associated potassium channel; prophage region 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome675,8190CT11.4% 93.0 / 3.3 35L13L (TTG→TTAydjJputative membrane associated potassium channel; prophage region 3
Reads supporting (aligned to +/- strand):  ref base C (17/14);  new base T (2/2);  total (19/16)
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.00e+00

ACTTAGAACCTTTGCGTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTATTTCCTTTTATCCTGAGTCTGAAACAATTGAACGAATGTTTGATTTTACAG  >  BS168_reference_genome/675684‑675953
                                                                                                                                       |                                                                                                                                      
aCTTAGAACCTTTGCGTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCaa                                                                                                                                      >  4:54995/1‑138 (MQ=255)
     gAACCTTTGCGTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAACTCAATCATCCCCGCTGTCAAtaataa                                                                                                                                <  7:121018/139‑1 (MQ=255)
      aaCCTTTGCGTGGTCACCGCTGATTTTGCGGCCGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAAtaataa                                                                                                                                >  8:132649/1‑138 (MQ=255)
       tcgtttGCGTGGTCAGCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGAGAGGCTTCTTGTGATCGGGCCGGCACCGCTCTCTAACCTTAACGTGGTTCATAAATAATCAATCATCCCCGCTGTCAATAATaaa                                                                                                                               <  8:134167/135‑1 (MQ=255)
           ttGCGTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGcc                                                                                                                          >  3:194380/1‑139 (MQ=255)
           ttGCGTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGcc                                                                                                                          >  5:161220/1‑139 (MQ=255)
           ttGCGTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACACGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGcc                                                                                                                          >  6:83517/1‑139 (MQ=255)
            tGCGTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGcc                                                                                                                          <  6:85620/138‑1 (MQ=255)
            tGCGTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCt                                                                                                                         >  1:183850/1‑139 (MQ=255)
               gTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCa                                                                                                                      <  4:194380/139‑1 (MQ=255)
                tGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCaa                                                                                                                     <  6:161220/139‑1 (MQ=255)
                tGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACATGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCaa                                                                                                                     <  1:113588/139‑1 (MQ=255)
                   tCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAAtaataa                                                                                                                                <  8:120487/125‑1 (MQ=255)
                   tCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAAtaataa                                                                                                                                >  7:120487/1‑125 (MQ=255)
                            aTTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTTAATAATAAAAGTCCTGCAATGATATAACCCt                                                                                                         >  8:105390/1‑139 (MQ=255)
                            aTTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATACCCGCTGTTAATAATAAAAGGCCTGCAATGATATAACCCt                                                                                                         >  6:188465/1‑139 (MQ=255)
                                 gCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCtcatt                                                                                                     <  2:183850/138‑1 (MQ=255)
                                                   ccAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCt                                                                                   >  7:50932/1‑138 (MQ=255)
                                                      ggCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATACCGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCa                                                                                                                      <  5:179600/100‑1 (MQ=255)
                                                      ggCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATACCGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCa                                                                                                                      >  6:179600/1‑100 (MQ=255)
                                                      ggCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAg                                                                               >  8:150456/1‑139 (MQ=255)
                                                           aCAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGc                                                                          >  1:46711/1‑139 (MQ=255)
                                                           aCAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGc                                                                          <  2:46711/139‑1 (MQ=255)
                                                                     tgtgATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTAt                                                                >  2:78753/1‑139 (MQ=255)
                                                                           cGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTttctt                                                          >  3:31634/1‑139 (MQ=255)
                                                                                  gCACCGCTCTCTAACCATAACGTTGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTAtt                                                   <  4:31634/139‑1 (MQ=255)
                                                                                  gCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTAtt                                                   <  8:50932/139‑1 (MQ=255)
                                                                                        ctctctAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTAATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTATTGTCTATTTACttg                                             >  4:128926/1‑138 (MQ=255)
                                                                                                     aCGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTATTTCCTTTTATCCTGAGTCt                                 >  7:58821/1‑138 (MQ=255)
                                                                                                       gTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCATCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTATATTTCATTTTATCCTGAGTCTGa                               >  6:68734/1‑138 (MQ=255)
                                                                                                        tGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTATTTCCTTTTATCCTGAGGCTGaaa                             <  3:54995/139‑1 (MQ=255)
                                                                                                                                   cTGTTAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTATTTCCTTTTATCCTGAGTCTGAAACAATTGAACGAATGTTTGATTTTACag  <  7:105390/139‑1 (MQ=255)
                                                                                                                                   cTGTTAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTATTTCCTTTTATCCTGAGTCTGAAACAAATGAACGAATGTTTGATTTTACag  <  5:188465/139‑1 (MQ=255)
                                                                                                                                   cTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTATTTCCTTTTATCCTGAGTCTGAAACAATTGAACGAATGTTTGATTTTa     >  7:219650/1‑136 (MQ=255)
                                                                                                                                   cTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTATTTCCTTTTATCCTGAGTCTGAAACAATTGAACGAATGTTTGATTTTACag  <  1:78753/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
ACTTAGAACCTTTGCGTGGTCACCGCTGATTTTGCGGCAGCCCTTCCACAGCCAGGCAGACAGGCATCTTGTGATCGGGCCGGCACCGCTCTCTAACCATAACGTGGTCCATAAAAAATCAATCATCCCCGCTGTCAATAATAAAAGGCCTGCAATGATATAACCCTCATTCATTTTGGTTGCCTCCCTTAAGTCAGCCTTCTTTTATTTTCTTCTCTATTTCCTTTTATCCTGAGTCTGAAACAATTGAACGAATGTTTGATTTTACAG  >  BS168_reference_genome/675684‑675953

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: