Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,345,064 C→T 6.1% G117G (GGG→GGA nth ← endonuclease III, apurinic apyrimidic DNA lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,345,0640CT 6.1% 211.7 / 2.7 65G117G (GGG→GGAnthendonuclease III, apurinic apyrimidic DNA lyase
Reads supporting (aligned to +/- strand):  ref base C (35/26);  new base T (2/2);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

TTTTCAACTTCCAGAACCGAGTCCTTCCACCGGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCGATTGATTTAATATC  >  BS168_reference_genome/2344927‑2345189
                                                                                                                                         |                                                                                                                             
ttttCAACTTCCAGAACCGAGTCCTTCCACCGGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGAccc                                                                                                                              >  3:88861/1‑139 (MQ=255)
      aCTTCCAGAACAGAGTCCTTCCACCGGCAAATGCCCATTCGATTGCTGACTCTCTCCTCATGTGTATCTACGGGTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCa                                                                                                                        <  8:128509/139‑1 (MQ=255)
        atccccaaCCCAGCCCATCCCCCCGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGCGTATCTACGGCTAAGGCCGGCACGCCAAACCCAACGGAAACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGt                                                                                                                      <  6:296765/133‑1 (MQ=255)
          ccAGAACCGAGTCCTTCCACCGGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGttt                                                                                                                    >  5:296765/1‑139 (MQ=255)
             gAACCGAGTCCTTCCACCGGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGAc                                                                                                                 <  4:106241/139‑1 (MQ=255)
                   aGTCCTTCCACCGGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCt                                                                                                            >  6:304258/1‑138 (MQ=255)
                          ccACCGGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATcgc                                                                                                      <  8:276938/137‑1 (MQ=255)
                          ccACCGGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATcgc                                                                                                      <  2:244015/137‑1 (MQ=255)
                          ccACCGGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCgg                                                                                                    >  7:276938/1‑139 (MQ=255)
                                                            tcCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCaataa                                                                  >  1:50776/1‑139 (MQ=255)
                                                              cacaTGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCgg                                                                                            <  7:290260/111‑1 (MQ=255)
                                                              cacaTGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCgg                                                                                            >  8:290260/1‑111 (MQ=255)
                                                              cacaTGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCaataa                                                                  <  4:161872/137‑1 (MQ=255)
                                                                    ggTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTATCCCGACCCCTGGCAGTATGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCAtatt                                                           >  2:113204/1‑138 (MQ=255)
                                                                        tCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTc                                                       <  7:126594/138‑1 (MQ=255)
                                                                               cTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCtttt                                                >  2:251362/1‑138 (MQ=255)
                                                                                 aTGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAAACTTCAATAATCATTTTAATAAGCATTTTaa                                             >  2:158240/1‑139 (MQ=255)
                                                                                  tGGCCGGCACGCCAAACGCAACGGATCCCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGaa                                             >  2:20667/1‑138 (MQ=255)
                                                                                  tGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCgg                                                                                            >  4:36517/1‑91 (MQ=255)
                                                                                  tGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCgg                                                                                            <  3:36517/91‑1 (MQ=255)
                                                                                  tGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAat                                            >  2:148238/1‑139 (MQ=255)
                                                                                  tGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATATTAATAATCTTTTGaa                                             >  2:89116/1‑138 (MQ=255)
                                                                                  tGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATATTCAATAATCATTAAACTCAGCTTTTGaa                                             >  8:130974/1‑138 (MQ=255)
                                                                                    gCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGCCCCCTGGCAGTTTGCCAAGCTCATCGCGTTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAatat                                          >  5:68982/1‑139 (MQ=255)
                                                                                    gCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCCCTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAatat                                          >  2:52749/1‑139 (MQ=255)
                                                                                         cACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTccg                                                                                  >  2:163626/1‑94 (MQ=255)
                                                                                         cACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTccg                                                                                  <  1:163626/94‑1 (MQ=255)
                                                                                         cACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCaataa                                                                  >  6:98157/1‑110 (MQ=255)
                                                                                         cACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCaataa                                                                  <  5:98157/110‑1 (MQ=255)
                                                                                         cACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTc                                     >  3:49769/1‑139 (MQ=255)
                                                                                         cACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTc                                     >  1:250422/1‑139 (MQ=255)
                                                                                         cACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCGTTTGAATATTTTTc                                     >  5:232226/1‑139 (MQ=255)
                                                                                         cACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAAACTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTc                                     >  4:88294/1‑139 (MQ=255)
                                                                                                  cggaaCGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTttatta                            <  2:175735/136‑1 (MQ=255)
                                                                                                  cGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTttatta                            <  4:59171/139‑1 (MQ=255)
                                                                                                  cGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTttatta                            <  4:88861/139‑1 (MQ=255)
                                                                                                  cGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTttatta                            <  1:84353/139‑1 (MQ=255)
                                                                                                  cGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTttatta                            <  5:304258/139‑1 (MQ=255)
                                                                                                  cGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTttatta                            <  2:250422/139‑1 (MQ=255)
                                                                                                      aCGGATACCACTACGTTTGCGGTCTTTCTCCCGACTCCTGGCAGTTTGAc                                                                                                                 >  5:28580/1‑50 (MQ=255)
                                                                                                      aCGGATACCACTACGTTTGCGGTCTTTCTCCCGACTCCTGGCAGTTTGAc                                                                                                                 >  5:32775/1‑50 (MQ=255)
                                                                                                      aCGGATACCACTACGTTTGCGGTCTTTCTCCCGACTCCTGGCAGTTTGAc                                                                                                                 <  6:28580/50‑1 (MQ=255)
                                                                                                      aCGGATACCACTACGTTTGCGGTCTTTCTCCCGACTCCTGGCAGTTTGAc                                                                                                                 <  6:32775/50‑1 (MQ=255)
                                                                                                          aTACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATAttttt                                      >  2:230137/1‑121 (MQ=255)
                                                                                                          aTACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATAttttt                                      <  1:230137/121‑1 (MQ=255)
                                                                                                          aTACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGAtata                     >  3:43220/1‑138 (MQ=255)
                                                                                                             ccACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAAcc                 >  7:115886/1‑139 (MQ=255)
                                                                                                                cTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCgat              <  6:244403/139‑1 (MQ=255)
                                                                                                                cTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCgat              <  1:158240/139‑1 (MQ=255)
                                                                                                                 tACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCgatt             >  2:203902/1‑139 (MQ=255)
                                                                                                                 tACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCgatt             >  2:269893/1‑139 (MQ=255)
                                                                                                                   cGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGGTTTATTACGATATTATCCgattga           >  6:122334/1‑139 (MQ=255)
                                                                                                                   cGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCgat              <  8:35390/136‑1 (MQ=255)
                                                                                                                   cGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCgat              >  7:35390/1‑136 (MQ=255)
                                                                                                                   cGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCGTTTGAATATTTTTCGCTTTATTACGATATAAACCgattga           >  6:297537/1‑139 (MQ=255)
                                                                                                                        gCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCaataa                                                                  >  8:137041/1‑79 (MQ=255)
                                                                                                                        gCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCaataa                                                                  <  7:137041/79‑1 (MQ=255)
                                                                                                                        gCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCGATTGATTTaa      >  8:74077/1‑139 (MQ=255)
                                                                                                                        gCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCGATTGATTTaa      >  4:48618/1‑139 (MQ=255)
                                                                                                                        gCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCGCCGCCGGAATCTTCAATAATCATTTTACTCAGCGTTTGAATATTTTTAGCTTTATTACGATATAAACAGATTGATTTaa      >  3:215702/1‑139 (MQ=255)
                                                                                                                          ggTCTTTCTCCAGACCCCTGGCAGTTTGACAAGCTCATCTCGGTCTCTCTGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCGATTGATTTAata    <  7:130974/139‑1 (MQ=255)
                                                                                                                            tCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCGATTGATTTAATATc  <  1:203902/139‑1 (MQ=255)
                                                                                                                             cTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTtct                                                                                     <  2:285267/55‑1 (MQ=255)
                                                                                                                             cTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTtct                                                                                     >  1:285267/1‑55 (MQ=255)
                                                                                                                                      gACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCgat              >  4:136182/1‑117 (MQ=255)
                                                                                                                                      gACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCgat              <  3:136182/117‑1 (MQ=255)
                                                                                                                                         |                                                                                                                             
TTTTCAACTTCCAGAACCGAGTCCTTCCACCGGCAAATGCCCAATCGTTTGCTGACTCTCTCCACATGGGTATCTACGGCTATGGCCGGCACGCCAAACGCAACGGATACCACTACGTTTGCGGTCTTTCTCCCGACCCCTGGCAGTTTGACAAGCTCATCGCGGTCTCTCGGCACTTCTCCGCCGTAATCTTCAATAATCATTTTACTCAGCTTTTGAATATTTTTCGCTTTATTACGATATAAACCGATTGATTTAATATC  >  BS168_reference_genome/2344927‑2345189

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: