Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,694,030	C→A	7.9%	R713L (CGA→CTA)	ligD ←	DNA ligase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,694,030	0	C	A	7.9%	185.8 / 3.9	59	R713L (CGA→CTA)	ligD	DNA ligase D
Reads supporting (aligned to +/- strand):  ref base C (20/34);  new base A (2/3);  total (22/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGG  >  NC_002947/3693884‑3694163                                                                                                                                                   |                                                                                                                                      cgcCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGc                                                                                                                                     <  3:18240/149‑1 (MQ=255)    ccTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCgg                                                                                                                                  <  4:220185/149‑1 (MQ=255)     cTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCggg                                                                                                                                 >  6:73327/1‑149 (MQ=255)        gCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCGAGCCGGGtg                                                                                                                               <  7:198236/148‑2 (MQ=255)        gCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCGAGCCGGGt                                                                                                                                <  3:163077/147‑1 (MQ=255)           gTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAgg                                                                                                                           >  4:27563/1‑149 (MQ=255)           gTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCACGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGAAGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAgg                                                                                                                           >  4:22520/1‑149 (MQ=255)                  aGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAAc                                                                                                                    <  2:301976/149‑1 (MQ=255)                     tAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCAc                                                                                                                 <  4:336752/149‑1 (MQ=255)                       gTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGACAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                               >  7:309287/1‑149 (MQ=255)                          gATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                               <  1:64139/146‑1 (MQ=255)                          gATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGt                                                                                                               >  2:64139/1‑146 (MQ=255)                          gATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGaa                                                                                                            <  8:264120/149‑1 (MQ=255)                          gATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGa                                                                                                             >  8:325240/1‑148 (MQ=255)                          gATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGa                                                                                                             >  8:112588/1‑148 (MQ=255)                            tGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATa                                                                                                          <  5:73327/149‑1 (MQ=255)                               aaaTGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATAccc                                                                                                       <  7:298217/149‑1 (MQ=255)                                 aTGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     <  7:112588/149‑1 (MQ=255)                                 aTGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     <  2:17820/149‑1 (MQ=255)                                 aTGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     <  1:95268/149‑1 (MQ=255)                                 aTGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     <  8:22187/149‑1 (MQ=255)                                 aTGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     <  3:22520/149‑1 (MQ=255)                                 aTGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     <  3:27563/149‑1 (MQ=255)                                     gtaccACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCATTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGcca                                                                                                 >  2:150427/4‑149 (MQ=255)                                     gTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGc                                                                                                   >  7:37134/1‑147 (MQ=255)                                     gTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGc                                                                                                   <  8:37134/147‑1 (MQ=255)                                      tCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGccac                                                                                                <  8:309287/149‑1 (MQ=255)                                        ccACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     >  1:94966/1‑142 (MQ=255)                                        ccACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     >  1:99900/1‑142 (MQ=255)                                        ccACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     <  2:94966/142‑1 (MQ=255)                                        ccACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCtt                                                                                                     <  2:99900/142‑1 (MQ=255)                                        ccACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGccacca                                                                                              <  2:49315/149‑1 (MQ=255)                                         cACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGCGCACAACCACGTGAATACCCTTGCcaccac                                                                                             <  1:30488/149‑1 (MQ=255)                                          aCGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCGAGCCGGGtgg                                                                                                                              >  8:296322/1‑113 (MQ=255)                                          aCGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCGAGCCGGGtgg                                                                                                                              <  7:296322/115‑3 (MQ=255)                                            gcgATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTgg                                                                                          >  2:99351/1‑149 (MQ=255)                                                tttttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGtttt                                                                                      <  6:317170/149‑1 (MQ=255)                                                   ttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCa                                                                                    <  7:325240/148‑1 (MQ=255)                                                   ttCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAg                                                                                   >  3:352139/1‑149 (MQ=255)                                                       gCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCTTGGCTGAAGTCATTGACCTCGTCCCACACGGCTCGACGGGTCAAGTGCACAACCACTAGAATACCCTTCCCAGCACTGGTTTTCAGGaa                                                                                >  4:205551/1‑148 (MQ=255)                                                       gCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGaa                                                                                >  2:364978/1‑148 (MQ=255)                                                         cccGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAaca                                                                              >  4:171911/1‑148 (MQ=255)                                                                  gcgcTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGcc                                                                    <  4:352139/149‑1 (MQ=255)                                                                                      cAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGa                                                <  3:209643/149‑1 (MQ=255)                                                                                        gCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAgg                                              <  1:364978/149‑1 (MQ=255)                                                                                            ggCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGt                                            <  1:99351/147‑1 (MQ=255)                                                                                                aGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTc                                                                                     >  6:33337/1‑102 (MQ=255)                                                                                                aGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTc                                                                                     <  5:33337/102‑1 (MQ=255)                                                                                                         aCGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGtt                                                                                        <  7:24327/90‑1 (MQ=255)                                                                                                         aCGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGCGCACAACCACGTGAATACCCTTGCCACCACTGGtt                                                                                        >  8:24327/1‑90 (MQ=255)                                                                                                          cGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTcc                            <  3:171911/149‑1 (MQ=255)                                                                                                              ggCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCcagca                                                      <  5:8038/119‑1 (MQ=255)                                                                                                              ggCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTAGCCGGGTCAAGGGCACAACCACGTGAATACCCTAGCCACCACTGGATTTCAGGAACACCTTCAGGCCCAGTACGTCcagca                                                      >  6:8038/1‑119 (MQ=255)                                                                                                                                gACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAgg                                                    >  7:240001/1‑103 (MQ=255)                                                                                                                                gACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAgg                                                    <  8:240001/103‑1 (MQ=255)                                                                                                                                      gTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagg  <  2:51506/147‑1 (MQ=255)                                                                                                                                      gTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGGCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGcagg  <  1:150427/147‑1 (MQ=255)                                                                                                                                            cccGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCATCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTccacc                         <  1:371727/118‑1 (MQ=255)                                                                                                                                            cccGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCATCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTccacc                         >  2:371727/1‑118 (MQ=255)                                                                                                                                                   |                                                                                                                                      CGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGATGGCGTGGCTGAAGTCCTTGACCTCGTCCCACCCGGCTCGCCGGGTCAAGGGCACAACCACGTGAATACCCTTGCCACCACTGGTTTTCAGGAACACCTTCAGGCCCAGTTCGTCCAGCAGGGTGAGGGTCAGTTGGGTCGCCTCCACCATGGCCTTCCAGGGCAGCGCAGG  >  NC_002947/3693884‑3694163

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