Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,737,599	G→C	7.4%	intergenic (‑154/+48)	hutU ← / ← PP_5034	urocanate hydratase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,737,599	0	G	C	7.4%	164.7 / 3.5	54	intergenic (‑154/+48)	hutU/PP_5034	urocanate hydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base C (2/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

TGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCC  >  NC_002947/5737459‑5737742                                                                                                                                             |                                                                                                                                                 tGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGggcggc                                                                                                                                          >  6:222060/1‑149 (MQ=255) tGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTAATGTTGGGGTATAGATCGAGCGCCGCGCGggcggc                                                                                                                                          >  7:137999/1‑149 (MQ=255)     gtcgtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                      >  7:83775/1‑149 (MQ=255)     gtcgtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                      <  5:222060/149‑1 (MQ=255)        gtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                    <  1:22007/148‑1 (MQ=255)        gtcCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTAATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                   <  8:137999/149‑1 (MQ=255)        gtcCGCGGGCGGGTATGTCGTGCTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGAAAAACAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCTAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                   <  2:142722/149‑1 (MQ=255)           cgcgGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGAtc‑t                                                                                                                                >  6:261732/1‑148 (MQ=255)                       tGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCaa                                                                                                                    >  6:184996/1‑148 (MQ=255)                           gTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGCGCGGCGCTCGATC‑TCACAGGCGCCAAAAATc                                                                                                               >  5:56011/1‑149 (MQ=255)                            tGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCaaaaa                                                                                                                 <  4:195090/146‑1 (MQ=255)                            tGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCaaaaa                                                                                                                 >  3:195090/1‑146 (MQ=255)                              gTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTa                                                                                                             <  1:139837/148‑1 (MQ=255)                              gTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTa                                                                                                             <  1:144758/148‑1 (MQ=255)                                      aCAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                       >  4:129269/1‑115 (MQ=255)                                           cacacCTTTACTTGTATGTACAAGCATGGACAACCAAAATAAACGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCtt                                                                                                 >  6:150847/1‑147 (MQ=255)                                              acCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAgcg                                                                                            >  7:111038/1‑149 (MQ=255)                                              acCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAgcg                                                                                            <  8:111038/149‑1 (MQ=255)                                                cTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGcc                                                                                          >  4:195993/1‑149 (MQ=255)                                                   tACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGccc                                                                                         <  5:311827/147‑1 (MQ=255)                                                     cTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATc                                                                                      >  2:305745/1‑148 (MQ=255)                                                     cTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATc                                                                                      <  1:305745/148‑1 (MQ=255)                                                     cTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCa                                                                                     <  8:83775/149‑1 (MQ=255)                                                        gtatgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACATGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGt                                                                                  >  5:236916/1‑149 (MQ=255)                                                          atgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACATGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGtt                                                                                 <  6:236916/148‑1 (MQ=255)                                                           tgtaCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTc                                                                                >  3:209732/1‑148 (MQ=255)                                                             taCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGa                                                                              <  3:195993/148‑1 (MQ=255)                                                             taCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                             <  6:56011/149‑1 (MQ=255)                                                              aCAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                             >  2:12875/1‑148 (MQ=255)                                                              aCAAGCATGGACAACCAAAATAATCCGACAATACCGTTCCTGCCTTCAGATTTATGTTGGGTTATAGATCCAGCTCCACGCGCGCGTCGCTCCATA‑TCACATGCGCCAAAAATATACAGCCGAACCTTCAGCGCCCATCAAGTTCGAt                                                                             >  2:25479/1‑148 (MQ=255)                                                               cAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATc                                                                            <  5:261732/148‑1 (MQ=255)                                                               cAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATc                                                                            <  6:305717/148‑1 (MQ=255)                                                                aaGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAg                                                                          >  3:57892/1‑149 (MQ=255)                                                                       gACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTCGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCAc                                                                   >  3:68993/1‑149 (MQ=255)                                                                          aaCCCAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGCGGTATAGATCGAGCGCCCCGCGTGCGGCGCTCGAAC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGAGCACATCAACTTCGAACAGGCAGCACCGc                                                                >  5:314310/1‑149 (MQ=255)                                                                          aaCCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGgcag                                                                      >  1:220914/1‑143 (MQ=255)                                                                          aaCCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGgcag                                                                      <  2:220914/143‑1 (MQ=255)                                                                             cAAAATAATCGTACCTTACCGATACTGGCTTGAGATTTATGTTGGCGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGcccc                                                             >  2:287408/1‑149 (MQ=255)                                                                              aaaaTAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                         <  5:82928/136‑1 (MQ=255)                                                                              aaaaTAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAgg                                                                         >  6:82928/1‑136 (MQ=255)                                                                               aaaTAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTg                                                           <  4:209732/149‑1 (MQ=255)                                                                                 ataatCGGACCTTACCGTTCCTGGCTTGAGATTTATGTCGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGca                                                         <  4:68993/149‑1 (MQ=255)                                                                                        gACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAg                                                  <  5:184996/149‑1 (MQ=255)                                                                                        gACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTACGATCAGCCAGCACAGCCCCTGCACAACCAg                                                  >  6:58536/1‑149 (MQ=255)                                                                                                      tGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTTGACCGGCGCCAAAAATGTACAGCCGAACCTTCAGCGCCCATAAATTTCGAACAGGCAGCCCCGCCCCTGCAAAACCACCCCCAGCAATtca                                     >  1:90564/1‑149 (MQ=255)                                                                                                             agaTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGc                              <  4:57892/148‑1 (MQ=255)                                                                                                               aTTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCtt                           <  1:12875/149‑1 (MQ=255)                                                                                                                 ttATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCATAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTc                          <  5:58536/148‑1 (MQ=255)                                                                                                                                gATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTc           <  1:287408/148‑1 (MQ=255)                                                                                                                                gATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTc           <  5:150847/148‑1 (MQ=255)                                                                                                                                gATCGAGCGCCGCCCGCGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTc           >  6:316485/1‑148 (MQ=255)                                                                                                                                  tCGAGCGCCGCCCGCGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTg                               <  7:101177/126‑1 (MQ=255)                                                                                                                                  tCGAGCGCCGCCCGCGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTg                               >  8:101177/1‑126 (MQ=255)                                                                                                                                         ccgcGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCACCCCCACCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAAccc  >  6:58117/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                                 TGCGGTCGTCCGCGGGCGGGTATGTCGTGGTCAATATACAAACACACCTTTACTTGTATGTACAAGCATGGACAACCAAAATAATCGGACCTTACCGTTCCTGGCTTGAGATTTATGTTGGGGTATAGATCGAGCGCCGCGCGGGCGGCGCTCGATC‑TCACAGGCGCCAAAAATCTACAGCCGAACCTTCAGCGCCCATCAAGTTCGATCAGGCAGCACCGCCCCTGCACAACCAGCCCCAGCAATTCATCATTGCCTTCCAGCCGCAAGCAGTCGAACAACCC  >  NC_002947/5737459‑5737742

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