Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,023,150 C→A 34.8% intergenic (+50/‑200) yhdF → / → bcaP putative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,023,1500CA34.8% 82.0 / 8.0 51intergenic (+50/‑200)yhdF/bcaPputative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter
Reads supporting (aligned to +/- strand):  ref base C (9/21);  new base A (10/6);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 5.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAC  >  BS168_reference_genome/1023024‑1023279
                                                                                                                                 |                                                                                                                                  
tgctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGGTTTATTTTAAAGTAATAAACGAAAAAACCCCGAAAGAGCGGGTTTTTTTTGTGGTGGAAAATTTGAAAAAAAAAAAt                                                                                                                        >  7:174690/1‑139 (MQ=255)
tgctgGCGTCTGACGAATCATCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTTAACGTTATGAACGAAAAACCAACACAAGTGCTGGGGTTTTGTTTGGTGGAAAATTTGAAAATAAAAAAt                                                                                                                        >  6:157662/1‑139 (MQ=255)
  ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTTAACGTAATAAACCAAAAAAAAGCTCAATAGCGGGTTTTTTTTTTTGTTTTAAAATTAAAAATAAAAAAatt                                                                                                                      >  3:103233/1‑139 (MQ=255)
  ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTTAACCTAAACAAACAAAAAACACCACAAAACGGGGGGTTTTGTGTGGGGTTAAAATTTAAAAAAAAAAAAta                                                                                                                       >  7:218339/1‑138 (MQ=255)
               aaTCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAACGTAATCAACGAAAAAACACCTCAAGAGCGGGGTTTTTGTGTGGGGTAAAATTTAAAAATAAAAAAATTTTGGAttttttt                                                                                                         >  8:61132/1‑136 (MQ=255)
                 tCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTGAACGGAATAAACCAAAAACCACCACAACATCTGGTGTTTTTTGTGGTGTTAAATTTAAAAATAAAAAAATTTTAGAtttttttta                                                                                                       >  4:182144/1‑134 (MQ=255)
                               aGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAAAGAAAAAACACCACAAAAGCTGGGTTTTTGTGTGGTGTAAAATTTAAAAAAAAAAAAATTTTCGATTTTTGTAAAAATTTTTTTTTTAt                                                                                        >  7:227034/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTTAAAGTTATAAAACAAAAACCAACACAAAAGCAGGGTTTTTGTGTGGTGAAAAAATTAAAAAAAAAAAAATTTTACATTTTTTTAAAAAAttttttttaaatt                                                                                      >  6:62474/1‑134 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTTTTTAAACGTAATAAACGAAAAAACACCTCAAGCGCTGGTTTTTTTTGTGGTGAAAAATTTAAAAATAAAAAATTTTTGGATTTTTATAAAAATTTTTTTTTGAttt                                                                                      >  5:96073/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAAAGTAATAAACCAAAAAACCGCTCAAGAGCGGGGTTTTTGTGTGGTGTAAAATTTTAAAATAAAAAATTTTTGGGTTTTTGTAACAATttttttttgtttt                                                                                      >  7:150898/1‑134 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTTAAAGTTATAAAAAAAAAAACACCACAAAAGCAGGTTTTTTGTTGTGGTTAAAATTTTAAAAAAAAAAATTTTTTGTTTTTTTTAAAAAATTTTTTTTTAttt                                                                                      >  4:54265/1‑139 (MQ=255)
                                 ggCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTAAAAATAATAAAACAAAAAACACCACAAGAGCAGGGTTTTTTTTTTGTGTAAAATTTAAAAAAAAAAAAATTTTAGATTTTTATTAAAAATTTTTTTTTAttt                                                                                      >  4:188081/1‑139 (MQ=255)
                                   cAGACCATTCATGTGAATGGCGGCCGTTTTTTTTGAACGTATTAAACGAAAAAACACCGCAAGACGGGTTTTTTTGTTTGGTGAAAAATTTTAAAATAAAAAAATATTGGATTTTTTTTACAAATTTTTTTTGAttttt                                                                                    >  6:251464/1‑139 (MQ=255)
                                    aGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                   <  1:93730/139‑1 (MQ=255)
                                      aCCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGTAAAATTTCAAAATAAAAAATTTTTCGATTTTTATAACAATTTTTTTTTGtttttttg                                                                                  >  6:246226/1‑136 (MQ=255)
                                      aCCATTCATGTGAATGGCGGCCGTGTTTTTTTAACGTAAACAAACAAAAAACACCGCAAGAGAGGGTGTTTTTGTTTGGGTAAAATTTAAAAAAAAAAAAATTTTGTGATTTTTTTTATAATTTTTTTTTAttttttt                                                                                   >  2:113702/1‑137 (MQ=255)
                                         aTTGATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                              <  8:218339/139‑1 (MQ=255)
                                         aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                              <  3:188081/139‑1 (MQ=255)
                                         aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                              <  1:22617/139‑1 (MQ=255)
                                         aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                              <  6:96073/139‑1 (MQ=255)
                                         aTTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtt                                                                              <  1:113702/139‑1 (MQ=255)
                                          ttCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAttt                                                                             <  3:54265/139‑1 (MQ=255)
                                            cATGTGAATGGCGGCCGTGTTTTTTAAAACTAATAAAACAAAAAACATCACAAGAGCAGGGTTTTTTTGTGGTTTAAAATTTCAAAAAAAAAAAATTGTGCAATTTTTTTAAAATTTAGTTTTTAATTTTTAAAtttt                                                                            >  8:164/1‑138 (MQ=255)
                                             aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTTGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                          <  5:7780/139‑1 (MQ=255)
                                             aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                          <  2:37235/139‑1 (MQ=255)
                                             aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                          <  3:33731/139‑1 (MQ=255)
                                                tgAATGGCGGCCGTTTTATTTCAACGTAATCAAAGAAAAACCAGCTCAAGAGCTGGGTTTTTGTGTGGTGTTAAATTTAAAAATAAAAAATTTTTTGATTTTTATAAAAATTTTTTTTTGATTTTTTAAATTTtatatt                                                                       >  3:225679/1‑138 (MQ=255)
                                                 gAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                       <  3:264206/138‑1 (MQ=255)
                                                      gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                  <  3:181493/138‑1 (MQ=255)
                                                          ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                             <  5:62474/139‑1 (MQ=255)
                                                            gTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTa                                                           <  7:189432/139‑1 (MQ=255)
                                                            gTTTTATTTCAACGTAATCAACGAAAAACCACCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTa                                                           <  7:193968/139‑1 (MQ=255)
                                                                        cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAAt                                               >  3:156054/1‑139 (MQ=255)
                                                                            aTCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGAtt                                           >  8:103189/1‑139 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGCGCTGGTTTTTTGTGTGGTTAAAAATTTAAAAATAAAAATTTTGTCTATTTTTATAAAAATTTTTTTTTCATTTTTGAAATTTTATTTATATTTTTGAGTTGGATAAAGGATTGATTCGCGAcc                                    <  8:174690/139‑1 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTTAAAAATTTCAAAAAAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                    <  4:103233/139‑1 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTTAAAAATTTAAAAATAAAAAATTTGTCGATTTTTCTAACATTTTTGTTTTTATTTTGGCAATTTTATATTTAATATTGAGTAGTATAAAGTAATGATTCGCGAcc                                    <  2:61631/139‑1 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTTTTTTCATTTTTGCAATTTTATATATAATATTGGGTAGGACAAAGGAATGATTCGCGAcc                                    <  3:182144/139‑1 (MQ=255)
                                                                                   aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                    <  4:156054/139‑1 (MQ=255)
                                                                                    aaaaCCAGCTAACGAGCTGTTTTTTTGTGTGGTTAAAAATTTAAAAATAAAAAATTTGTGGATTTTTATAAAAATTTTTTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGTATAAAGTAATTATTCGCGAc                                     <  4:225679/137‑1 (MQ=255)
                                                                                        ccAGCTCAAGAGATGGTTTTTTGTGTGGTTAAAAATTTCAAAATAAAAAATTTGTCGATTTTTATAACAATTTTGTTTTCATTTTTGAAATTTTATATATAATATTGAGTAGGATAAAGGAATGATTCGCGACCAGTAc                               <  8:150898/139‑1 (MQ=255)
                                                                                           gCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                   <  3:176050/84‑1 (MQ=255)
                                                                                           gCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTg                                                                                   >  4:176050/1‑84 (MQ=255)
                                                                                            cTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTTACAATTTTGTTTTGATTTTTGGAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAat                           >  3:398/1‑139 (MQ=255)
                                                                                                      ggTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGaa                  >  3:287252/1‑138 (MQ=255)
                                                                                                      ggTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTAATTTTTGCAATTTTATATATAATATTGATTAGGAGAAAGTAATGATTCGAGACCAGTACAAATATAGTATaa                  >  7:288476/1‑138 (MQ=255)
                                                                                                                 tgtgGTGAAAAATTTAAAAATAAAAAATTTGTCGATTTTTATAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaat        <  7:61132/139‑1 (MQ=255)
                                                                                                                   tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataa      >  8:41551/1‑139 (MQ=255)
                                                                                                                   tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaata       <  7:103189/138‑1 (MQ=255)
                                                                                                                         aaaaaTTTAAAAATAAAAAATTTGTCGTTTTTTATAACATTTTTTTTTTCATTTTTGAAATTTTATATATAATATTGATTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAc  <  1:179315/138‑1 (MQ=255)
                                                                                                                         aaaaaTTTAAAAAAAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGAAATTTTAAATATAATATTGATTAGGCTAAAGTAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAc  <  8:227034/138‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                  
TGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATGAC  >  BS168_reference_genome/1023024‑1023279

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: