Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,023,150 C→A 20.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,1500CA20.8% 89.8 / 2.2 39intergenic (+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/20);  new base A (4/3);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATG  >  BS168_reference_genome/1023026‑1023277
                                                                                                                              |                                                                                                                               
ctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTTAACGTTATCAAAGAAAAACAACCACAACAGCGGGTTTTTTTTGTGGTGATAAAATTAAAAAATAAAAAAtt                                                                                                                    >  5:60752/1‑139 (MQ=255)
ctgGCGGCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAAttt                                                                                                                   <  6:175793/139‑1 (MQ=255)
   gCGTCTGACGAATATTACTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTGTTATTTTAAAGAAATAAAACAAAAAACACCACAAGAGCGGGTTTTTTTTGTGGTGGAAAAATTTAAAAAAAAAAAAatttg                                                                                                                 >  1:134950/1‑137 (MQ=255)
              aTCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTaa                                                                                                     <  1:37031/139‑1 (MQ=255)
                               ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCAttt                                                                                    >  1:116673/1‑139 (MQ=255)
                                 cAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATCAACGAAAAAAAAGCTCAAGAGCTGGTTTTTTTTGTGGTGAAAAAATTCAAAATAAAAAATTTTTCGATTTTTATAAAAATTTTTTTTTAAttttt                                                                                  >  5:154448/1‑139 (MQ=255)
                                    aCCATTCATGTGAATGGCGGCCGTGTTTTTTGAAACGAAACAAAGAAAAACCAACCCAAACGCAGGTTTTTTTTGTGGTGATAAAATTCAAAAAAAAAAAATTTTAGATTTTTGTAAAAAATTTTTTTTGAttttttg                                                                                >  6:122631/1‑136 (MQ=255)
                                         tCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAAAGAAAAACCAACTCAAGAGCTGGGTTTTTGTGTGGGGGAAAATTTTAAAATAAAAAAATTTTCGATTTTTTTAAAAAATTTTTTTTAATTTTTTGAAtttt                                                                          >  4:127299/1‑139 (MQ=255)
                                         tCATGTGAATGGCGGCCGTGTTTTTTAAAAGTAATAAAACAAAAAAAAACGCAAAAGCGGGGTTTTTTTGTGGTGTAAAATTTTAAAATAAAAAAATTTTGGATTTTTTGAAAAAATTTTTTTTGATTTTTTGAAtttt                                                                          >  2:206040/1‑139 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                          <  1:206040/138‑1 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                          <  4:104809/138‑1 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                          <  2:134950/138‑1 (MQ=255)
                                          cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                          <  2:78977/138‑1 (MQ=255)
                                           aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                        <  8:219874/139‑1 (MQ=255)
                                           aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                        <  3:209577/139‑1 (MQ=255)
                                              tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                     <  3:127299/139‑1 (MQ=255)
                                              tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                     <  1:174737/139‑1 (MQ=255)
                                              tgAATGGCGGCCGTTTTATTTAAACGTAATAAAACAAAAAACAGCGAAAAAGCAGGGTTTTTTTGTGGTGTAAAAATTAAAAATAAAAAAATTTTCGATTTTTGTATAAAATTTTTTTTTATTTTTTAAAtttttttta                                                                     >  2:174737/1‑134 (MQ=255)
                                               gAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTTGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                     <  7:64448/138‑1 (MQ=255)
                                                aaTGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatatata                                                                   <  1:13494/139‑1 (MQ=255)
                                                  tGGCGGCCGTTTTATTTAAACGTAATAAACGAAAAAACAGCTCAATAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTTTCGATTTTTCTAACAATTTTTTTTTGATTTTTTGAATTTTTTTTataa                                                                  >  7:219874/1‑138 (MQ=255)
                                                  tGGCGGCCGTGTTATTTAAACGTAATAAAACAAAAACCAGCTCAAGAGCTGGTTTTTTTTTTGGTGAAAAATTTGAAAATTAAAAAATTTTGGATTTTTTTAACAAATTTTTTTTCATTTTTTAAATTTTTTATttaa                                                                  >  4:209577/1‑138 (MQ=255)
                                                        ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                           <  2:116673/139‑1 (MQ=255)
                                                        ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                           <  7:130706/139‑1 (MQ=255)
                                                              tATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGga                                                      <  5:122631/138‑1 (MQ=255)
                                                                          aTCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGAGTTTTGTAACAATTTTTTTTTAATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGAtt                                         >  5:32716/1‑139 (MQ=255)
                                                                                gAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                          <  1:103636/100‑1 (MQ=255)
                                                                                gAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                          >  2:103636/1‑100 (MQ=255)
                                                                                 aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                  <  3:21973/139‑1 (MQ=255)
                                                                                                    ggTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTGTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGTATGATTCGCGACCAGTACAAATATAGTAGaa                >  4:182858/1‑138 (MQ=255)
                                                                                                     gTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                >  5:223189/1‑89 (MQ=255)
                                                                                                     gTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                                <  6:223189/89‑1 (MQ=255)
                                                                                                        ttttttGTGGTTAAAAATTTCAAAATAAAAAATTTGTCTATTTTTCTAACATTTTTGTTTTAATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATTATTCGCGACCAGTACAAATATAGTAGaaaaaa            <  8:90253/138‑1 (MQ=255)
                                                                                                         tttgttttgtTAAAAATTTAAAAATAAAAAATTTTGCTATTTTTATAAAAATTTTTTTTTATTTTTTGCAATTTTATATATAATATTGAGTAGGACAAAGTAATCATTCGCGACCAGTACAAATATAGTAGaaaaaaa           <  1:98551/138‑1 (MQ=255)
                                                                                                             tgggttgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTTCAATTTTAAATATAATATTGAGTCGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaat      <  1:137843/134‑1 (MQ=255)
                                                                                                               tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGTAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGtataa    >  1:31642/1‑139 (MQ=255)
                                                                                                               tggtgAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataa    >  7:40915/1‑139 (MQ=255)
                                                                                                                tttttaaaaaTTTAAAAAAAAAAAATTTGTCTTTTTTTATAAAATTTTTGTTTTAATTTTTTAAATTTTATATATAATTTTGAGTAGTATAAAGTAATTATCCGCGACCAGTACAAATATAGTAGAAAAAAAAGaata     <  6:32716/133‑1 (MQ=255)
                                                                                                                  gttaaaaaTTTAAAAATAAAAAATTTTTCGATTTTTCTAAAATTTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGAATCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATg  <  1:34170/136‑1 (MQ=255)
                                                                                                                              |                                                                                                                               
CTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGATTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAATG  >  BS168_reference_genome/1023026‑1023277

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: