Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,861,943	C→T	7.4%	A18T (GCC→ACC)	purM ←	phosphoribosylformylglycinamidine cyclo‑ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,861,943	0	C	T	7.4%	157.7 / 2.7	54	A18T (GCC→ACC)	purM	phosphoribosylformylglycinamidine cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base C (24/26);  new base T (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CCGTCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTGG  >  NC_002947/1861798‑1862085                                                                                                                                                    |                                                                                                                                                ccGTCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTc                                                                                                                                              <  4:69856/149‑1 (MQ=255)   gTCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGa                                                                                                                                            >  5:269441/1‑149 (MQ=255)    tCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGAt                                                                                                                                           <  3:176412/149‑1 (MQ=255)      ggTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGt                                                                                                                                         >  4:74248/1‑149 (MQ=255)          ccGGAGACCAACACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTa                                                                                                                                      <  5:50636/148‑1 (MQ=255)             gagaCCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGAt                                                                                                                                           <  7:128644/140‑1 (MQ=255)             gagaCCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGAt                                                                                                                                           >  8:128644/1‑140 (MQ=255)                  cAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTacaagatag                                                                                                                              >  3:182441/1‑142 (MQ=255)                  cAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTac                                                                                                                                     <  4:182441/141‑1 (MQ=255)                  cAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGt                                                                                                                             >  4:99051/1‑149 (MQ=255)                  cAGCACCGGCTGCTTGTAGCCGGCAGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGAGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCAGGCGTCGATCTCTACACCGGCga                                                                                                                             >  7:255499/1‑148 (MQ=255)                    gCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTcc                                                                                                                           <  6:269441/149‑1 (MQ=255)                    gCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATCACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTCCGCCGGCGTCGATGTCTACACCGGCGTcc                                                                                                                           <  1:24326/149‑1 (MQ=255)                       ccGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATTACTTCAGGGCGTGCGGTGCGTT‑TCAGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCtt                                                                                                                         <  1:96389/149‑1 (MQ=255)                       ccGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTg                                                                                                                        <  3:299635/149‑1 (MQ=255)                        cGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGt                                                                                                                       <  6:191315/149‑1 (MQ=255)                         ggCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTa                                                                                                                      <  2:264732/149‑1 (MQ=255)                              ctTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCATT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCa                                                                                                                 <  2:60201/149‑1 (MQ=255)                                          cggGATCTCGCAGAGGGCGCCGAATCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgct                                                                                                     >  7:47069/1‑149 (MQ=255)                                             gATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgcttgc                                                                                                  <  3:74248/149‑1 (MQ=255)                                             gATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGACGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgcttgc                                                                                                  >  1:75510/1‑149 (MQ=255)                                                ctcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCg                                                                                                                              <  3:240474/118‑1 (MQ=255)                                                ctcGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCg                                                                                                                              >  4:240474/1‑118 (MQ=255)                                                   gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCAtt                                                                                             <  4:183466/148‑1 (MQ=255)                                                   gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCAt                                                                                              <  8:47069/147‑1 (MQ=255)                                                   gCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTg                                                                                            <  8:173320/149‑1 (MQ=255)                                                          gcgcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGagctaca                                                                                     <  3:99051/149‑7 (MQ=255)                                                            gcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGTGTAGACATCGACGCCgg                                                                                                                                >  5:150227/1‑104 (MQ=255)                                                            gcCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGTGTAGACATCGACGCCgg                                                                                                                                <  6:150227/104‑1 (MQ=255)                                                                  gccgccCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTag                                                                              >  1:30421/1‑148 (MQ=255)                                                                     gccCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGA‑CAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAggaggg                                                                         >  4:157206/1‑149 (MQ=255)                                                                                ccATGACTTCAGGGCGTGCGGTGCGTTGT‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgc                                                                                                      >  1:177875/1‑111 (MQ=255)                                                                                ccATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGCCGAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCt                                                                >  3:247337/1‑148 (MQ=255)                                                                                ccATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGttgc                                                                                                      <  2:177875/110‑1 (MQ=255)                                                                                ccATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCt                                                                >  8:66294/1‑148 (MQ=255)                                                                                       ttCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAAc                                                         >  3:91978/1‑148 (MQ=255)                                                                                          aGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGATTCGCCGGCGTCGATGTCTACACCGGCGTCATTGTAGCTCAGGGAGGGTTGCTTGCTAATTGATCCAGGCCTGTAGGAGGGAGGGATTCTTAAAAACGAc                                                      >  5:175876/1‑148 (MQ=255)                                                                                             gCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCAt                                                   >  3:309876/1‑148 (MQ=255)                                                                                                                  aCGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGa                                                       >  5:146766/1‑125 (MQ=255)                                                                                                                  aCGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGa                                                       <  6:146766/125‑1 (MQ=255)                                                                                                                    gCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATaa                             <  7:66294/149‑1 (MQ=255)                                                                                                                               gTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGt                  >  6:153309/1‑149 (MQ=255)                                                                                                                                   gACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTg              >  8:186852/1‑148 (MQ=255)                                                                                                                                    aCCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGaa            <  4:309876/149‑1 (MQ=255)                                                                                                                                     ccAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGaaa           <  5:153309/149‑1 (MQ=255)                                                                                                                                     ccAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGaaa           <  2:75510/149‑1 (MQ=255)                                                                                                                                      cAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAAc          >  8:178033/1‑149 (MQ=255)                                                                                                                                        gTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGt        >  3:195564/1‑149 (MQ=255)                                                                                                                                         aaggacGCCGGTGTAGACATCGACGCCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCt                                  >  7:296156/6‑123 (MQ=38)                                                                                                                                         aaggacGCCGGTGTAGACATCGACGCCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCt                                  <  8:296156/118‑1 (MQ=38)                                                                                                                                           cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGGTTGGTCATTGATCCAGGGCGTTAGGAGGGAGGGATTCTTAAAATCGTCCATGACGGCCAAGGCCGGCTGATAAGCGGCGGGGGTGCTGaagggt        >  1:272816/1‑143 (MQ=255)                                                                                                                                           cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAg     <  7:186852/149‑1 (MQ=255)                                                                                                                                           cTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAg     <  4:91978/149‑1 (MQ=255)                                                                                                                                              cGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTgg  >  7:247000/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                CCGTCGGTGCCGGAGACCAGCACCGGCTGCTTGTAGCCGGCCGGGATCTCGCAGAGGGCGCCGAAGCCGCCCAGGCCACCCATGACTTCAGGGCGTGCGGTGCGTT‑T‑TGCCACGCCCTTGATGCGTTCGACCAGTGCTTCGCCGGCGTCGATGTCTACACCGGCGTCCTTGTAGCTCAGGGAGGGTTGCTTGCTCATTGATCCAGGCCTTTAGGAGGGAGGGATTCTTAAAAACGACCATGACGGCCAAGGCCGGCTGATAAGCGGCGGCGGT‑CTGAAACGTCAGTGG  >  NC_002947/1861798‑1862085

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