Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,128,521 C→T 8.3% Q182Q (CAG→CAA)  yodD ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,128,5210CT 8.3% 153.3 / 3.8 48Q182Q (CAG→CAA) yodDputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base T (2/2);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CCGTCTGCCGCTTTCAAAAGTATCAACTTGGCTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAGGGAGCCCCG  >  BS168_reference_genome/2128385‑2128654
                                                                                                                                        |                                                                                                                                     
ccGTCTGCCGCTTTCAAAAGTATCAACTTGGCTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCt                                                                                                                                      >  7:325623/1‑138 (MQ=255)
ccGTCTGCCGCTTTCAAAAGTATCAACTTGGCTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCt                                                                                                                                      >  6:66406/1‑138 (MQ=255)
  gTCTGCCGCTTTCAAAAGTATCAACTTGGCTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGAt                                                                                                                                   >  4:298231/1‑139 (MQ=255)
     tGCCGCTTTCAAAAGTATCAACTTGGCTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGc                                                                                                                                >  7:303848/1‑139 (MQ=255)
            ttCAAAAGTATCAACTTGGCTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCATCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTc                                                                                                                         >  6:364318/1‑139 (MQ=255)
                         aCTTGGCTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCGCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACg                                                                                                            <  5:88699/139‑1 (MQ=255)
                             ggCTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGa                                                                                                          <  6:355328/137‑1 (MQ=255)
                               cTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCt                                                                                                      <  3:298231/139‑1 (MQ=255)
                                        tatGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGGTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCg                                                                                             <  8:126719/139‑1 (MQ=255)
                                                       tCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAg                                                                              <  6:230478/139‑1 (MQ=255)
                                                                       aGGTTTTTTTATACAATCGCCTCTGTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCtcttct                                                               <  8:196233/138‑1 (MQ=255)
                                                                             ttttATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACg                                                                                                            <  8:208955/87‑1 (MQ=255)
                                                                             ttttATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACg                                                                                                            >  7:208955/1‑87 (MQ=255)
                                                                                 ataCAATCGCCTCTTTGTACCCTTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACTGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTcttcct                                                     <  5:97368/138‑1 (MQ=255)
                                                                                  taCAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATccttcctt                                                    <  2:148861/138‑1 (MQ=255)
                                                                                       tCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGc                                               >  7:47195/1‑138 (MQ=255)
                                                                                       tCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCa                                              <  8:325623/139‑1 (MQ=255)
                                                                                         gCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCaaa                                            <  8:303848/139‑1 (MQ=255)
                                                                                            tctTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAagag                                         >  1:395277/1‑139 (MQ=255)
                                                                                              tttGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGa                                       >  7:407265/1‑139 (MQ=255)
                                                                                                    ccATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTa                                  <  4:296299/138‑1 (MQ=255)
                                                                                                    ccATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTa                                  <  2:395277/138‑1 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGTTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTc                              >  5:55672/1‑139 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGTTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTc                              >  5:50396/1‑139 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTcc                                                                                                                        <  6:402348/49‑1 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTcc                                                                                                                        >  5:402348/1‑49 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAATTAGCGGTAACGCTCTTCTGATTCTTCCTTCGTGCAAAGATGATCGTAtttg                              >  8:103653/1‑138 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAg                                                                      >  6:157621/1‑99 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAg                                                                      <  5:157621/99‑1 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTc                              >  6:321409/1‑139 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTc                              >  8:187684/1‑139 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTc                              >  6:9714/1‑139 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGAGTCTTCCTTCGGGCAAAGAGGATCGTATTTc                              >  6:180260/1‑139 (MQ=255)
                                                                                                       ttCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGAAACGCAAATATCGGTAAAGCTCGTCTGATTCTACCTTCGTGCAATGACGATCGTATTTc                              >  2:176621/1‑139 (MQ=255)
                                                                                                                   ttCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCa                                                                                   <  6:342647/74‑1 (MQ=255)
                                                                                                                   ttCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCa                                                                                   >  5:342647/1‑74 (MQ=255)
                                                                                                                   ttCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATa                  <  7:187684/139‑1 (MQ=255)
                                                                                                                          tCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAg           >  1:126221/1‑139 (MQ=255)
                                                                                                                             tGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGAGGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTGCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAgg          >  2:194905/1‑137 (MQ=255)
                                                                                                                               cTGAGTCAGTTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAGGGAg       <  6:55672/138‑1 (MQ=255)
                                                                                                                               cTGAGTCAGTTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAGGGAg       <  6:50396/138‑1 (MQ=255)
                                                                                                                                  aGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGa                                       <  8:240459/103‑1 (MQ=255)
                                                                                                                                  aGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGa                                       >  7:240459/1‑103 (MQ=255)
                                                                                                                                  aGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAGGGAGcccc   >  6:229803/1‑139 (MQ=255)
                                                                                                                                  aGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAGGGAGcccc   >  6:159432/1‑139 (MQ=255)
                                                                                                                                  aGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGAGAAAAACATGGAGcccc   >  7:356628/1‑139 (MQ=255)
                                                                                                                                   gTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAGGGAGCCCCg  <  5:9714/139‑1 (MQ=255)
                                                                                                                                   gTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAGGGAGCCCCg  <  8:47195/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CCGTCTGCCGCTTTCAAAAGTATCAACTTGGCTGTAAGGATATGAAAAACCTTAGTCCGAAATCCGGACTAAGGTTTTTTTATACAATCGCCTCTTTGTACCATTCACGAGCCTGTTCAGCTTCATGCTGAGTCAGCTGATGGCCCCCGTCCTGCCAGTAAACGGAGGCTGATGCCCCGCTGTCACGCAAATAGCGGTAAAGCTCTTCTGATTCTTCCTTCGTGCAAAGAGGATCGTATTTCCCGGCACCGATAAAAACAGGGAGCCCCG  >  BS168_reference_genome/2128385‑2128654

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: ‑