Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,884	C→G	7.7%	R157G (CGC→GGC)	recN →	DNA repair/recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,884	0	C	G	7.7%	162.6 / 4.4	52	R157G (CGC→GGC)	recN	DNA repair/recombination protein
Reads supporting (aligned to +/- strand):  ref base C (28/20);  new base G (2/2);  total (30/22)
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.23e-01

AACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCAC  >  NC_002947/5378737‑5379030                                                                                                                                                     |                                                                                                                                                   aaCGGCACGCCCTGCCCGCTCGGCGTCCTCAAGGCGGTGGGCGAACTGGTGATCGATATCCACAGCCAGCACCAGCACCAGTCCCTGGTCAAGACCGACCTCCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGccc                                                                                                                                                    <  1:75002/149‑1 (MQ=255)         gCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGGCGGCGCCGTCGACCTGGCCCGCCAGGTa                                                                                                                                            >  5:284906/1‑148 (MQ=255)         gCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTa                                                                                                                                            >  3:290735/1‑148 (MQ=255)          gccTGCCCGCTCGGCGACCTCAAGGCGCTGGGCTAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCCCCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCGGGCCAGGtcga                                                                                                                                          <  7:117838/148‑4 (MQ=255)            cTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGc                                                                                                                                        <  8:125519/149‑1 (MQ=255)                ccGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGc                                                                                                                                    <  5:112560/149‑1 (MQ=255)                  gCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAg                                                                                                                                  >  5:269501/1‑149 (MQ=255)                   cTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGc                                                                                                                                 <  6:284906/149‑1 (MQ=255)                            cTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCGGGCCAGGtcgac                                                                                                                                         <  1:207384/132‑5 (MQ=255)                            cTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCGGGCCAGGtcgac                                                                                                                                         >  2:207384/1‑128 (MQ=255)                            cTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTg                                                                                                                        >  8:193919/1‑149 (MQ=255)                               aaGGCGCTGGGCGAACTGCTGATCCATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGACGCCAGCTTGACGAGTAAGCCCGCGCCGTCGACCTCGCCCGCCAGGAACAGCTTGCAGCCAAGCGCTGGa                                                                                                                      >  3:76722/1‑148 (MQ=255)                                     cTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGAc                                                                                                               <  6:269501/149‑1 (MQ=255)                                            aaCTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctg                                                                                                        <  4:290735/149‑1 (MQ=255)                                                 ctgATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctggagc                                                                                                    <  2:275503/148‑1 (MQ=255)                                                    aTCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctggagctg                                                                                                  >  6:271150/1‑147 (MQ=255)                                                            ccACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctggagctggag                                                                                               >  6:59866/1‑142 (MQ=255)                                                            ccACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctggagctggag                                                                                               <  5:59866/142‑1 (MQ=255)                                                                  ccagcACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACt                                                                                  >  8:91293/1‑149 (MQ=255)                                                                          aGCACCAGTCCCTGCTCATGACCGA‑CTCCCACCGGCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCGAGGTACAGCTTGCAGCCATGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGcgac                                                                          <  4:76722/149‑1 (MQ=255)                                                                              ccAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTCCTTGACGAGTATGCAGGCGCCGTCGACCTGGCCCGCCAGGAACAGCTTGCACCCAAGCGCTGGAATCAGACCAGCCTGGAGCAGGAGCGCCTCTCCAACTCGGGcgacgat                                                                       >  6:303534/1‑147 (MQ=255)                                                                               cAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGagca                                                                     <  7:91293/149‑1 (MQ=255)                                                                                    ccTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGc                                                                                                                                        <  1:106384/77‑1 (MQ=255)                                                                                    ccTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGc                                                                                                                                        >  2:106384/1‑77 (MQ=255)                                                                                    ccTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctggag                                                                                                     >  4:305590/1‑112 (MQ=255)                                                                                    ccTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctggag                                                                                                     <  3:305590/112‑1 (MQ=255)                                                                                    ccTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAgcgcg                                                                >  8:90666/1‑149 (MQ=255)                                                                                    ccTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAgcgcg                                                                >  7:137340/1‑149 (MQ=255)                                                                                        cTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGa                                                                                                                <  1:83813/97‑1 (MQ=255)                                                                                        cTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGa                                                                                                                >  2:83813/1‑97 (MQ=255)                                                                                        cTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGcgcc                                                              <  5:271150/147‑1 (MQ=255)                                                                                        cTCAACACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGcgccc                                                             <  4:163167/148‑1 (MQ=255)                                                                                                cGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAgc                                                    >  7:85796/1‑149 (MQ=255)                                                                                                 gA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAgct                                                   >  6:100953/1‑149 (MQ=255)                                                                                                         accGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGag                                                                       >  8:109334/1‑122 (MQ=255)                                                                                                         accGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGag                                                                       <  7:109334/122‑1 (MQ=255)                                                                                                          ccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGtt                                           <  5:100953/149‑1 (MQ=255)                                                                                                                 cTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCt                                    >  5:74211/1‑149 (MQ=255)                                                                                                                               tGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGAc                       >  2:234024/1‑148 (MQ=255)                                                                                                                               tGCCGGCGCCGGCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACa                      <  7:222628/149‑1 (MQ=255)                                                                                                                                 ccGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAAc                    >  3:55704/1‑149 (MQ=255)                                                                                                                                      gcCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCgg               >  6:8139/1‑149 (MQ=255)                                                                                                                                       ccGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCACCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCa                                              <  4:203730/117‑1 (MQ=255)                                                                                                                                       ccGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCACCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCa                                              >  3:203730/1‑117 (MQ=255)                                                                                                                                         gTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGcgcccgc                                                           >  6:303786/1‑102 (MQ=255)                                                                                                                                         gTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGcgcccgc                                                           <  5:303786/102‑1 (MQ=255)                                                                                                                                            taCCTGGCGGGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTggg         >  2:147163/2‑149 (MQ=255)                                                                                                                                             aCCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGc        >  6:163611/1‑149 (MQ=255)                                                                                                                                             aCCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGc        >  6:145334/1‑149 (MQ=255)                                                                                                                                              ccTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGc        >  4:233258/1‑148 (MQ=255)                                                                                                                                               cTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGc        >  5:266602/1‑147 (MQ=255)                                                                                                                                                  gCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCa   <  7:90666/149‑1 (MQ=255)                                                                                                                                                  gCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCa   >  6:168234/1‑149 (MQ=255)                                                                                                                                                    ccGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCAc  >  5:308500/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGA‑CACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGAGCAC  >  NC_002947/5378737‑5379030

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