Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,315,315 T→A 7.6% H272L (CAT→CTT)  hom ← homoserine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,315,3150TA 7.6% 149.9 / 3.9 53H272L (CAT→CTT) homhomoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base A (2/2);  total (26/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 = 6.89e-01

ACAGCGACGAGGTCAGAAACAACGGATGTCGCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGATGTCCTCGTCTGTAATTTGGGAGATCCCCT  >  BS168_reference_genome/3315180‑3315453
                                                                                                                                       |                                                                                                                                          
aCAGCGACGAGGTCAGAAACAACGGATGTCGCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGgatgat                                                                                                                                         >  6:42380/1‑139 (MQ=255)
       cgaGGTCAGAAACAACGGATGTCGCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAgg                                                                                                                                  <  4:359830/139‑1 (MQ=255)
              aGAAACAACGGATGTCGCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAgt                                                                                                                           >  3:389569/1‑139 (MQ=255)
              aGAAACAACGGATGGCGCTGTCGGCATGCTTCCGGCTCCCGGCCCGGAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAACGCGGATGATGGTCAGGAAGCAgt                                                                                                                           <  8:175784/139‑1 (MQ=255)
                 aacaacGGATGTCGCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACCTCGTCTCACCGACAGCCTCGACGTATTCATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAgtgt                                                                                                                         <  2:212686/138‑1 (MQ=255)
                      cGGATGTCGCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTg                                                                                                                   >  6:141542/1‑139 (MQ=255)
                              gCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGa                                                                                                           >  6:188153/1‑139 (MQ=255)
                              gCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCGGa                                                                                                           >  2:179834/1‑139 (MQ=255)
                              gCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCCTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGa                                                                                                           >  8:225486/1‑139 (MQ=255)
                                   cGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCGGTACGCTGACCTCg                                                                                                      >  1:123014/1‑139 (MQ=255)
                                                    ccGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACg                                                                                     <  6:87790/139‑1 (MQ=255)
                                                    ccGGACCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGTATGATGGTCAGGAAGCTGTGTAGGCTGTACTCTGACCTCTATTTTGCTGCCGTCACg                                                                                     <  6:349785/139‑1 (MQ=255)
                                                     cGGCCCGTAGCACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGGCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGc                                                                                    >  3:400014/1‑139 (MQ=255)
                                                            tAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCaa                                                                             >  6:207731/1‑139 (MQ=255)
                                                                 cATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCTTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCg                                                                        <  6:270600/139‑1 (MQ=255)
                                                                   tCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGAt                                                                      <  7:225486/139‑1 (MQ=255)
                                                                    cGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCg                                                                                          <  2:129917/118‑1 (MQ=255)
                                                                    cGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCg                                                                                          >  1:129917/1‑118 (MQ=255)
                                                                      tctcACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAg                                                                   <  5:42380/139‑1 (MQ=255)
                                                                         cACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTCTGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCtt                                                                <  4:97638/139‑1 (MQ=255)
                                                                         cACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCtt                                                                >  6:194065/1‑139 (MQ=255)
                                                                           ccGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCa                                                              <  5:141542/139‑1 (MQ=255)
                                                                           ccGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCa                                                              <  3:355000/139‑1 (MQ=255)
                                                                           ccGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCa                                                              >  7:30654/1‑139 (MQ=255)
                                                                           ccGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGGCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCa                                                              <  3:346876/139‑1 (MQ=255)
                                                                                aGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGta                                                         <  5:188153/139‑1 (MQ=255)
                                                                                aGCCTCGCCGTATACATAAACAGCGGTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGta                                                         <  6:104613/139‑1 (MQ=255)
                                                                                            taCTTAAACTCCGTTAAACTCATTATGAACAGCTTAAAGGGCATGATGGTCAGGAAGCTGTGGCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTTTGCAATCCCGATCAGCTTCTTTTTTTAGCCGAgggt                                              <  4:389569/138‑3 (MQ=255)
                                                                                            taCATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACg                                                                                                               >  7:257904/1‑73 (MQ=255)
                                                                                            taCATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACg                                                                                                               <  8:257904/73‑1 (MQ=255)
                                                                                            taCATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGCCCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGt                                             >  2:28001/1‑139 (MQ=255)
                                                                                                 aaaCAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGACGGCCGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCt                                        >  6:76325/1‑139 (MQ=255)
                                                                                                      gCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGa                                                                                                     >  6:105376/1‑73 (MQ=255)
                                                                                                      gCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGa                                                                                                     <  5:105376/73‑1 (MQ=255)
                                                                                                      gCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGCGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGc                                   >  7:286989/1‑139 (MQ=255)
                                                                                                      gCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGc                                   <  5:181498/139‑1 (MQ=255)
                                                                                                       cGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCa                                                              >  8:210823/1‑111 (MQ=255)
                                                                                                       cGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCa                                                              <  7:210823/111‑1 (MQ=255)
                                                                                                       cGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCt                                  >  7:317632/1‑139 (MQ=255)
                                                                                                       cGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCt                                  >  4:223536/1‑139 (MQ=255)
                                                                                                        gTTAAACTCATTATGAACAGCAGAAAGCGGCTGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTg                                 <  4:400014/139‑1 (MQ=255)
                                                                                                         ttAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGCCCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGa                                >  7:190768/1‑139 (MQ=255)
                                                                                                         ttAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCTGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGa                                >  3:126042/1‑139 (MQ=255)
                                                                                                                    attaaCAGCAGAAAGCGGATGATTGTCAGGAAGCAGTGTCGGCTGTACGCTGAACTCGATTTTGCAGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGATGTCCTCGTc                      <  4:217055/135‑1 (MQ=255)
                                                                                                                    aTGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGATGTCCTCGTc                      <  8:286989/138‑1 (MQ=255)
                                                                                                                         cagcagAAAGCGGAAGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGta                                                         <  4:362680/98‑1 (MQ=255)
                                                                                                                         cagcagAAAGCGGAAGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGta                                                         >  3:356741/1‑98 (MQ=255)
                                                                                                                         cagcagAAAGCGGAAGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGta                                                         >  3:362680/1‑98 (MQ=255)
                                                                                                                         cagcagAAAGCGGAAGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGta                                                         <  4:356741/98‑1 (MQ=255)
                                                                                                                           gcagAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTGATTGTATAGCCGCGGCGTGTTCTGGAGCTGATGTCCTCGAAGGTAAttt              >  6:228580/1‑139 (MQ=255)
                                                                                                                           gcagAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGATGTCCTCGTCTGTAAttt              >  7:69862/1‑139 (MQ=255)
                                                                                                                           gcagAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTGTGCTGAAGCTGCTGTCCTCGTCTGTAAttt              >  2:282799/1‑139 (MQ=255)
                                                                                                                                aaGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGATGTCCTCGTCTGTAATTTGGgag         <  1:259688/139‑1 (MQ=255)
                                                                                                                                  gCGGATGATGGTCAGGAAGCAGTGACGGCCGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGATGTCCTCGTCTGTAATTTGGGAGAt       <  5:76325/139‑1 (MQ=255)
                                                                                                                                      atgatgGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTtgctg                                                                                             >  7:359777/1‑49 (MQ=255)
                                                                                                                                      atgatgGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTtgctg                                                                                             <  8:359777/49‑1 (MQ=255)
                                                                                                                                       tgatgGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGATGTCCTCGTCTGTAATTTGGGAGATCCCCt  <  3:48684/139‑1 (MQ=255)
                                                                                                                                       tgatgGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGATGTCCTCGTCTGTAATTTGGGAGATCCCCt  <  5:207731/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
ACAGCGACGAGGTCAGAAACAACGGATGTCGCTGTCGGCATGCTTCCGGCTCCCGGCCCGTAGAACATCGTCTCACCGACAGCCTCGCCGTATACATAAACAGCGTTAAACTCATTATGAACAGCAGAAAGCGGATGATGGTCAGGAAGCAGTGTCGGCTGTACGCTGACCTCGATTTTGCTGCCGTCACGCTGAGCAATCCCGATCAGCTTCATTGTATAGCCGAGGCGTTTGCTGAAGCTGATGTCCTCGTCTGTAATTTGGGAGATCCCCT  >  BS168_reference_genome/3315180‑3315453

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