Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	616,094	C→T	7.1%	T75M (ACG→ATG)	ydhC →	putative transcriptional regulator (GntR family)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	616,094	0	C	T	7.1%	182.7 / 3.2	56	T75M (ACG→ATG)	ydhC	putative transcriptional regulator (GntR family)
Reads supporting (aligned to +/- strand):  ref base C (23/29);  new base T (2/2);  total (25/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 = 8.39e-01

GGGGAACGGATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACGAAGCGGACGAAGAAGAGCTGACCATACTCGAAAAACAGCTTGAAGAAACAGAGA  >  BS168_reference_genome/615964‑616231                                                                                                                                   |                                                                                                                                          ggggAACGGATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTaa                                                                                                                                   >  7:123471/1‑139 (MQ=255) ggggAACGGATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTaa                                                                                                                                   >  3:249038/1‑139 (MQ=255) ggggAACGGATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAAGGGGTTTTCGATCCCGTCGTTaa                                                                                                                                   >  8:22773/1‑139 (MQ=255)  gggAACGGATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTa                                                                                                                                    <  5:109761/137‑1 (MQ=255)  gggAACGGATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTa                                                                                                                                    <  3:288119/137‑1 (MQ=255)      aCGGATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCa                                                                                                                              <  5:59824/139‑1 (MQ=255)      aCGGATAAACGAAACACAGCTGGCGAAAAGCGTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCa                                                                                                                              <  6:36848/139‑1 (MQ=255)        ggATAAACGAAACACAGCTGGCGAAACGCTTTGGCGTCAGCCGCTCTCCGATAAGTGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGGCGTTAACGGCaaa                                                                                                                            <  1:411752/139‑1 (MQ=255)        ggATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCaaa                                                                                                                            <  5:187579/139‑1 (MQ=255)        ggATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCaaa                                                                                                                            <  2:292099/139‑1 (MQ=255)        ggATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCaa                                                                                                                             >  1:292099/1‑138 (MQ=255)         gATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCaa                                                                                                                             >  2:249274/1‑137 (MQ=255)          aTAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAg                                                                                                                          <  8:214646/139‑1 (MQ=255)             aaCGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATg                                                                                                                       >  4:91675/1‑139 (MQ=255)                  aaCACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGAt                                                                                                                  >  6:408558/1‑139 (MQ=255)                  aaCACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGAt                                                                                                                  <  7:61951/139‑1 (MQ=255)                            gCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGGTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTtat                                                                                                         <  8:181164/138‑1 (MQ=255)                                      ttGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAAt                                                                                              >  4:52930/1‑139 (MQ=255)                                              aGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTg                                                                                       >  2:348955/1‑138 (MQ=255)                                                   ctccCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGa                                                                                                                         <  1:266454/99‑1 (MQ=255)                                                   ctccCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGa                                                                                                                         >  2:266454/1‑99 (MQ=255)                                                           tAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGt                                                                         <  1:348955/139‑1 (MQ=255)                                                             aGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTc                                                                        >  4:390421/1‑138 (MQ=255)                                                                gAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTa                                                                                                                                    <  4:317604/75‑1 (MQ=255)                                                                gAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTa                                                                                                                                    >  3:317604/1‑75 (MQ=255)                                                                gAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    <  4:249038/139‑1 (MQ=255)                                                                gAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    <  3:390421/139‑1 (MQ=255)                                                                gAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    <  8:28996/139‑1 (MQ=255)                                                                gAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    <  3:91675/139‑1 (MQ=255)                                                                gAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    >  2:88704/1‑139 (MQ=255)                                                                gAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    <  8:123471/139‑1 (MQ=255)                                                                gAAGCCATGCGATTGCTTGAGAAAGATGGTTTGGTAATAGAGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAACATCAGAATTCCGCTGGAGCAGCTTGCAGTCGAAt                                                                    <  8:189756/139‑1 (MQ=255)                                                                  aGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGt                                                                         >  8:292752/1‑132 (MQ=255)                                                                  aGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGt                                                                         <  7:292752/132‑1 (MQ=255)                                                                  aGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    >  8:366374/1‑137 (MQ=255)                                                                  aGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    >  1:142494/1‑137 (MQ=255)                                                                        gCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCATGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    <  1:82738/131‑1 (MQ=255)                                                                        gCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCATGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    >  2:82738/1‑131 (MQ=255)                                                                        gCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCATGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    <  1:88990/131‑1 (MQ=255)                                                                        gCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCATGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    >  2:88990/1‑131 (MQ=255)                                                                         cGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCa                                                              >  7:102471/1‑136 (MQ=255)                                                                          gATTGCTTGAGATAGATTGTTTTTTAAAAGGGGATGACCTAATTGGGTTTTGGATCACGTCGTTAACGGCAAAAGATGTATAGCATATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGa                                                          <  2:113042/139‑1 (MQ=255)                                                                              gCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    <  2:69506/125‑1 (MQ=255)                                                                              gCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    >  1:69506/1‑125 (MQ=255)                                                                              gCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACg                                                        >  3:189575/1‑137 (MQ=255)                                                                                   agaAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACGAAGCGGa                                                 <  5:422481/139‑1 (MQ=255)                                                                                    gaAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACGAAGCGGa                                                 <  1:88704/138‑1 (MQ=255)                                                                                        gATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACGAAGCGGACgaag                                            <  6:200572/139‑1 (MQ=255)                                                                                         aTGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACGAAGCGGACgaag                                            >  6:331583/1‑138 (MQ=255)                                                                                                   aaaaGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCa                                                                                                                              <  7:315505/46‑1 (MQ=255)                                                                                                   aaaaGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCa                                                                                                                              >  8:315505/1‑46 (MQ=255)                                                                                                     aaGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACGAAGCGGACGAAGAAGAGCTGACCAt                               <  5:331583/139‑1 (MQ=255)                                                                                                                          tttCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    <  3:165626/81‑1 (MQ=255)                                                                                                                          tttCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAAt                                                                    >  4:165626/1‑81 (MQ=255)                                                                                                                          tttCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACGAAGCGGACGAAGAAGAGCTGCCCATACTCGAAAAACAGCTTgaag           >  6:7303/1‑138 (MQ=255)                                                                                                                                  aCGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACGAAGCGGACGAAGAAGAGCTGACCATACTCGAAAAACAGCTTGAAGAAACAgaga  <  4:189575/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                          GGGGAACGGATAAACGAAACACAGCTGGCGAAAAGCTTTGGCGTCAGCCGCTCTCCGATAAGGGAAGCCATGCGATTGCTTGAGAAAGATGGTTTGTTAAAAGCGGATGACCGAAATGGGTTTTCGATCACGTCGTTAACGGCAAAAGATGTAGATGAGATTTATAAAATCAGAATTCCGCTGGAGCAGCTGGCAGTCGAATTAGTCATAGACGAAGCGGACGAAGAAGAGCTGACCATACTCGAAAAACAGCTTGAAGAAACAGAGA  >  BS168_reference_genome/615964‑616231

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