Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 966,494 G→C 26.1% A59G (GCG→GGG)  ygaO ← putative integral inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome966,4940GC26.1% 40.3 / 2.7 29A59G (GCG→GGG) ygaOputative integral inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (9/12);  new base C (4/3);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATATATGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCAAGGCCA  >  BS168_reference_genome/966366‑966628
                                                                                                                                 |                                                                                                                                      
cGATATATGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCTGGGCGGCTGTTTTTAGGTTAGGTGAAAAAAACCGCCGGcggcgc                                                                                                                              >  1:67664/1‑139 (MQ=255)
       ttaCGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGATTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCaacaaac                                                                                                                       <  7:24062/137‑1 (MQ=255)
       tgccgcgcgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCaacaaac                                                                                                                       <  6:285233/136‑1 (MQ=255)
       tGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCaacaaac                                                                                                                       <  5:107417/139‑1 (MQ=255)
        gACGCGCTACATGCTTTTGCCTTCATCATTTGCTCTCTTTCATCCTGATCTTCAAAATCATAAGGCAGAAGTAGCTGCTTTCTGGAAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCaacaaac                                                                                                                       <  3:183937/138‑1 (MQ=255)
        gACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTAAAAATAAACCGCAGACGCCGCaacaaac                                                                                                                       <  2:5618/138‑1 (MQ=255)
          cgcgcgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTTAGTACAAATAAACCCCAGACGCCGCaacaaacat                                                                                                                     >  3:132358/1‑137 (MQ=255)
          cgcgcgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTATGGTAGTTAAAATAAAACGCACACGCCCCcacaaacca                                                                                                                     >  4:103417/1‑136 (MQ=255)
          cgcgcgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTAGGTTAATACAAATAAAACCCACACGCCGCcaccaaaac                                                                                                                     >  4:15290/1‑135 (MQ=255)
            cgcgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTGTTTTATTACAAAATAAACGCAAAAGCCCCAACAAACaaaa                                                                                                                   >  7:111547/1‑138 (MQ=255)
             gcgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCGTTTTTTCTGTATGTACAAATAAACCGCAGCCGCCGCCACAAACCGAAt                                                                                                                  >  5:300821/1‑138 (MQ=255)
             gcgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGAAGAAGTAGCTGCTTTCTGGCAGGGCGGCGTTTTTTATGTAAGTACAAATAAACCACAGACGCCGCCACAAAAAAAAt                                                                                                                  >  8:24062/1‑138 (MQ=255)
                       tttGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCTTCAGGCAGAAGTAGCTGCTTTCTGGCTGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCCCCGACGCCGCAACAAACAGAATAAAGGGAGTg                                                                                                        <  3:15290/138‑1 (MQ=255)
                            ggTCATCTTCTGCTCTCTTTAATCCTGCTCTTCAAACTCATCAGGCCGACGTAGCTGCTTTGTGGCAGGGCGGCTTTTTTTCTGAAAGTACAATTAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGttt                                                                                                  <  8:77928/138‑1 (MQ=255)
                            cgacatcatCTGCTCTCTTTCATCCTGATCTTAAAACTCATCTGGAACATGTCGCGGCTTTATGGCCGGGCGGTTTTTTTTCTTTAATAACAAATAAACCCCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGttt                                                                                                  <  2:271666/136‑1 (MQ=255)
                                atcatcTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCGTTTTTTATGTAATTACAAAAAAACCCCAGACCCCCCACCAAACACAAAAAAAGGAGTGAGGtgttaa                                                                                               >  1:55577/1‑133 (MQ=255)
                                     cTGATCTATTTCATCAGGCGCTACAAACTCATCTGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTAAAAATAAACCGCCGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCaaa                                                                                          <  8:111547/138‑1 (MQ=255)
                                        ctctctTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCGTTTTTTGTGTTATGTAAAATAAACCGCCGACGGCGCCACAAACAGAATAAAGGGTGTGGGGTTTAATTAAAATAc                                                                                       >  2:102240/1‑138 (MQ=255)
                                           tatttCATCGTGCTCTTCAAACTCATCAGGCACAAGTAGCGGCTTTCTGGCTGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCACACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCa                                                                                     <  5:155630/135‑1 (MQ=255)
                                           ccatccATCCTCCTCTACAACCACATCCGCCAGAAGTCGCTGCTTCGGGGCTGGGCGGCTTTTTTTCTGTATGACAAAATAAACCGCCCACGCCGCAACAAACAGATTAAAGGGAGGGAGGTTAAAATCAAATACCa                                                                                     <  2:67664/132‑1 (MQ=255)
                                                  tCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTGTTTTTGTGTAATTACAAAAAAACCCCAGGCCCCCCAACAAACAGAATAAAAGGAGTGAGGTTTAAATAAAATACCAACCcgccgc                                                                            >  4:209276/1‑139 (MQ=255)
                                                  tCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCGTTTTTTGTGTTAGTAAAAAAAAACCGCCGACGCCGCAACCAACACAAAAAAGGGCGTGGGGTTTAAATAAAAAACCAACcagcccc                                                                            >  6:78668/1‑137 (MQ=255)
                                                                   atcaGGCAGATGGAGGTGGTTTCTGGCAGGGGGTCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGaaaaa                                                            <  2:55577/138‑1 (MQ=255)
                                                                            ctggtgctgcgTTCTGGGGGGGGGGTTTTTTTTGTTAAATTAAAAAAAACCCCCCCACGCCGCAACAAACAGAATAAAGGGATTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTAt                                                   <  2:119689/133‑1 (MQ=255)
                                                                             aGTAGCTGCTTTCTGGCAGGGCGGGTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATg                                                  <  3:103417/138‑1 (MQ=255)
                                                                                      tttCTGGCAGGGCGGCGTTTTTTCGGGTAGTACAAATAAACCGCCGACGCCGCCACCAACAGAATAAAGGGTGTGAGGTTTAAATCAAATACCCTCCCGCCCCGGTTTGTGGCGGAAAAATAAGCGATGGCGTCGGAgg                                        >  2:69201/1‑139 (MQ=255)
                                                                                               gggCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATaa                                                                                                                <  5:296297/58‑1 (MQ=255)
                                                                                               gggCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATaa                                                                                                                >  6:296297/1‑58 (MQ=255)
                                                                                                                              aCCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTATATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCAAGGCCa  >  3:56163/1‑138 (MQ=255)
                                                                                                                                 |                                                                                                                                      
CGATATATGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCAAGGCCA  >  BS168_reference_genome/966366‑966628

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