Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,145,387	T→C	7.8%	I62T (ATC→ACC)	atnUDX →	NUDIX hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,145,387	0	T	C	7.8%	146.0 / 2.6	51	I62T (ATC→ACC)	atnUDX	NUDIX hydrolase
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGGAACTGGCCCGGGGCGTTCCCACAGGCGATACCCACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTGCAGACATCGTCTCGCCGTATTCGATCTTCAGCGT  >  NC_002947/4145240‑4145528                                                                                                                                                    |                                                                                                                                               aGGAACTGGCCCGGGGCGTTCCCACAGGCGATACCCACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATc                                                                                                                                               <  1:244114/149‑1 (MQ=255)            cGGGGCGTTCCCACAGGCGATACCCACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGccc                                                                                                                                    <  6:181344/149‑1 (MQ=255)                       cacaGGCGATACCCACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGa                                                                                                                         >  5:196448/1‑149 (MQ=255)                        acaGGCGATACCCACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCACCCCCCCTCGCCCCGGCACCTGGAc                                                                                                                        >  8:70902/1‑149 (MQ=255)                           ggCGATACCCACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCtt                                                                                                                     <  3:308796/149‑1 (MQ=255)                                tACCCACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTtgcctg                                                                                                                >  5:91496/1‑149 (MQ=255)                                   ccACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCg                                                                                                             >  8:224556/1‑149 (MQ=255)                                   ccACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCg                                                                                                             >  4:173649/1‑149 (MQ=255)                                    cACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCACCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCgg                                                                                                            >  5:67271/1‑149 (MQ=255)                                      cGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGt                                                                                                          <  4:45484/149‑1 (MQ=255)                                      cGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGt                                                                                                          <  4:126336/149‑1 (MQ=255)                                           gCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATg                                                                                                     <  1:186281/149‑1 (MQ=255)                                            ccTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATgg                                                                                                    <  3:173649/149‑1 (MQ=255)                                                  ttGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTg                                                                                              <  2:189150/149‑1 (MQ=255)                                                         ggCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAac                                                                                       <  5:75769/149‑1 (MQ=255)                                                         ggCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAac                                                                                       <  5:148361/149‑1 (MQ=255)                                                           cTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAacc‑a                                                                                    <  2:61251/149‑1 (MQ=255)                                                                   aCATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGaaa                                                                                        >  1:142287/1‑138 (MQ=255)                                                                   aCATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGaaa                                                                                        <  2:142287/138‑1 (MQ=255)                                                                     aTCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAggc                                                                          >  8:124584/1‑149 (MQ=255)                                                                      tCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAggcg                                                                         >  3:365005/1‑149 (MQ=255)                                                                      tCCATTACATCAACCCCAAGATCATCGCCAGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAggcg                                                                         >  3:238548/1‑149 (MQ=255)                                                                         aTTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGgcggcg                                                                      >  6:313568/1‑149 (MQ=255)                                                                          ttACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGgcgc                                                                     >  7:329658/1‑149 (MQ=255)                                                                                   ccccAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAgg                                                            <  7:224556/149‑1 (MQ=255)                                                                                   ccccAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAgg                                                            <  5:313568/149‑1 (MQ=255)                                                                                   ccccAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAgg                                                            >  5:155102/1‑149 (MQ=255)                                                                                    cccAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAgg                                                            >  6:202329/1‑148 (MQ=255)                                                                                    cccAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAgg                                                            >  3:185065/1‑148 (MQ=255)                                                                                     ccAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTa                                                          >  2:274340/1‑149 (MQ=255)                                                                                       aaGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAggc                                                                          >  2:273328/1‑131 (MQ=255)                                                                                       aaGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAggc                                                                          <  1:273328/131‑1 (MQ=255)                                                                                           tcatcGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCACCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGgaa                                                    <  7:70902/149‑1 (MQ=255)                                                                                              tcGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGgaaga                                                  >  2:331182/1‑148 (MQ=255)                                                                                              tcGCCGGCTGCATCATCGAGAGCGATGGAAACTACCTGCTGTGCCAGCGCCCCATCCCCCCTCGCCCCGGCACCTGGACCATGCCAGCCGGGTTCATGGCGGCTGGCGAAACA‑ACCGAACAGGCCGCGCAGCGCGTCGTATGGgtaga                                                  >  4:6675/1‑148 (MQ=255)                                                                                               cGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAg                                                             <  5:270211/136‑1 (MQ=255)                                                                                               cGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAg                                                             >  6:270211/1‑136 (MQ=255)                                                                                                 ccGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAgc                                              <  4:365005/149‑1 (MQ=255)                                                                                                     cTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAAc                                                                              >  1:113910/1‑113 (MQ=255)                                                                                                     cTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAAc                                                                              <  2:113910/113‑1 (MQ=255)                                                                                                       gCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGaa                                                                                         >  1:7422/1‑101 (MQ=255)                                                                                                       gCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGaa                                                                                         <  2:7422/101‑1 (MQ=255)                                                                                                        catcatCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCACCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTAc                                       <  6:67271/149‑1 (MQ=255)                                                                                                            atcGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACCTACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTg                                    <  1:274340/149‑1 (MQ=255)                                                                                                            atcGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTg                                    <  6:155102/148‑1 (MQ=255)                                                                                                                         gAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTGCAGACATCGTCTCg                      >  4:215459/1‑149 (MQ=255)                                                                                                                          aaaGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTGCAGACATCGTCTCGc                     >  7:258433/1‑149 (MQ=255)                                                                                                                                ccTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTGCAGACATCGTCTCGCCGTAtt               <  4:185065/149‑1 (MQ=255)                                                                                                                                          cAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTGCAGACATCGTCTCGCCGTATTCGATCTTCAg     >  6:287988/1‑149 (MQ=255)                                                                                                                                          cAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTGCAGACATCGTCTCGCCGTATTCGATCTTCAg     >  4:102218/1‑149 (MQ=255)                                                                                                                                             cgcgcCATCCCCCCTCGCCCCGTCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTGCAGACATCGTCTCGCCGTATTCGATCTTCAGCGt  <  6:196448/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               AGGAACTGGCCCGGGGCGTTCCCACAGGCGATACCCACGAACGCCTGCATTGTACCGGCTGCGGCTACATCCATTACATCAACCCCAAGATCATCGCCGGCTGCATCATCGAGCGCGATGGAAAGTACCTGCTGTGCCAGCGCGCCATCCCCCCTCGCCCCGGCACCTGGACCTTGCCTGCCGGGTTCATGGAGGCTGGCGAAACC‑ACCGAACAGGCGGCGCTGCGCGAGGTATGGGAAGAAAGCGGCGTACGTGCAGACATCGTCTCGCCGTATTCGATCTTCAGCGT  >  NC_002947/4145240‑4145528

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