Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,820,393 G→A 6.6% A28V (GCA→GTA)  dtd ← gly‑tRNA(Ala) deacylase / D‑Tyr‑tRNATyr deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,820,3930GA 6.6% 197.4 / 2.9 60A28V (GCA→GTA) dtdgly‑tRNA(Ala) deacylase / D‑Tyr‑tRNATyr deacylase
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base A (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

TTTTTCGTATCTCCATAAAGCGTAAACTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGAACAACTAATCTCTTTTCTAACCCCTT  >  BS168_reference_genome/2820258‑2820529
                                                                                                                                       |                                                                                                                                        
ttttttGTTTCTCCATAAAGCGTAAACTGTGAGCCGGATTAAAACTGTCCTCCTTTTTGTACAAGAGAAAGAATCAATTTACCTTCACTGTCATCCTAAAAACCGAGATTGACACCTTTAACGGCTAAATAACCTg                                                                                                                                          >  8:271661/1‑136 (MQ=255)
 ttttCGTATCTCCATAAAGCGTAAACTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGca                                                                                                                                        <  3:289159/137‑1 (MQ=255)
      gTATCTCCATAAAGCGTAAACTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCAtctt                                                                                                                                 <  3:28972/139‑1 (MQ=255)
      gTATCTCCATAAAGCGTAAACTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCAtctt                                                                                                                                 <  8:209043/139‑1 (MQ=255)
      gTATCTCCATAAAGCGTAAACTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCAtctt                                                                                                                                 <  7:21945/139‑1 (MQ=255)
         tctcCATAAAGCGTAAACTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCAtcttct                                                                                                                               <  4:167685/138‑1 (MQ=255)
                    cGTAAACTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCgt                                                                                                                     >  6:16613/1‑137 (MQ=255)
                       aaaCTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCgtgagt                                                                                                                 >  4:153072/1‑138 (MQ=255)
                       aaaCTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCgtgagt                                                                                                                 >  6:169398/1‑138 (MQ=255)
                             tgAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCGCTTCACTGACATCGAAAATACGGAGATTGACAACTTTTTCGGCTAAATAAGCTGCATCTTCTTCTGTTTCATCGTGAGTGATGcc                                                                                                           >  8:86421/1‑138 (MQ=255)
                             tgAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGcc                                                                                                           >  2:180087/1‑138 (MQ=255)
                              gagaCGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGa                                                                                                         >  8:30993/1‑139 (MQ=255)
                              gagaCGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTAATCGGCTAAATAAGCTGCAACATCTTCTGTATCATCGGGAGTGATGCCGa                                                                                                         >  5:187029/1‑139 (MQ=255)
                               agaCGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGAc                                                                                                        <  7:30993/139‑1 (MQ=255)
                                  cGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     <  5:169398/139‑1 (MQ=255)
                                  cGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     <  5:16613/139‑1 (MQ=255)
                                  cGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     <  5:11545/139‑1 (MQ=255)
                                  cGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     <  1:180087/139‑1 (MQ=255)
                                  cGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     <  3:153072/139‑1 (MQ=255)
                                  cGGATAAAATCTCTCCTCCTATATCTACAAGAGACACATTCATTTTCCCTTCACTGTCATCGAAAATACTGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTAGCATCGTGAGTGATGCCGACCAg                                                                                                     <  5:43513/139‑1 (MQ=255)
                                     aTAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCAc                                                                                                  >  2:35704/1‑139 (MQ=255)
                                               tcctccTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     >  3:283338/1‑126 (MQ=255)
                                               tcctccTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     <  4:283338/126‑1 (MQ=255)
                                               tcctccTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCt                                                                                        >  3:239925/1‑139 (MQ=255)
                                               acctcctATATCTACAAGATACAGTTTCATTTTCCCTTCACTGTGATCGAAAATACGGAGATTGACAACTATATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGAGATGCCGACCAGCACCATAAGGc                                                                                          <  7:86421/136‑1 (MQ=255)
                                                 ctcctATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGcc                                                                                         <  7:297099/136‑1 (MQ=255)
                                                 ctcctATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGcc                                                                                         <  2:5308/136‑1 (MQ=255)
                                                 ctcctATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTg                                                                                      <  7:271661/139‑1 (MQ=255)
                                                    ctATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTcc                                                                                   <  4:239925/139‑1 (MQ=255)
                                                       tatCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAAt                                                                                >  4:312674/1‑139 (MQ=255)
                                                       tatCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAAt                                                                                >  8:22164/1‑139 (MQ=255)
                                                         tCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGcatc                                                                                                                                      <  4:104888/83‑1 (MQ=255)
                                                         tCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGcatc                                                                                                                                      >  3:104888/1‑83 (MQ=255)
                                                         tCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAAttt                                                                              >  3:33674/1‑139 (MQ=255)
                                                           tACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTgg                                                                            >  8:195094/1‑139 (MQ=255)
                                                             cAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGcc                                                                          >  5:111463/1‑139 (MQ=255)
                                                             cAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGcc                                                                          >  1:90505/1‑139 (MQ=255)
                                                                   aCAGATTCATTTTCCCTTCACTGTCATAGAAAATACGGAGATTGACAACTTTATAGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATATGGCCgacgac                                                                    >  5:102352/1‑139 (MQ=255)
                                                                    cAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGAcc                                                                   <  3:312674/139‑1 (MQ=255)
                                                                       aTTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACctct                                                                >  8:7788/1‑139 (MQ=255)
                                                                       aTTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACctct                                                                >  5:225632/1‑139 (MQ=255)
                                                                                          tCATCTAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACgct                                               <  3:40446/137‑1 (MQ=255)
                                                                                          tCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATa                                                                                              <  3:232272/90‑1 (MQ=255)
                                                                                          tCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATa                                                                                              >  4:232272/1‑90 (MQ=255)
                                                                                           cATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     >  2:130688/1‑82 (MQ=255)
                                                                                           cATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     <  1:130688/82‑1 (MQ=255)
                                                                                           cATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATAATCTTCTGTATCATCGTGAGTCATGCCGACCAGCACCATAAGGCGTTGTCCAATTTGTCCGACGACCTCTTCGTCACCTTTAACgcttgc                                            >  2:27612/1‑139 (MQ=255)
                                                                                                      cGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGAcc                                                                   >  6:67810/1‑105 (MQ=255)
                                                                                                      cGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGAcc                                                                   <  5:67810/105‑1 (MQ=255)
                                                                                                         agaTTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     >  4:76415/1‑68 (MQ=255)
                                                                                                         agaTTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAg                                                                                                     <  3:76415/68‑1 (MQ=255)
                                                                                                            atGACAACTTTATCGGCTAAATAAGCTGCATGATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACAATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGaa                           <  6:102352/138‑1 (MQ=255)
                                                                                                                   cTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGAACAACTaa                    >  3:152913/1‑139 (MQ=255)
                                                                                                                     ttATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGaa                           <  7:161600/130‑1 (MQ=255)
                                                                                                                     ttATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGaa                           >  8:161600/1‑130 (MQ=255)
                                                                                                                     ttATCGGCTAAATAAGCTACATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGTCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGaa                           >  6:5688/1‑130 (MQ=255)
                                                                                                                     ttATCGGCTAAATAAGCTACATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGaa                           >  6:9982/1‑130 (MQ=255)
                                                                                                                     ttATCGGCTAAATAAGCTACATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGaa                           <  5:9982/130‑1 (MQ=255)
                                                                                                                     ttATCGGCTAAATAAGCTACATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGaa                           <  5:5688/130‑1 (MQ=255)
                                                                                                                               aaTAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGAACAACTAATCTCTTTTCTaa        <  1:35704/139‑1 (MQ=255)
                                                                                                                                     cTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGAACAACTAATCTCTTTTCTAACCCCtt  <  6:111463/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
TTTTTCGTATCTCCATAAAGCGTAAACTGTGAGACGGATAAAATCTCTCCTCCTATATCTACAAGAGACAGATTCATTTTCCCTTCACTGTCATCGAAAATACGGAGATTGACAACTTTATCGGCTAAATAAGCTGCATCATCTTCTGTATCATCGTGAGTGATGCCGACCAGCACCATAAGGCCTTGTCCAATTTGGCCGACGACCTCTTCGTCCACTGTAACGCTTGCTTCTGTTACTCGCTGAACAACTAATCTCTTTTCTAACCCCTT  >  BS168_reference_genome/2820258‑2820529

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