Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,369,488 T→G 6.4% G293G (GGA→GGC tyrA ← prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,369,4880TG 6.4% 174.3 / 3.3 58G293G (GGA→GGCtyrAprephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base G (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GTTAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCT  >  BS168_reference_genome/2369350‑2369623
                                                                                                                                          |                                                                                                                                       
gTTAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCt                                                                                                                                         <  8:254790/139‑1 (MQ=255)
 ttAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTc                                                                                                                                        >  1:159568/1‑139 (MQ=255)
 ttAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTc                                                                                                                                        >  1:6330/1‑139 (MQ=255)
           ctcGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGcc                                                                                                                               >  8:342345/1‑138 (MQ=255)
                 tttCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCgg                                                                                                                        >  6:29569/1‑139 (MQ=255)
                 tttCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGATCCCTTCTGCCGAAGCgg                                                                                                                        >  7:26859/1‑139 (MQ=255)
                      aTAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGcc                                                                                                                   >  2:40467/1‑139 (MQ=255)
                         aTGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGt                                                                                                                <  2:159568/139‑1 (MQ=255)
                              gATATTCGTGATACTGATGCGCTCCGCTGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGa                                                                                                           <  6:262072/139‑1 (MQ=255)
                              gATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGa                                                                                                           <  1:301912/139‑1 (MQ=255)
                              gATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGa                                                                                                           <  1:40467/139‑1 (MQ=255)
                                   tCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATc                                                                                                      >  8:137810/1‑139 (MQ=255)
                                      tGATACTGATGCGCTCCGCCGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGAGCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGCTaaaactt                                                                                                   >  8:11994/1‑134 (MQ=255)
                                         tACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGc                                                                                                <  2:6330/139‑1 (MQ=255)
                                             gATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGtt                                                                                            >  6:296701/1‑139 (MQ=255)
                                              aTGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTgg                                                                                                 <  1:39911/133‑1 (MQ=255)
                                              aTGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTgg                                                                                                 >  2:39911/1‑133 (MQ=255)
                                              aTGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGttt                                                                                           >  3:164553/1‑139 (MQ=255)
                                                   cTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCGGTTTTtaaa                                                                                      >  4:48299/1‑139 (MQ=255)
                                                                atagccgTTATCTCGGAAATCACACCGGGATGATCGGGAACATCCACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATAATCCTTGGCCGTTTTAAAATAACGAAATAGAt                                                                         >  8:338355/7‑139 (MQ=255)
                                                                aTCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTaaaa                                                                                      >  2:1863/1‑126 (MQ=255)
                                                                aTCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTaaaa                                                                                      <  1:1863/126‑1 (MQ=255)
                                                                 tCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGAtt                                                                        <  6:129740/139‑1 (MQ=255)
                                                                  cGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGAttt                                                                       >  5:296230/1‑139 (MQ=255)
                                                                       ttATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCg                                                                   <  5:29569/138‑1 (MQ=255)
                                                                             cGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGGATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc                                                            <  8:151748/139‑1 (MQ=255)
                                                                              ggATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc                                                            <  2:28048/138‑1 (MQ=255)
                                                                              ggATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc                                                            <  7:360178/138‑1 (MQ=255)
                                                                              ggATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtc                                                               <  8:67941/135‑1 (MQ=255)
                                                                               gATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCATGGCGGTTTTAAAAAAAAGAAATAGATTATACGCAtctaatt                                                          >  8:360178/1‑135 (MQ=255)
                                                                               gATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCAGTGGCTGTTTTAAAAAAACGAAATATATTTTTCGCATcttctt                                                          >  1:170815/1‑139 (MQ=255)
                                                                                atatTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTg                                                         <  4:226755/139‑1 (MQ=255)
                                                                                   ttACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGtt                                                       <  7:373822/138‑1 (MQ=255)
                                                                                   ttACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTc                                                      <  8:30378/139‑1 (MQ=255)
                                                                                        cccGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACat                                                 <  5:296701/139‑1 (MQ=255)
                                                                                        cccGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACat                                                 <  1:364590/139‑1 (MQ=255)
                                                                                        cccGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACat                                                 <  7:81909/139‑1 (MQ=255)
                                                                                           ggATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCGTGGCTGTTTTTAAATAACGAAATAGATTTTCCGCATATTATTGTTCTACatatt                                              >  6:329295/1‑138 (MQ=255)
                                                                                              tgatCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGt                                                                                                                >  6:127024/1‑70 (MQ=255)
                                                                                              tgatCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGt                                                                                                                <  5:127024/70‑1 (MQ=255)
                                                                                              tgatCGGGAACATCCACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATa                                                                                                         >  1:271903/1‑77 (MQ=255)
                                                                                              tgatCGGGAACATCCACATATAAATCGTAAAATGCAGGTATCGCGCCCTTCTGCCGAAGCGGCAGTCCGTCGCGATa                                                                                                         <  2:271903/77‑1 (MQ=255)
                                                                                                     gTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTgg                                                                                                 <  7:220841/78‑1 (MQ=255)
                                                                                                     gTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTgg                                                                                                 >  8:220841/1‑78 (MQ=255)
                                                                                                              caTGTAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAAttt                            >  3:397396/1‑138 (MQ=255)
                                                                                                               atatAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCa                          <  5:249861/139‑1 (MQ=255)
                                                                                                               atatAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAAACCTTGGGTGTTTTAAAATAACGAAATAGATTTTACGCATCGTCTTGTACTACAAAATTACGGAGATTGTGAAATTCa                          >  4:307587/1‑139 (MQ=255)
                                                                                                                     aTCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGaaa                                                                              >  7:116529/1‑81 (MQ=255)
                                                                                                                     aTCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGaaa                                                                              <  8:116529/81‑1 (MQ=255)
                                                                                                                     aTCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATc                    >  1:24845/1‑139 (MQ=255)
                                                                                                                         tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGCGTATCCGTTGGGGGTTTAAAAAAACGAAATAGATTTTACGCATCGTGTTGTTATACATATGTTCGGTTCTTGTCAATTTCACGACTCCCCt                >  8:276584/1‑139 (MQ=255)
                                                                                                                         tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc                                                            <  4:372027/95‑1 (MQ=255)
                                                                                                                         tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCAtcttc                                                            >  3:372027/1‑95 (MQ=255)
                                                                                                                         tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCa                  >  2:272227/1‑137 (MQ=255)
                                                                                                                         tAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCa                  <  1:272227/137‑1 (MQ=255)
                                                                                                                                cAGGTAGCGCGCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaa         <  6:49289/139‑1 (MQ=255)
                                                                                                                                 aGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCaaaa        <  7:342345/139‑1 (MQ=255)
                                                                                                                                  ggTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAAc       <  1:406003/139‑1 (MQ=255)
                                                                                                                                       cGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCt  <  1:379904/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                       
GTTAATATCCTCTCGTGTTTCGATAATGCGGATATTCGTGATACTGATGCGCTCCGCAGCTAAGATCGCTGTTATCTCGGATATTACACCCGGATGATCGGGTACATCCACATATAAATCATAAAATGCAGGTATCGCTCCCTTCTGCCGAAGCGGCAGCCCGTCGCGATAATCCTTGGCTGTTTTAAAATAACGAAATAGATTTTCCGCATCTTCTTGTTCTACATATGTACGGATCTTGTCAATTTCACGAATCCACTCATCAAAACGGTCT  >  BS168_reference_genome/2369350‑2369623

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: