Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	442,499	C→T	10.8%	T310M (ACG→ATG)	gabT →	4‑aminobutyrate aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	442,499	0	C	T	10.8%	106.2 / 3.2	37	T310M (ACG→ATG)	gabT	4‑aminobutyrate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (2/2);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

TTCGCAATTGAGCACTTTGATGTAGTGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGACAAGGCGTATGAGTGGAAACAAGAATTC  >  BS168_reference_genome/442363‑442635                                                                                                                                         |                                                                                                                                         ttCGCAATTGAGCACTTTGATGTAGTGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACg                                                                                                                                         >  2:106548/1‑138 (MQ=255)           aGCACTTTGATGTAGTGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGccg                                                                                                                              >  2:223977/1‑137 (MQ=255)             cACTTTGATGTAGTGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGcc                                                                                                                            >  6:208722/1‑139 (MQ=255)                  tGATGTAGTGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGTCGGCAGCCCGCt                                                                                                                        <  1:180699/138‑1 (MQ=255)                  tGATGTAGTGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCt                                                                                                                        <  7:75552/138‑1 (MQ=255)                          tGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCa                                                                                                                               >  7:263515/1‑123 (MQ=255)                          tGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCa                                                                                                                               <  8:263515/123‑1 (MQ=255)                          tGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTgcgc                                                                                                               >  1:146538/1‑139 (MQ=255)                          tGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTgcgc                                                                                                               >  5:132570/1‑139 (MQ=255)                          tGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTgcgc                                                                                                               >  7:82499/1‑139 (MQ=255)                           gCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGAGCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGcgcg                                                                                                              <  2:163266/139‑1 (MQ=255)                            cccGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCgcgg                                                                                                             >  8:214398/1‑139 (MQ=255)                                        aCGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCa                                                                                                  <  8:82499/138‑1 (MQ=255)                                        aCGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCa                                                                                                  <  1:223977/138‑1 (MQ=255)                                                          gCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCga                                                                                 <  7:155500/137‑1 (MQ=255)                                                             gCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCgaagaag                                                                            <  5:208722/139‑1 (MQ=255)                                                                            aGCGGTGTGAGCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAg                                                              <  8:34854/138‑1 (MQ=255)                                                                                gtgtGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATc                                                         >  1:45589/1‑139 (MQ=255)                                                                                 tgtgATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCa                                                                                                                               >  6:169577/1‑68 (MQ=255)                                                                                 tgtgATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCa                                                                                                                               <  5:169577/68‑1 (MQ=255)                                                                                 tgtgATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCGTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCt                                                        >  2:247576/1‑139 (MQ=255)                                                                                   tgATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTg                                                       >  8:70241/1‑138 (MQ=255)                                                                                         ggCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCGCGGCTGCGCGGCGGCGTTGGCAGTCGTGGATATTATCGAAGGAGCAGGACTGTATGAGCGATCTGCAGAAAt                                                >  6:278981/1‑139 (MQ=255)                                                                                               cGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCATGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTgg                                                                                                    <  5:75950/81‑1 (MQ=255)                                                                                               cGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCATGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTgg                                                                                                    >  6:80218/1‑81 (MQ=255)                                                                                               cGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCATGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTgg                                                                                                    >  6:75950/1‑81 (MQ=255)                                                                                               cGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCATGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTgg                                                                                                    <  5:80218/81‑1 (MQ=255)                                                                                                 gAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCGGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCaaaa                                        <  2:15544/139‑1 (MQ=255)                                                                                                            tGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGa                              >  2:243095/1‑138 (MQ=255)                                                                                                             gCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGACa                            <  1:243095/139‑1 (MQ=255)                                                                                                             gCGGCAGCGCAAGGGGAGCGGGGCGGCACGTATGCCGGCAGCCCGCGCGGCGGCGCGGCGGAGTTGGCCGTCTTTGCTATTATCGCACAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGACa                            <  7:166763/139‑1 (MQ=255)                                                                                                              cGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGTACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGACa                            <  6:132570/138‑1 (MQ=255)                                                                                                                               cTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATAGAAGAATAAGGACTGATTGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGACAAGGCGTATGAGTGGaaa           >  7:233042/1‑138 (MQ=255)                                                                                                                                  ggcggcACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGACAAGGCGTATGAGTGGAAACAAg       <  2:45589/139‑1 (MQ=255)                                                                                                                                       cACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGATAAGGCGTAt                    <  2:180594/121‑1 (MQ=255)                                                                                                                                       cACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGATAAGGCGTAt                    >  1:180594/1‑121 (MQ=255)                                                                                                                                       cACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGACAAGGCGTATGAGTGGAAACAAGAATTc  >  4:68752/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                         TTCGCAATTGAGCACTTTGATGTAGTGCCCGATCTGATCACGGTTTCTAAATCGCTTGCGGCAGGCCTGCCATTAAGCGGTGTGATCGGGCGTGCGGAAATGCTTGATGCGGCAGCGCCAGGGGAGCTGGGCGGCACGTATGCCGGCAGCCCGCTCGGCTGCGCGGCGGCTTTGGCAGTCTTGGATATTATCGAAGAAGAAGGACTGAATGAGCGATCTGAAGAAATTGGCAAAATCATTGAAGACAAGGCGTATGAGTGGAAACAAGAATTC  >  BS168_reference_genome/442363‑442635

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