Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,622,361	T→C	L43P (CTT→CCT)	purN →	phosphoribosylglycinamide formyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,622,361	0	T	C	98.2%	184.5 / NA	54	L43P (CTT→CCT)	purN	phosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (28/25);  minor base G (1/0);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGG  >  NC_000913/2622227‑2622489                                                                                                                                       |                                                                                                                                 tCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAAc                                                                                                                             >  1:455141/1‑140 (MQ=255)    aataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAAcgcg                                                                                                                          <  2:23512/140‑1 (MQ=255)       aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAAcgcg                                                                                                                          <  2:459729/137‑1 (MQ=255)       aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAAcgcg                                                                                                                          >  1:459729/1‑137 (MQ=255)       aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGcgccc                                                                                                                       <  2:598464/140‑1 (MQ=255)       aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGcgccc                                                                                                                       <  2:704526/140‑1 (MQ=255)        atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGcgcccg                                                                                                                      >  2:177378/1‑140 (MQ=255)           aaTATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCGATAAGGCGGACACGTTCGGCCGTGAACGCGCCCGCCa                                                                                                                   >  1:222518/1‑140 (MQ=255)           aaTATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCa                                                                                                                   >  2:482829/1‑140 (MQ=255)              atTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGc                                                                                                                >  1:31886/1‑140 (MQ=255)                  tggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGt                                                                                                            >  2:121705/1‑140 (MQ=255)                   ggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTa                                                                                                           >  2:411/1‑140 (MQ=255)                     tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTAtt                                                                                                         >  1:203838/1‑140 (MQ=255)                       cTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                       <  1:411/140‑1 (MQ=255)                       cTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                       >  1:495692/1‑140 (MQ=255)                       cTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                       >  1:611360/1‑140 (MQ=255)                          aTTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                    >  1:38614/1‑140 (MQ=255)                                              tttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgcg                                                                                <  2:455141/140‑1 (MQ=255)                                                tACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGcgcg                                                                              >  2:67158/1‑140 (MQ=255)                                                  cAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACTCGCCCGCCCGGCAGGTATTGCAACGTATACCCTCATCGCCAGCGCGtt                                                                            >  2:274751/1‑140 (MQ=255)                                                      cAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCTCGCCAGGCGGGTATTGTAACGCATAAGCTCTTCGCCAGCGCGTTTGAc                                                                        >  2:297175/1‑140 (MQ=255)                                                       aaTTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTggcc                                                                       >  1:103114/1‑137 (MQ=255)                                                            ttGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                  >  1:558614/1‑140 (MQ=255)                                                             tGACGCCTGTAAAACCAGCAAATTTAAAGGCACCGTACGCGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTg                                                                 <  2:222518/140‑1 (MQ=255)                                                             tGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTg                                                                 <  1:177378/140‑1 (MQ=255)                                                             tGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTg                                                                 <  2:31886/140‑1 (MQ=255)                                                               aCGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCCGGTATTGCTACGTATAAGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                               >  2:246329/1‑140 (MQ=255)                                                                 gCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCgcctgac                                                             >  2:434274/1‑134 (MQ=255)                                                                   cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                           >  2:612366/1‑140 (MQ=255)                                                                   cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                           >  2:474678/1‑140 (MQ=255)                                                                   cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                           >  1:492249/1‑140 (MQ=255)                                                                    tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTa                                                          <  2:495692/140‑1 (MQ=255)                                                                    tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTa                                                          <  1:529219/140‑1 (MQ=255)                                                                         aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAcc                                                     <  1:482829/140‑1 (MQ=255)                                                                             aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                 >  1:678948/1‑140 (MQ=255)                                                                              acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAg                                                >  1:547114/1‑140 (MQ=255)                                                                                aaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGtt                                              <  2:38614/140‑1 (MQ=255)                                                                                        aGGCACCGTACGGGCAGTTTTCCGCCATAAGGCCGACACGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATg                                      <  1:274751/140‑1 (MQ=255)                                                                                        aGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATg                                      <  2:547114/140‑1 (MQ=255)                                                                                            aCCGTACGGACAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAAt                                  >  1:658748/1‑140 (MQ=255)                                                                                                tACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGAc                              >  2:762863/1‑140 (MQ=255)                                                                                                tACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGAc                              <  1:67158/140‑1 (MQ=255)                                                                                                     gCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                         <  1:592968/140‑1 (MQ=255)                                                                                                     gCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                         >  2:592968/1‑140 (MQ=255)                                                                                                     gCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCAGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                         <  1:759376/140‑1 (MQ=255)                                                                                                       aGTTTTCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACg                       <  2:658748/140‑1 (MQ=255)                                                                                                         ttttCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCa                     <  2:492249/140‑1 (MQ=255)                                                                                                         ttttCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCa                     <  1:9766/140‑1 (MQ=255)                                                                                                           ttCAGCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCAcc                   <  1:121705/140‑1 (MQ=255)                                                                                                               gCAATAAGGCCGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAt               >  2:91855/1‑140 (MQ=255)                                                                                                                        ccGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTg      <  2:611360/140‑1 (MQ=255)                                                                                                                        ccGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTg      <  2:678948/140‑1 (MQ=255)                                                                                                                        ccGACGCGTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTg      <  1:474678/140‑1 (MQ=255)                                                                                                                            cgcgTTCGGCCCTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg  >  2:482073/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                 TCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGG  >  NC_000913/2622227‑2622489

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