Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	2,865,943	A→T	8.2%	F221Y (TTC→TAC)	folC ←	folyl‑polyglutamate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,865,943	0	A	T	8.2%	143.8 / 4.6	49	F221Y (TTC→TAC)	folC	folyl‑polyglutamate synthase
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base T (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

TCTTTGATGGGTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTGCAATTTCTTCAATGGTGTTTCCCAAA  >  BS168_reference_genome/2865813‑2866080                                                                                                                                   |                                                                                                                                          tCTTTGATGGGTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCgca                                                                                                                                   <  4:95347/139‑1 (MQ=255)    ttGATGGGTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCgcagca                                                                                                                                >  3:81259/1‑139 (MQ=255)    ttGATGGGTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCgcagca                                                                                                                                >  4:257058/1‑139 (MQ=255)           gTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGTACGGCGCAGCATGCCGtt                                                                                                                         >  5:170583/1‑139 (MQ=255)           gTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGTACGGCGCAGCATGCCGtg                                                                                                                         >  5:174629/1‑138 (MQ=255)           gTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCTGCATGCCGtt                                                                                                                         >  7:125393/1‑139 (MQ=255)           gTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGtt                                                                                                                         <  4:222614/139‑1 (MQ=255)           gTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGtt                                                                                                                         <  1:23416/139‑1 (MQ=255)           gTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGtt                                                                                                                         >  7:48601/1‑139 (MQ=255)                      agaTGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATgg                                                                                                              >  8:165158/1‑139 (MQ=255)                           tCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         >  2:237358/1‑139 (MQ=255)                           tCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         <  2:226108/139‑1 (MQ=255)                           tCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         <  2:78109/139‑1 (MQ=255)                           tCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         <  3:257058/139‑1 (MQ=255)                           tCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         >  2:12010/1‑139 (MQ=255)                           tCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         <  1:222131/139‑1 (MQ=255)                             cGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGa                                                                                                       <  3:106469/139‑1 (MQ=255)                                atatCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCt                                                                                                    <  4:60092/139‑1 (MQ=255)                                             aTTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCgg                                                                                       <  8:48601/139‑1 (MQ=255)                                               ttCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         >  8:152260/1‑119 (MQ=255)                                               ttCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         <  7:152260/119‑1 (MQ=255)                                                      ccGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGt                                                                                 <  7:165158/136‑1 (MQ=255)                                                         ttttGAATGAAAACTGTTCTCCTGCAGGCAGAGCCACTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAActgc                                                                           >  7:273620/1‑139 (MQ=255)                                                              tatgaaAACTGTTCTACTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCTTCATTCAAAGACTGGAACGGCGCATCATGCCGTTCGGCTTCATGGAGGATTACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTa                                                                       <  8:288191/137‑1 (MQ=255)                                                                  aaaaCTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGa                                                                   <  3:34997/138‑1 (MQ=255)                                                                      cTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGTACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGa                                                              <  6:174629/139‑1 (MQ=255)                                                                      cTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGTACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGa                                                              <  6:170583/139‑1 (MQ=255)                                                                            tcCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCgg                                                                                       >  6:233683/1‑108 (MQ=255)                                                                            tcCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCgg                                                                                       <  5:233683/108‑1 (MQ=255)                                                                              cTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATAcc                                                         <  1:12010/136‑1 (MQ=255)                                                                              cTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACcttc                                                      >  8:281837/1‑139 (MQ=255)                                                                                   ggcagAGCCGCTTCATTGAACAAAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTtaa                                                 <  8:125393/139‑1 (MQ=255)                                                                                        agCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCTGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATg                                            >  3:71882/1‑139 (MQ=255)                                                                                        agCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGt                                                                                 >  2:244634/1‑102 (MQ=255)                                                                                        agCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGt                                                                                 <  1:244634/102‑1 (MQ=255)                                                                                        agCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATg                                            >  4:168785/1‑139 (MQ=255)                                                                                             ctTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACg                                                                    >  6:123240/1‑110 (MQ=255)                                                                                             ctTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACg                                                                    <  5:123240/110‑1 (MQ=255)                                                                                                       aaTAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCt                             >  2:155630/1‑139 (MQ=255)                                                                                                        aTAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTg                            >  5:135447/1‑139 (MQ=255)                                                                                                            cacaTGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTGCAAt                        <  4:127270/139‑1 (MQ=255)                                                                                                               aTGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGTAGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTCTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTGCAATttc                     >  2:89194/1‑139 (MQ=255)                                                                                                                       tGCAATGACTTGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTGCAATTTCTTCAATgg             >  2:183412/1‑139 (MQ=255)                                                                                                                       tGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         >  4:81801/1‑47 (MQ=255)                                                                                                                       tGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGAt                                                                                                         <  3:81801/47‑1 (MQ=255)                                                                                                                       tGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTGCAATTTCTTCAATgg             >  2:14302/1‑139 (MQ=255)                                                                                                                                  gAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTGCAATTTCTTCAATGGTGTTTCCaaa   >  7:74012/1‑138 (MQ=255)                                                                                                                                  gAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTGCAATTTCTTCAATGGTGTTTCCCaaa  <  7:281837/139‑1 (MQ=255)                                                                                                                                  gAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATAGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTCCAATTTCTTCATTGGAGTATCCCaa   >  1:143582/1‑138 (MQ=255)                                                                                                                                   |                                                                                                                                          TCTTTGATGGGTGCCAATCAGAGATGTCCGAATATCTTCATAGCATTTCTCTTCCGTTTTGAATGAAAACTGTTCTCCTGCAGGCAGAGCCTCTTCATTGAAAATAACACATGCATCATGCAATGACTGGAACGGCGCAGCATGCCGTTCGGCTTCATGGCGGATGACCTGTAAAGCTTCCGGCTGGGTAACTGCTGTAACGATTGGAATACCTTCTTTAATAATGCCGGCCTTTTCTCCTGCAATTTCTTCAATGGTGTTTCCCAAA  >  BS168_reference_genome/2865813‑2866080

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