Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,887,078	Δ1 bp	16.7%	coding (1167/1356 nt)	uacT →	uric acid permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,887,078	0	G	.	16.7%	108.9 / 16.4	36	coding (1167/1356 nt)	uacT	uric acid permease
Reads supporting (aligned to +/- strand):  ref base G (10/20);  new base . (3/3);  total (13/23)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCGCCCGCAGTACTGGGCGGCGCGTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTC  >  NC_002947/4886957‑4887210                                                                                                                          |                                                                                                                                     tCCCGCCCGCAGTACTGGGCGGCGCGTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCgg                                                                                                                        <  1:126039/136‑1 (MQ=255) tCCCGCCCGCAGTACTGGGCGGCGCGTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCgg                                                                                                                        <  1:51167/136‑1 (MQ=255)   ccgcccgcAGTACTGGGCGGCGCGGCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGt                                                                                                                      <  2:279141/136‑1 (MQ=255)             ttcttggcggcTCGTCCATCTCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGTCGGACATCTGCGACCGCAGCATCCTGCTGTTGGTTGCCGTGAGCTTCTGCTTGTGGTTGATCCACgt                                                                                                            <  2:40359/131‑1 (MQ=255)                  ggcggcGCGTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGCCATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGT‑AGCGTCGGCTTCGGGGTGATGCCCCTGGTGCg                                                                                                      >  1:317825/1‑136 (MQ=255)                  ggcggcGCGTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGT‑AGCGTCGGCTTCGGGTTGATCCCCGTGGTGCg                                                                                                      >  1:171591/1‑136 (MQ=255)                  ggcggcGCGTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGT‑AGCGTCGGCTTCGGGTTGATCCCCGTGGTGCg                                                                                                      >  1:433824/1‑136 (MQ=255)                     ggcgCGTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTTATACCCCTTGTGCGTc                                                                                                    >  2:59653/1‑136 (MQ=255)                      gcgcGTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTcc                                                                                                   <  2:431986/136‑1 (MQ=255)                          gTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAg                                                                                               >  2:427421/1‑136 (MQ=255)                          gTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAg                                                                                               >  2:341793/1‑136 (MQ=255)                          gTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAg                                                                                               >  2:331130/1‑136 (MQ=255)                          gTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAg                                                                                               >  2:27222/1‑136 (MQ=255)                          gTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAg                                                                                               >  1:165128/1‑136 (MQ=255)                          gTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAg                                                                                               >  1:52774/1‑136 (MQ=255)                                 gCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTg                                                                                        <  2:197002/136‑1 (MQ=255)                                 gCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTg                                                                                        <  1:210745/136‑1 (MQ=255)                                    aTGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGcca                                                                                     <  2:34707/136‑2 (MQ=255)                                           gCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGcc                                                                              <  1:237352/136‑1 (MQ=255)                                           gCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGcc                                                                              <  1:356834/136‑1 (MQ=255)                                               ggTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAg                                                                          >  2:377383/1‑136 (MQ=255)                                                         aCCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGc                                                                <  1:27222/136‑1 (MQ=255)                                                                tCAAGATTCTTCAGGAGGCGGAGATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTTATCCCCGTGGTGCGTCCGTAGTTCTTTGCCGAGATGCCACAGTGGATGGAGCCCATCAc                                                         <  1:96229/136‑1 (MQ=255)                                                                                   ggACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATg                                      <  1:469099/136‑1 (MQ=255)                                                                                                cgccgcAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCtgtg                         <  1:377383/136‑1 (MQ=255)                                                                                                 gccgcAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTcc                                                                                                   >  1:34373/1‑61 (MQ=255)                                                                                                 gccgcAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTcc                                                                                                   <  2:34373/61‑1 (MQ=255)                                                                                                           gCTGTTGGTTGCCGT‑AGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGa             <  2:317825/136‑1 (MQ=255)                                                                                                           gCTGTTGGTTGCCGT‑AGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGa             <  2:171591/136‑1 (MQ=255)                                                                                                           gCTGTTGGTTGCCGT‑AGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGa             <  2:433824/136‑1 (MQ=255)                                                                                                              gttggttgCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCTGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAAc           <  1:341793/136‑1 (MQ=255)                                                                                                              gttggttgCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAAc           <  1:163117/136‑1 (MQ=255)                                                                                                                      ccGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGtt   <  2:478769/136‑1 (MQ=255)                                                                                                                       cGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTc  <  2:52774/136‑1 (MQ=255)                                                                                                                         tGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGt                                                                         >  2:165698/1‑63 (MQ=255)                                                                                                                         tGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGt                                                                         <  1:165698/63‑1 (MQ=255)                                                                                                                          |                                                                                                                                     TCCCGCCCGCAGTACTGGGCGGCGCGTCCATCGCCATGTTCGGCATGGTCACCGCCACCGGGATCAAGATTCTTCAGGAGGCGGACATCGGCGACCGCCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTC  >  NC_002947/4886957‑4887210

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