Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,493,138	G→A	P39P (CCC→CCT)	frc ←	formyl‑CoA transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,493,138	0	G	A	93.8%	93.2 / ‑0.8	32	P39P (CCC→CCT)	frc	formyl‑CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (13/17);  total (14/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 = 2.22e-01

GTGGATTGATTTCTTGAATATGCTCCCAGGTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAATTTAGCGTTTATTTAAGGCGTAATATCTTCTATAAATAACCCGCAAAATATCGTCGTAGCTGCACTCGATATTTTTAAAATCCTGTATCAAGT  >  NC_000913/2492903‑2493381                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    gTGGATTGATTTCTTGAATATGCTCCCAGGTGAAGCCCCTGTGATCAATGGCACCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATCTCTCGCAGCTCGTGACGCGTTACGTCACCAACGCCAGGCCGTTCAAtttt                                                                                                                                                                                                                                       >  1:224564/1‑250 (MQ=255)                gAATATGCTCCCAGGTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCATCTTTTCCATTTCCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCa                                                                                                                                                                                                                       <  2:212830/251‑1 (MQ=255)                 aaTATGCTCCCAGGTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCa                                                                                                                                                                                                                       >  2:111822/1‑250 (MQ=255)                                         gtgATCAATGGCCCCTGGATGAAAGTTCGCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACaa                                                                                                                                                                                               <  1:340795/250‑1 (MQ=255)                                                              aaaGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTTCCTCTTTGCCTTCCGCTGTTTTAGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCTGCGATATCAGGAATTTCTCGCAGCTGGTGACGCGTTACGGCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAATCCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCgg                                                                                                                                                                         <  1:355409/251‑1 (MQ=255)                                                                 gTTCTCGACTAAGATATCAGCTTCGCGGATCTGCTTTTCCATCACCTCTTTGCCCTCCACTGTTTTGGTATTTAACTCAATAGAACGTTTGATACTGTTCAAAATGGCGCAGTAAAGCGCATCGATCTCAGGAATATCTCGCAGCTCGTGCAGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGAGCCAGATTGCACACCGGTg                                                                                                                                                                       >  1:145387/1‑250 (MQ=255)                                                                       gACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCg                                                                                                                                                                 >  2:277032/1‑250 (MQ=255)                                                                               tatCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGACCttt                                                                                                                                                         >  1:365491/1‑250 (MQ=255)                                                                                          cgGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCACTAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCAACAACGCCAGGACGTTCAATTTTAATAACGTCCGCGCCAAACCAGGCCAGCATATGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGa                                                                                                                                               >  2:38399/1‑249 (MQ=255)                                                                                                   ttttCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTg                                                                                                                                     >  1:164940/1‑250 (MQ=255)                                                                                                           tACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCa                                                                                                                             <  1:31532/250‑1 (MQ=255)                                                                                                             ccTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGGGCAGTAAATCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTACTAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAAt                                                                                                                           >  2:197302/1‑250 (MQ=255)                                                                                                                     ccTTCCGATGTTTTGGTATTTAACTCAATAGAACGTTTGGTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGTCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATTCCTATTAg                                                                                                                  <  2:176023/251‑1 (MQ=255)                                                                                                                         ccGCTGTTTTGGTATTTAACACAATAGATCGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGttaa                                                                                                              <  1:79501/251‑1 (MQ=255)                                                                                                                                 ttGGTATTTAACTCAATAGAACGTTTGTTACTGGTAAGCATGTTGAAGTAAAGCGCATCGAGATCAGGAATATCTCGCAGCTGGTGACGCGGTACGTCTCCAACGCCAGGACGTGCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAGCTGTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAAGTAAata                                                                                                       <  1:321035/250‑1 (MQ=255)                                                                                                                                  tGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCAGCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACTCCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATTCCTATTAGTTAATTAAatat                                                                                                      <  1:124788/250‑1 (MQ=255)                                                                                                                                    gTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATAtgtg                                                                                                   <  1:277032/251‑1 (MQ=255)                                                                                                                                      aTTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGtt                                                                                                 >  1:7639/1‑251 (MQ=255)                                                                                                                                         tAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAAt                                                                                              >  2:207923/1‑251 (MQ=255)                                                                                                                                                      cGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGTATTTGAGAACAAGATGGGCCAGATTGTACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAATTTAGCGtttattt                                                                                 >  1:61393/1‑251 (MQ=255)                                                                                                                                                                tGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGCAGTGGAGtt                                                                                                                                      >  2:82245/1‑188 (MQ=255)                                                                                                                                                                tGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGtt                                                                                                                                      <  1:82245/188‑1 (MQ=255)                                                                                                                                                                                   aaGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAATTTAGCGTTTATTTAAGGCGTAATATCTTCTATAAATAACCCg                                                    <  2:100500/251‑1 (MQ=255)                                                                                                                                                                                               atCAGAAATATCTCGCAGCTGATGACGCGTTAGGACACCAACGCCAGGTCGTTTAATTTTAATAACGTCAGCGCCAAACCAGGATAACATTTGAGTACAAGATGGGCTAGATTGCTCACCTGTGAAATCGCTAACTTTAATTCCTGGAAGTGGAGTGGACATAGCAATACCTATTTGTTAATTAAATATGTGTTAATTTAGCGTTTATTTAAGGCGTAATATCTCCGACTACTAGCCCGCAAAATAtcgtc                                        <  2:86678/251‑1 (MQ=255)                                                                                                                                                                                                   ggAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAATTTAGCGTTTATTTAAGGCGTAATATCTTCTATAAATAACCCGCAAAATATcgtcgt                                      >  2:197409/1‑249 (MQ=255)                                                                                                                                                                                                          cactcAGCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAATTGAGCGTTTAGTTAAGGCGTAATATCTTCTATAAATAACCCGCAAAATATCGTCGTAGCTGCACt                             <  2:30463/247‑1 (MQ=255)                                                                                                                                                                                                                gCTGGTGACGCGTTACGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAATTTAGCGTTTATTTAAGGCGTAATATCTTCTATAAATAACCCGCAAAATATCGTCGTAGCTGCACTCGatat                       <  2:176607/251‑1 (MQ=255)                                                                                                                                                                                                                    gTGACGCGTTACGGCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTAGTAGTTAATTAAATATGTGTTAATTTAGCGTTTATTTAAGGCGTAATATCTTCTATAAATAACCCGCAAAATATCGTCGTAGCTGCACGCGATAttttt                   <  2:342090/251‑1 (MQ=255)                                                                                                                                                                                                                            ttACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTg                                                                                                                                                <  2:165350/118‑1 (MQ=255)                                                                                                                                                                                                                            ttACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTg                                                                                                                                                >  1:165350/1‑118 (MQ=255)                                                                                                                                                                                                                              aCGTCACCAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAATTTAGCGTTTATTTAAGGCGTAATATCTTCTATAAATAACCCGCAAAATATCGTCGTAGCTGCACTCGATATTTTTAAAATCCTg          <  1:111822/250‑1 (MQ=255)                                                                                                                                                                                                                                     cAACGCCAGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAATTTAGCGTTTATTTAAGGCGTAATATCTTCTATAAATAACCCGCAAAATATCGTCGTAGCTGCACTCGATATTTTTAAAATCCTGTATCAAGt  <  2:117739/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    GTGGATTGATTTCTTGAATATGCTCCCAGGTGAAGCCCATGTGATCAATGGCCCCTGGATGAAAGTTCTCGACTAAGATATCAGCTTCGCGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAATACCTATTAGTTAATTAAATATGTGTTAATTTAGCGTTTATTTAAGGCGTAATATCTTCTATAAATAACCCGCAAAATATCGTCGTAGCTGCACTCGATATTTTTAAAATCCTGTATCAAGT  >  NC_000913/2492903‑2493381

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