Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,692,772	T→C	E924E (GAA→GAG)	purL ←	phosphoribosylformylglycinamide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,692,772	0	T	C	97.8%	149.5 / ‑4.6	45	E924E (GAA→GAG)	purL	phosphoribosylformylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (21/23);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGTTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGGCCGCCATCAGAGCGGTCGTGATACGCCAGC  >  NC_000913/2692635‑2692907                                                                                                                                          |                                                                                                                                        ccGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGctc                                                                                                                                       >  2:428007/1‑140 (MQ=255)      ggCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCGTGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGt                                                                                                                                  >  1:147601/1‑140 (MQ=255)      ggCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGt                                                                                                                                  >  1:808403/1‑140 (MQ=255)               cacaAAACGGTCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAAcg                                                                                                                         <  2:613984/140‑1 (MQ=255)                        gTCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTATGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGc                                                                                                                <  2:608857/140‑1 (MQ=255)                         tCACCGGAAACCGCCTGCCCTTCATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGc                                                                                                                <  1:586021/139‑1 (MQ=255)                         tCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGc                                                                                                                >  2:723030/1‑139 (MQ=255)                         tCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGc                                                                                                                <  1:723030/139‑1 (MQ=255)                         tCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGc                                                                                                                >  2:586021/1‑139 (MQ=255)                             cGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCg                                                                                                           <  1:734357/140‑1 (MQ=255)                             cGGAAACCACCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCg                                                                                                           <  2:195943/140‑1 (MQ=255)                               gAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTc                                                                                                         >  2:506024/1‑140 (MQ=255)                                    cGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCa                                                                                                    <  2:135785/140‑1 (MQ=255)                                           ccTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGc                                                                                             <  1:501343/140‑1 (MQ=255)                                           ccTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGc                                                                                             <  1:506024/140‑1 (MQ=255)                                                    tGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCcgc                                                                                    >  2:677188/1‑140 (MQ=255)                                                    tGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCcgc                                                                                    >  2:686837/1‑140 (MQ=255)                                                    tGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCcgc                                                                                    >  1:66002/1‑140 (MQ=255)                                                    tGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCcgc                                                                                    >  1:267858/1‑140 (MQ=255)                                                    tGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCcgc                                                                                    >  1:779064/1‑140 (MQ=255)                                                        cacaATCAGCAAGCCCGTGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCa                                                                                <  2:147601/140‑1 (MQ=255)                                                         acaATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCaa                                                                               >  2:571971/1‑140 (MQ=255)                                                         acaATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCaa                                                                               >  1:695166/1‑140 (MQ=255)                                                          caATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAAt                                                                              <  1:677188/140‑1 (MQ=255)                                                                gCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCcaca                                                                        <  1:428007/140‑1 (MQ=255)                                                                     cccATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGAc                                                                   <  2:695166/140‑1 (MQ=255)                                                                      ccATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGAcc                                                                  >  2:190860/1‑140 (MQ=255)                                                                       cATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCa                                                                 >  2:448310/1‑140 (MQ=255)                                                                               gCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGc                                                         >  2:124898/1‑140 (MQ=255)                                                                                  aGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCAt                                                      <  1:686837/140‑1 (MQ=255)                                                                                    tACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGTTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCAttt                                                    >  2:829540/1‑140 (MQ=255)                                                                                    tACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCAttt                                                    >  1:457046/1‑140 (MQ=255)                                                                                      cGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTcc                                                  <  2:779064/140‑1 (MQ=255)                                                                                      cGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTcc                                                  >  1:35820/1‑140 (MQ=255)                                                                                        gACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGc                                                <  2:457046/140‑1 (MQ=255)                                                                                        gACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGc                                                >  1:454919/1‑140 (MQ=255)                                                                                                     cACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACcagc                                   >  1:32392/1‑140 (MQ=255)                                                                                                         gTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGgc                               >  2:590078/1‑140 (MQ=255)                                                                                                         gTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGgc                               <  1:571971/140‑1 (MQ=255)                                                                                                         gTCAGCGGCACGAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGgc                               <  2:32392/140‑1 (MQ=255)                                                                                                                    gAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGGCCGCCATCagag                    <  1:491518/140‑1 (MQ=255)                                                                                                                    gAACCTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGGCCGCCATCagag                    <  1:124898/140‑1 (MQ=255)                                                                                                                        cTGAATCACCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGGCCGCCATCAGAGCGGt                <  1:590078/140‑1 (MQ=255)                                                                                                                               aCCGCACCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGGCCGCCATCAGAGCGGTCGTGATa         >  2:706998/1‑140 (MQ=255)                                                                                                                                    aCCCAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGGCCGCCATCAGAGCGGTCGTGATACGCca    <  2:808403/140‑1 (MQ=255)                                                                                                                                      ccAGCTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGGCCGCCATCAGAGCGGTCGTGATACGCcagc  <  1:743513/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        CCGTTGGCGGTAATCACAAAACGGTCACCGGAAACCGCCTGCCCTACATAATGGACACAATCAGCAAGCCCATGCTGTGCCAGTACGGACTCGACCGCTTCACGGTCAGCGGCACGAACCTGAATCACCGCACCCAGTTCTTCGTTAAACAACGCCGCCAGGCGATCGTCACCCAGAGTGGCGATATCCGCGTCAATGCCACAATGACCAGCAAAGGCCATTTCCGCCAGCGTTACCAGCAGGCCGCCATCAGAGCGGTCGTGATACGCCAGC  >  NC_000913/2692635‑2692907

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