Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,886,835	C→A	11.4%	F308L (TTC→TTA)	uacT →	uric acid permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,886,835	0	C	A	11.4%	116.1 / 2.2	35	F308L (TTC→TTA)	uacT	uric acid permease
Reads supporting (aligned to +/- strand):  ref base C (15/16);  new base A (2/2);  total (17/18)
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

CTGGCCCTGGGCAAGGTGACTGATCGTGAAGTAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGATCCCGC  >  NC_002947/4886713‑4886962                                                                                                                           |                                                                                                                                cTGGCCCTGGGCAAGGTGACTGATCGTGAAGTAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCg                                                                                                                    <  1:148690/136‑1 (MQ=255)           gCAAGGTGACTGATCGTGAAGTAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCa                                                                                                          >  1:447501/1‑136 (MQ=255)                gTGACTGATCGTGAAGTAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGa                                                                                                     <  1:179158/136‑1 (MQ=255)                       aTCGTGAAGTAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGt                                                                                              <  2:331234/136‑1 (MQ=255)                         cGTGAAGTAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTAGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTgc                                                                                            >  1:244356/1‑136 (MQ=255)                         cGTGAAGTAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTAGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTgc                                                                                            >  1:404591/1‑136 (MQ=255)                                tAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCg                                                                                     <  2:447501/136‑1 (MQ=255)                                tAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCg                                                                                     <  2:429707/136‑1 (MQ=255)                                tAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCg                                                                                     >  2:126039/1‑136 (MQ=255)                                    cccTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTAt                                                                                 >  1:270887/1‑136 (MQ=255)                                     ccTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATg                                                                                >  1:52134/1‑136 (MQ=255)                                     ccTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATg                                                                                >  2:357791/1‑136 (MQ=255)                                           aTGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTAATTGCGGGCTTCTTCAAAACCTTCACACACTCCTCGTTCGCCCAGAAAATCGGCATTGTGCAGATGACCGGGGTGAGCTGCCGATATGTCACCa                                                                          >  1:40359/1‑136 (MQ=255)                                                 cGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTgg                                                                    >  2:237352/1‑136 (MQ=255)                                                 cGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTgg                                                                    >  2:356834/1‑136 (MQ=255)                                                    cgcgGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCg                                                                 <  1:270968/136‑1 (MQ=255)                                                     gcgGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGTCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCgg                                                                <  2:399198/136‑1 (MQ=255)                                                          ctgctgTGCGATGCCGGTGCGTCGTTCATTGCGGTCTTCTTCAACACCTTCACCCACTCCTCTTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCAc                                                           <  2:277786/136‑1 (MQ=255)                                                                             cgtcgtTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCGCTGCTGATCCTGCTCAGCct                                        >  1:197002/1‑136 (MQ=255)                                                                             cgtcgtTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCct                                        >  2:451927/1‑136 (MQ=255)                                                                             cgtcgtTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCct                                        >  2:321359/1‑136 (MQ=255)                                                                                        gCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAgg                             <  1:451927/136‑1 (MQ=255)                                                                                        gCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAgg                             <  1:321359/136‑1 (MQ=255)                                                                                         cGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAggc                            <  1:357791/136‑1 (MQ=255)                                                                                          gggCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAggcg                           >  2:255344/1‑136 (MQ=255)                                                                                            gCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAggcggc                         >  1:145237/1‑136 (MQ=255)                                                                                                              cccACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGAt       <  2:270887/136‑1 (MQ=255)                                                                                                              cccACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGAt       <  1:255344/136‑1 (MQ=255)                                                                                                              cccACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGAt       >  1:459196/1‑136 (MQ=255)                                                                                                              cccACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGAt       >  1:34707/1‑136 (MQ=255)                                                                                                              cccACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGAt       <  2:459196/136‑1 (MQ=255)                                                                                                                cACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGATcc     <  1:26460/136‑1 (MQ=255)                                                                                                                   tcctcGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGATcccgc  <  2:52134/136‑1 (MQ=255)                                                                                                                   tcctcGTTAGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGATcccgc  <  2:244356/136‑1 (MQ=255)                                                                                                                   tcctcGTTAGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGATcccgc  <  2:404591/136‑1 (MQ=255)                                                                                                                           |                                                                                                                                CTGGCCCTGGGCAAGGTGACTGATCGTGAAGTAACCCCTGGCATGCTGCGTCGCGGCCTGCTGTGCGATGCCGGTGCGTCGTTCATTGCGGGCTTCTTCAACACCTTCACCCACTCCTCGTTCGCCCAGAACATCGGCCTGGTGCAGATGACCGGGGTGCGCTGCCGCTATGTCACCATCGTGGCCGGTGCACTGCTGATCCTGCTCAGCCTGCTGCCCAAGGCGGCCTTCCTGATTGCCTCGATCCCGC  >  NC_002947/4886713‑4886962

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