Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,958,235	C→A	8.5%	A282E (GCG→GAG)	ilvK →	branched‑chain amino acid aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,958,235	0	C	A	8.5%	149.6 / 4.0	47	A282E (GCG→GAG)	ilvK	branched‑chain amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base A (2/2);  total (25/22)
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.24e-01

AACATCTTCTTCAAAATCAACGGTGAAATCGTAACACCGATGCTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGATGAAACACTTTCGATCAACAACGGTGAAA  >  BS168_reference_genome/3958102‑3958369                                                                                                                                      |                                                                                                                                       aaCATCTTCTTCAAAATCAACGGTGAAATCGTAACACCGATGCTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGcgatcg                                                                                                                                   <  2:174777/139‑1 (MQ=255)         cttcAAAATCAACGGTGAAATCGTAACACCGATGCTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAgg                                                                                                                             <  3:21596/137‑1 (MQ=255)         cttcAAAATCAACGGTGAAATCGTAACACCGATGCTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAg                                                                                                                              <  5:296634/136‑1 (MQ=255)              aaaTCAACGGTGAAATCGTAACACCGATGCTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCa                                                                                                                      <  3:91004/139‑1 (MQ=255)              aaaTCAACGGTGAAATCGTAACACCGATGCTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCa                                                                                                                      <  6:253385/139‑1 (MQ=255)                                        aTGCTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATcc                                                                                                  >  8:129026/1‑133 (MQ=255)                                        aTGCTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATcc                                                                                                  <  7:129026/133‑1 (MQ=255)                                           cTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACTAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGgaagaa                                                                                          <  6:113296/138‑1 (MQ=255)                                              aaCGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACgg                                                                                                       <  5:62891/122‑1 (MQ=255)                                              aaCGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACgg                                                                                                       >  6:62891/1‑122 (MQ=255)                                              aaCGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCt                                                                                      >  7:276823/1‑139 (MQ=255)                                                   gAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAgg                                                                                                                             <  1:21490/95‑1 (MQ=255)                                                   gAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAgg                                                                                                                             >  2:21490/1‑95 (MQ=255)                                                   gAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAATGACGGCATCCTGGAAGAAGCCTTCGGa                                                                                 <  5:136748/139‑1 (MQ=255)                                                   gAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCgg                                                                                  >  2:40867/1‑138 (MQ=255)                                                   gAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGa                                                                                 <  4:220088/139‑1 (MQ=255)                                                             gAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTAc                                                                         <  1:40867/137‑1 (MQ=255)                                                                   aTTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGGCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTGCGGAACAGGTACAGCAGCTg                                                                 >  2:181797/1‑139 (MQ=255)                                                                      aCGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCtt                                                                                     >  3:140328/1‑116 (MQ=255)                                                                      aCGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCtt                                                                                     <  4:140328/116‑1 (MQ=255)                                                                       cgcgCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTAt                                                             >  1:73640/1‑139 (MQ=255)                                                                             aTTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGCACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTcccc                                                       <  4:218482/139‑1 (MQ=255)                                                                                   tCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCgg                                                 >  7:267080/1‑139 (MQ=255)                                                                                      tCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGAGATCGATGAgg                                                                                                                             >  7:163459/1‑60 (MQ=255)                                                                                      tCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGAGATCGATGAgg                                                                                                                             >  7:168017/1‑60 (MQ=255)                                                                                      tCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGAGATCGATGAgg                                                                                                                             <  8:163459/60‑1 (MQ=255)                                                                                      tCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGAGATCGATGAgg                                                                                                                             <  8:168017/60‑1 (MQ=255)                                                                                        gCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGc                                            >  3:152372/1‑139 (MQ=255)                                                                                             gCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCa                                                                                                                      >  1:277436/1‑60 (MQ=255)                                                                                             gCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCa                                                                                                                      <  2:277436/60‑1 (MQ=255)                                                                                             gCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGgaagaa                                                                                          >  1:107503/1‑88 (MQ=255)                                                                                             gCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGgaagaa                                                                                          <  2:107503/88‑1 (MQ=255)                                                                                                 aaGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTgg                                    >  6:299491/1‑138 (MQ=255)                                                                                                  aGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTgg                                    >  4:152520/1‑137 (MQ=255)                                                                                                     aTTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGa                               <  4:152372/139‑1 (MQ=255)                                                                                                       tGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGACGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGAt                              >  7:190635/1‑138 (MQ=255)                                                                                                        ggggCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGa                                                                                 >  7:11992/1‑86 (MQ=255)                                                                                                        ggggCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGa                                                                                 <  8:11992/86‑1 (MQ=255)                                                                                                        ggggCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGATg                             >  3:297538/1‑138 (MQ=255)                                                                                                              ttCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCgg                                                                                  <  2:310715/79‑1 (MQ=255)                                                                                                              ttCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCgg                                                                                  <  4:211686/79‑1 (MQ=255)                                                                                                              ttCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCgg                                                                                  >  3:211686/1‑79 (MQ=255)                                                                                                              ttCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGATGAACCACt                      >  1:94948/1‑139 (MQ=255)                                                                                                              ttCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGATGAAACACt                      >  3:194384/1‑139 (MQ=255)                                                                                                              ttCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAACACGGCATCCTGGAAGAAGCCTTCgg                                                                                  >  1:310715/1‑79 (MQ=255)                                                                                                               tCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGATGAAACACtt                     >  4:230214/1‑139 (MQ=255)                                                                                                                                  aTTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGATGAAACACTTTCGATCAACAACGGTGaaa  <  4:194384/139‑1 (MQ=255)                                                                                                                                    tGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGATGAAACACTTTCGATCAACAACGGTGaaa  >  1:80091/1‑137 (MQ=255)                                                                                                                                      |                                                                                                                                       AACATCTTCTTCAAAATCAACGGTGAAATCGTAACACCGATGCTGAACGGGAGCATCCTGGAAGGCATTACGCGCAATTCAGTCATCGCCTTGCTTAAGCATTGGGGCCTTCAAGTTTCAGAACGAAAAATTGCGATCGATGAGGTCATCCAAGCCCATAAAGACGGCATCCTGGAAGAAGCCTTCGGAACAGGTACAGCAGCTGTTATTTCCCCAGTCGGCGAGCTGATCTGGCAGGATGAAACACTTTCGATCAACAACGGTGAAA  >  BS168_reference_genome/3958102‑3958369

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: ‑