Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,528,041	A→G	8.7%	K363E (AAA→GAA)	purB →	adenylosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,528,041	0	A	G	8.7%	124.0 / 3.2	46	K363E (AAA→GAA)	purB	adenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base G (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCTCAAAGAGCTGACCCGTG  >  NC_002947/4527896‑4528188                                                                                                                                                  |                                                                                                                                                     ccTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAg                                                                                                                                                   <  4:166732/149‑1 (MQ=255)          cGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTTCAGGCTGGCCTCGAAGGcgatg                                                                                                                                             >  8:60713/1‑142 (MQ=255)          cGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTTCAGGCTGGCCTCGAAGGcgatg                                                                                                                                             <  7:60713/146‑5 (MQ=255)           gCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGc                                                                                                                                           >  3:273790/1‑147 (MQ=255)                  cgcTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGaa                                                                                                                                  <  3:110008/149‑1 (MQ=255)                    cTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCAAACCGTAATCGCCTTCGACGCCAGCCTCAAACGAATCCGCAAACTGGaat                                                                                                                                 >  8:241333/1‑147 (MQ=255)                       ttCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTTCAGGCTGGCCTCGAAGGcgatga                                                                                                                                            <  4:295685/134‑6 (MQ=255)                       ttCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTTCAGGCTGGCCTCGAAGGcgatga                                                                                                                                            >  3:295685/1‑129 (MQ=255)                       ttCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCaa                                                                                                                             <  3:18885/149‑1 (MQ=255)                          cAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCAt                                                                                                                                               >  1:207938/1‑128 (MQ=255)                          cAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCAt                                                                                                                                               <  2:207938/128‑1 (MQ=255)                               ccTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTc                                                                                                                     <  2:122116/149‑1 (MQ=255)                                              gCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGAcc                                                                                                      >  2:219780/1‑149 (MQ=255)                                                       gcgcTTGCATCGTGACCTGACCGACTACACCGTGCTGCGCAAAATTGGCGTGGGCTTTGCCCAAAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGa                                                                                                                                   <  1:315522/111‑1 (MQ=255)                                                       gcgcTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTTGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTgg                                                                                                   <  4:226116/143‑1 (MQ=255)                                                        cgcTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCct‑g                                                                                           >  5:226084/1‑149 (MQ=255)                                                         gcTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCct‑gc                                                                                          <  4:273790/149‑1 (MQ=255)                                                           tGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGACTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCct‑gctg                                                                                        >  2:234477/1‑149 (MQ=255)                                                                 cGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCct‑gctg                                                                                        <  7:303135/143‑1 (MQ=255)                                                                 cGTGACATGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACCGCGTAATCGCCTTCGAGGCCAGCCTGAAACGGCTCGGCAAGCTGGAAGTCAACGAAGCTCGCCTCTACGCCGCCCTGGCCGCct‑ggag                                                                                        >  8:303135/1‑140 (MQ=255)                                                                  gTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGt                                                                                  <  6:226084/148‑1 (MQ=255)                                                                   tgacctgaccGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAACGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGc                                                                                >  4:239368/1‑149 (MQ=255)                                                                      cctgaccGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTgg                                                                             <  6:269506/149‑1 (MQ=255)                                                                            cGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGgccgagc                                                                       <  3:239368/149‑1 (MQ=255)                                                                            cGACTCCACCGTGCTGCGCAACCTGGGCGTGGACTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGgccgagc                                                                       <  1:234477/149‑1 (MQ=255)                                                                                       tgctgcGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGAcc                                                            >  7:41203/1‑149 (MQ=255)                                                                                             gcAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTCGCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGAt                                                       <  1:219780/149‑1 (MQ=255)                                                                                                     ggCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGcgccgc                                                <  8:41203/147‑1 (MQ=255)                                                                                                           ggCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACa                                         <  6:214836/148‑1 (MQ=255)                                                                                                                    cacaGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAAccc                               >  5:71433/1‑149 (MQ=255)                                                                                                                     acaGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGaaaa                                  >  2:52005/1‑145 (MQ=255)                                                                                                                     acaGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGaaaa                                  <  1:52005/145‑1 (MQ=255)                                                                                                                     acaGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAAcccc                              >  5:39882/1‑149 (MQ=255)                                                                                                                             aTCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACgag                        >  8:267931/1‑147 (MQ=255)                                                                                                                                gCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCTGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCt                   >  3:29375/1‑149 (MQ=255)                                                                                                                                  cTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATc                                       >  5:60996/1‑127 (MQ=255)                                                                                                                                  cTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATc                                       <  6:60996/127‑1 (MQ=255)                                                                                                                                   ttCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACgaga                       >  7:5534/1‑142 (MQ=255)                                                                                                                                   ttCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACgaga                       <  8:5534/142‑1 (MQ=255)                                                                                                                                      gAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAATTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCTCAAag              <  7:267931/148‑1 (MQ=255)                                                                                                                                           cAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCTCAAAGAGCTGa        >  8:147240/1‑149 (MQ=255)                                                                                                                                              ccTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCTCAAAGAGCTGAccc     >  1:175837/1‑149 (MQ=255)                                                                                                                                              ccTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCTCAAAGAGCTGAccc     >  8:269805/1‑149 (MQ=255)                                                                                                                                                 gAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCTc                  >  8:151477/1‑133 (MQ=255)                                                                                                                                                 gAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCTc                  <  7:151477/133‑1 (MQ=255)                                                                                                                                                 gAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCTCAAAGAGCTGACCCGTg  >  5:169728/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCTCAAAGAGCTGACCCGTG  >  NC_002947/4527896‑4528188

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