Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 238,265 G→T 21.1% pseudogene (184/285 nt) ybfEc ← conserved protein of unknown function (C‑terminal end)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome238,2650GT21.1% 65.6 / 3.1 35pseudogene (184/285 nt)ybfEcconserved protein of unknown function (C‑terminal end)
Reads supporting (aligned to +/- strand):  ref base G (6/21);  new base T (6/2);  total (12/23)
Fisher's exact test for biased strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGAGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTAATGTCACAG  >  BS168_reference_genome/238129‑238402
                                                                                                                                        |                                                                                                                                         
atgaGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGTTCATACCGACTATGCGGGCGTGTTCGGATTTTTGTTTGAAGCGGTCAGAAATTTGa                                                                                                                                          >  7:9401/1‑138 (MQ=255)
atgaGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTGTTAGGTTTTTTTTTTTAAACGTTGAGAAATTTTaa                                                                                                                                         >  7:124191/1‑139 (MQ=255)
atgaGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGAGCGCGTGTGGAGGATTTTGTTTTTAGAGTTTAGAAATTTg                                                                                                                                           >  7:124187/1‑137 (MQ=255)
          aGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTGTTGGGTTTTTTGTTTTAAGCGGTCAGGAAATTTAACTTttgttg                                                                                                                                >  1:56480/1‑138 (MQ=255)
               aaGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTGTTGGGGGTTTTTTTTGTAGCGTTGTGGAATTTTAAGTTTTGGTGGAAATg                                                                                                                          >  3:192765/1‑139 (MQ=255)
                   aTCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGa                                                                                                                      <  5:97357/139‑1 (MQ=255)
                        tGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGaa                                                                                                                     <  5:275351/135‑1 (MQ=255)
                             tCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTGGTGGGGTTTTTGTTTGTAGCGGTGAGAAATTTGAACTTTTGTGGTAAATGGGGGTACAATaaa                                                                                                             >  8:155347/1‑138 (MQ=255)
                                                cTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTGTTGGGTTTTTTGTTTTAATCGTTCAGAAATTTGAACTTTTATTTTAAAAGGGGAAAAAATAAAAACCCAAAAAATGCCGAta                                                                                          >  7:140067/1‑137 (MQ=255)
                                                     tttGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTGTTGGGAGTTTTGTTTGAAGCGTTCAGAAATTTTAACTTTTCTTGAAAATGGGGGAACAATAAAAACCCAAAAAAGGCCGATCAAGGc                                                                                     >  7:154421/1‑138 (MQ=255)
                                                         tATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCaaaa                                                                                 <  1:144671/138‑1 (MQ=255)
                                                         tATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCaaaa                                                                                 <  3:120487/138‑1 (MQ=255)
                                                                  tCTAATACGAATTCGATCATACCTACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTTAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGt                                                                        <  1:136725/138‑1 (MQ=255)
                                                                                    aTACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAAt                                                      <  4:28793/138‑1 (MQ=255)
                                                                                          aCTATGCGGGCGTGTTCGGATTTTTTTTTGAAGCGTTGAGAAATTTTAACTTTTGTTGTAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACAATTGTACAGCTAATTATattaattaatta                                                >  7:128211/1‑138 (MQ=255)
                                                                                                cGGGCGTGTTGGGGGTTTTTTTTGAAGCGGTCAGAAAATTGAACTTTTCTTGAAAATTGGGAAACAAAAAAAACCCAAAAAATGCCGCTCAAGGCAAAAAACAATTGTTCAGCGAATTATATTTATTAATTAAAACGAg                                         >  8:207918/1‑139 (MQ=255)
                                                                                                  ggCGTTTTGGGATTTTTGTTTGAAGCGTTAAGAAATTTTAAATTATATTGAAAAAGGGGAAACAAAAAAAACCCAAAAAAGGCCGCACAAGGCAAAAAACAATTTTACATCTAATTCTATTAATTCATTAAAAAGAATg                                       >  8:239488/1‑139 (MQ=255)
                                                                                                              tttttGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAg                           <  5:29446/139‑1 (MQ=255)
                                                                                                               ttttGTTTGAAGCGTTCATAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAg                           <  1:127422/138‑1 (MQ=255)
                                                                                                                 ttgtttgAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGaa                         <  2:170260/138‑1 (MQ=255)
                                                                                                                    tttgAAGCGTTCAGAAATTTGAACTTTTCATGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGta                      <  4:24455/138‑1 (MQ=255)
                                                                                                                    tttgAAGCCTTCAGAAATTTGAAATTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGtat                     <  5:197687/139‑1 (MQ=255)
                                                                                                                       gAAGCGTTCATAAATTTGAACTTTTCTTGAAATTGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATaa                   <  7:129083/138‑1 (MQ=255)
                                                                                                                        aaGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACg                 <  7:239488/139‑1 (MQ=255)
                                                                                                                            gTTCAGAAATTTTAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGttat             <  7:6882/139‑1 (MQ=255)
                                                                                                                             ttCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAgt                                                                                                <  8:278380/55‑2 (MQ=255)
                                                                                                                             ttCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAgg                                                                                                >  7:278380/1‑55 (MQ=255)
                                                                                                                               cAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGttatta           <  4:275454/138‑1 (MQ=255)
                                                                                                                                aGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTaa          <  4:251349/138‑1 (MQ=255)
                                                                                                                                    aTTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTAATGTca     <  3:137316/139‑1 (MQ=255)
                                                                                                                                    aTTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTAATGTca     <  8:124191/139‑1 (MQ=255)
                                                                                                                                    aTTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTAATGTca     <  8:128211/139‑1 (MQ=255)
                                                                                                                                    aTTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTAATGTca     <  3:53415/139‑1 (MQ=255)
                                                                                                                                    aTTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTAATGTca     <  8:9401/139‑1 (MQ=255)
                                                                                                                                    aTTTGAACTTTTCATGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTAATGTca     <  8:124187/139‑1 (MQ=255)
                                                                                                                                       tGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTAATGTCACAg  >  1:248594/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                         
ATGAGGACAAAGGGAAAGCATCTCTGATTTCGAGATTATTCAGAGATGCTTTTTTTGTATTTTTGTTCTAATACGAAGTCGATCATACCGACTATGCGGGCGTTTTCGGATTTTTGTTTGAAGCGTTCAGAAATTTGAACTTTTCTTGAAAATGGGGAAACAATAAAAACCCAAAAAAGGCCGATCAAGGCAAAAAACATTTGTACAGCTAATTCTATTAATTCATTAAAACGAGTGCTGAAAAAATAGAAGTATAACGTTATTAATGTCACAG  >  BS168_reference_genome/238129‑238402

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: