Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,660,083	C→A	9.1%	P65T (CCG→ACG)	nudL →	NUDIX hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,660,083	0	C	A	9.1%	127.2 / 4.3	44	P65T (CCG→ACG)	nudL	NUDIX hydrolase
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base A (2/2);  total (23/21)
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.48e-01

ATTGGAAACCGACAGGCGGTTCCCCGAAGCGGCGGTCCTCTTGCCCATTACCCGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCGCTACATGGCCTGAAAGT  >  NC_002947/1659938‑1660207                                                                                                                                                   |                                                                                                                             aTTGGAAACCGACAGGCGGTTCCCCGAAGCGGCGGTCCTCTTGCCCATTACCCGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCg                                                                                                                            >  7:264351/1‑148 (MQ=255) aTTGGAAACCGACAGGCGGTTCCCCGAAGCGGCGGTCCTCTTGCCCATTACCCGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCg                                                                                                                            >  4:345010/1‑148 (MQ=255)                gCGGTTCCCCGAAGCGGCGGTCCTCTTGCCCATTACCCGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTgg                                                                                                            >  1:287143/1‑149 (MQ=255)                   gTTCCCCGAAGCGGCGGTCCTCTTGCCCATTACCCGCAGCGAAGCGCACGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCt                                                                                                         >  3:52673/1‑149 (MQ=255)                    ttCCCCGAAGCGGCGGTCCTCTTGCCCATTACCCGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCtt                                                                                                        <  7:197077/149‑1 (MQ=255)                      ccccGAAGCGGCGGTCCTCTTGCCCATTACCCGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCttt                                                                                                       <  1:185805/148‑1 (MQ=255)                                             ccATTACCCGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTa                                                                                                      <  4:181494/126‑1 (MQ=255)                                             ccATTACCCGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTa                                                                                                      >  3:181494/1‑126 (MQ=255)                                                    ccGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCg                                                                        >  8:62043/1‑149 (MQ=255)                                                    ccGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCg                                                                        >  6:304199/1‑149 (MQ=255)                                                    ccGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCg                                                                        <  3:345010/149‑1 (MQ=255)                                                    ccGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCg                                                                        >  2:302606/1‑149 (MQ=255)                                                     cGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCgg                                                                       <  5:304199/149‑1 (MQ=255)                                                                gcCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCcctg                                                            >  2:15961/1‑149 (MQ=255)                                                                 cccGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCCGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCACcctgg                                                           >  7:57976/1‑149 (MQ=255)                                                                     aaCTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCcctggcctg                                                       <  7:62043/149‑1 (MQ=255)                                                                        tGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGACCCCTGGCCtggtg                                                    >  1:143455/1‑149 (MQ=255)                                                                          gTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCt                                                                    <  2:73870/131‑1 (MQ=255)                                                                          gTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCt                                                                    >  1:73870/1‑131 (MQ=255)                                                                          gTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCtggtg                                                    >  8:442/1‑147 (MQ=255)                                                                          gTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCtggtg                                                    <  2:287143/147‑1 (MQ=255)                                                                          gTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACACGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCtggtg                                                    >  6:110801/1‑147 (MQ=255)                                                                          gTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACACGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCtggtg                                                    >  6:116017/1‑147 (MQ=255)                                                                            cTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAgg                                                >  5:297166/1‑149 (MQ=255)                                                                            cTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAgg                                                >  7:268509/1‑149 (MQ=255)                                                                                  ccTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTg                                          <  6:297166/149‑1 (MQ=255)                                                                                  catGCTCGCCAAGGGCCTGGCTACCCACGGTGGCGAAGTTGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTTGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGAtt                                           <  8:57976/146‑1 (MQ=255)                                                                                         gcCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCtggt                                                     >  8:298/1‑131 (MQ=255)                                                                                           cAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTTCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCa                                 <  8:264351/149‑1 (MQ=255)                                                                                            aaGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAg                                <  7:442/149‑1 (MQ=255)                                                                                                    gTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGCAGCCCCGGAAGACCCGGACCTGGACTTTACTCCCCTGCGCGAGCCTGAACAAGAAATCGCCATGCCCCCTGGCAGGGAGGAGGTGATTGGTCCGCTCCGCCCGCt                          >  1:250312/1‑148 (MQ=255)                                                                                                         ccccCGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGAGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCg                   <  2:250312/149‑1 (MQ=255)                                                                                                             cGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCGCTAc               <  1:302606/149‑1 (MQ=255)                                                                                                                    gAAGTAGCGTTTCCCGGTGGCCGCCGC‑GACCCGGAACACCCGGACCTGTTCTTTACTGCCCTGCGCGAGGGTGATGAAGAAATCGGCCTGCCCCCTGGCCTTGTGGATGTGATTGGTCCGCTCATCGCGCTGATTTCGCTACATGGCCt        <  4:52673/149‑1 (MQ=255)                                                                                                                    gAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGAt                                            <  7:311437/113‑1 (MQ=255)                                                                                                                    gAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGAt                                            >  8:311437/1‑113 (MQ=255)                                                                                                                       gTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCGCTACATGGCCTGa      >  2:232673/1‑148 (MQ=255)                                                                                                                        tAGCCTTTCCCGGTGGCCGCCGC‑GACCGGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGGTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCGCTACATGGCCTGaaa    <  6:264857/149‑1 (MQ=255)                                                                                                                          gCCTTTCCCGGTGGCCGCCGC‑GACACGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCGCTACATGGCCTGAAAGt  <  5:116017/149‑1 (MQ=255)                                                                                                                          gCCTTTCCCGGTGGCCGCCGC‑GACACGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCGCTACATGGCCTGAAAGt  <  5:110801/149‑1 (MQ=255)                                                                                                                                       gccgccGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCg                            >  3:41737/1‑110 (MQ=255)                                                                                                                                       gccgccGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCg                            <  4:41737/110‑1 (MQ=255)                                                                                                                                             gc‑gACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCGCTAc               <  7:319763/117‑1 (MQ=255)                                                                                                                                             gc‑gACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCGCTAc               >  8:319763/1‑117 (MQ=255)                                                                                                                                                   |                                                                                                                             ATTGGAAACCGACAGGCGGTTCCCCGAAGCGGCGGTCCTCTTGCCCATTACCCGCAGCGAAGCGCCCGAACTGGTCTTGACCCTGCGCGCCAAGGGCCTGTCTACCCACGGTGGCGAAGTAGCCTTTCCCGGTGGCCGCCGC‑GACCCGGAAGACCCGGACCTGGTCTTTACTGCCCTGCGCGAGGCTGAAGAAGAAATCGGCCTGCCCCCTGGCCTGGTGGAGGTGATTGGTCCGCTCAGCCCGCTGATTTCGCTACATGGCCTGAAAGT  >  NC_002947/1659938‑1660207

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