Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,528,028	C→T	9.3%	F358F (TTC→TTT)	purB →	adenylosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,528,028	0	C	T	9.3%	118.9 / 2.7	43	F358F (TTC→TTT)	purB	adenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base C (17/22);  new base T (2/2);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

AAAACTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCT  >  NC_002947/4527880‑4528171                                                                                                                                                     |                                                                                                                                                 aaaaCTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTc                                                                                                                                                  >  5:269506/1‑149 (MQ=255)        cGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCTATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGc                                                                                                                                             <  7:250436/147‑1 (MQ=255)        cGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCa                                                                                                                                           <  4:283427/149‑1 (MQ=255)                 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                                                                                                                           >  3:295685/1‑129 (MQ=255)                                       ttCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTTCAGGCTGGCCTCGAAGGcgatga                                                                                                                           <  4:295685/134‑6 (MQ=255)                                       ttCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCaa                                                                                                            <  3:18885/149‑1 (MQ=255)                                          cAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCAt                                                                                                                              <  2:207938/128‑1 (MQ=255)                                          cAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCAt                                                                                                                              >  1:207938/1‑128 (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)                                                                       gcgcTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGa                                                                                                                                                >  8:183847/1‑81 (MQ=255)                                                                       gcgcTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGa                                                                                                                                                <  7:183847/81‑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                      <  6:60996/127‑1 (MQ=255)                                                                                                                                                  cTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATc                      >  5:60996/1‑127 (MQ=255)                                                                                                                                                   ttCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACgaga      >  7:5534/1‑142 (MQ=255)                                                                                                                                                   ttCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACgaga      <  8:5534/142‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AAAACTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCT‑GCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAAAACCCCTACGAGAAGCT  >  NC_002947/4527880‑4528171

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