Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,486,218	G→T	13.7%	V310V (GTG→GTT)	purD →	phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,486,218	0	G	T	13.7%	98.5 / 4.0	35	V310V (GTG→GTT)	purD	phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (18/12);  new base T (3/2);  total (21/14)
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

GCAACGTCTACACCGGCTTCCTCTACGCCGGCCTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCGGGCGACTACGCCAAGGGCGACGTGATCAACGGCCTGGACGCTGCCGCCAAGAT  >  NC_002947/5486073‑5486364                                                                                                                                                  |                                                                                                                                                   gCAACGTCTACACCGGCTTCCTCTACGCCGGCCTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAg                                                                                                                                                >  2:750029/1‑150 (MQ=255)    aCGTCTACACCGGCTTCCTCTACGCCGGCCTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCg                                                                                                                                             >  1:146250/1‑150 (MQ=255)      gTCTACACCGGCTTCCTCTACGCCGGCCTGATGATCGACAAGGCGGGCAACCCCAAGGGCATCGAGTTCAAATATCGCTTCGGCGACCCGTAAACCCAGGCGGTCATGGTGCGCCTGGACTTAGGCTTAGTCTTGCTGGTTGAAAGCGCa                                                                                                                                           >  2:341739/1‑150 (MQ=255)           cacCGGCTTCCTCTACGCCGGCCTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGc                                                                                                                                      <  1:88101/150‑1 (MQ=255)                     ctctACGCCGGCCTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTg                                                                                                                            >  2:758729/1‑150 (MQ=255)                      tctACGCCGGCCTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTgg                                                                                                                           >  1:415456/1‑150 (MQ=255)                      tctACGCCGGCCTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTgg                                                                                                                           >  2:515244/1‑150 (MQ=255)                              ggCCTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTc                                                                                                                   <  2:146250/150‑1 (MQ=255)                                 cTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGaa                                                                                                                <  2:206514/150‑1 (MQ=255)                                  tgatgaTCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAg                                                                                                               <  1:515244/150‑1 (MQ=255)                                        tCGACAGGGCGGGCAGCCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGt                                                                                                         >  2:212626/1‑150 (MQ=255)                                        tCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGt                                                                                                         >  1:798937/1‑150 (MQ=255)                                             aaGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAgg                                                                                                                     >  1:536552/1‑133 (MQ=255)                                             aaGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAgg                                                                                                                     <  2:536552/133‑1 (MQ=255)                                              aGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGAcc                                                                                                   >  2:602857/1‑150 (MQ=255)                                               ggCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGACGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGAccc                                                                                                  >  2:325723/1‑150 (MQ=255)                                                          ccAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTTGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGAccc                                                                                                  <  1:658676/139‑1 (MQ=255)                                                          ccAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTTGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGAccc                                                                                                  >  2:658676/1‑139 (MQ=255)                                                            aaGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTg                                                                                     <  1:325723/150‑1 (MQ=255)                                                                                ccGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAgcggg                                                                 >  1:15750/1‑150 (MQ=255)                                                                                           aCCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCgg                                                      <  1:212626/150‑1 (MQ=255)                                                                                             ccGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAGGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCgcgc                                                                                              >  1:715288/1‑108 (MQ=255)                                                                                             ccGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAGGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCgcgc                                                                                              <  2:715288/108‑1 (MQ=255)                                                                                                            gTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTTGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAg                                                                                                               >  2:354143/1‑76 (MQ=255)                                                                                                            gTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTTGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAg                                                                                                               >  2:354167/1‑76 (MQ=255)                                                                                                            gTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTTGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAg                                                                                                               <  1:354167/76‑1 (MQ=255)                                                                                                            gTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTTGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAg                                                                                                               <  1:354143/76‑1 (MQ=255)                                                                                                                   tgcGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCGGGCGACTATGCCAAGGGCGACGTGa                              >  1:274992/1‑150 (MQ=255)                                                                                                                   tgcGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCGGGCGACTACGCCAAGGGCGACGTGa                              >  2:371580/1‑150 (MQ=255)                                                                                                                         tGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTgg                                                                        >  2:510385/1‑102 (MQ=255)                                                                                                                         tGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTgg                                                                        <  1:510385/102‑1 (MQ=255)                                                                                                                                gAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCGGGCGACTACGCCAAGGGCGACGTGATCAACGGCCTGGa                 >  2:433877/1‑150 (MQ=255)                                                                                                                                    cTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCGGGCGACTACGCCAAGGGCGACGTGATCAACGGCCTGGACGCt             <  2:367701/150‑1 (MQ=255)                                                                                                                                         ccTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCGGGCGACTACGCCAAGGGCGACGTGATCAACGGCCTGGACGCTgccgc        <  1:602857/150‑1 (MQ=255)                                                                                                                                             ctgGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCGGGCGACTACGCCAAGGGCGACGTGATCAACGGCCTGGACGCTGCCGCCAAg    >  1:432683/1‑150 (MQ=255)                                                                                                                                               ggtggAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCGGGCGACTACGCCAAGGGCGACGTGATTAACGGCCTGGACGCTGCCGCCAAGAt  <  2:798937/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GCAACGTCTACACCGGCTTCCTCTACGCCGGCCTGATGATCGACAAGGCGGGCAACCCCAAGGTCATCGAGTTCAACTGCCGCTTCGGCGACCCGGAAACCCAGCCGGTCATGCTGCGCCTGGAGTCGAGCCTGGTCCTGCTGGTGGAAGCCGCATTCGCCAAGGCACTGGACAAGGTCGAAGCACAGTGGGACCCGCGCCCGAGCCTGGGCGTGGTCCTGGCAGCGGGCGGTTACCCGGGCGACTACGCCAAGGGCGACGTGATCAACGGCCTGGACGCTGCCGCCAAGAT  >  NC_002947/5486073‑5486364

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑