Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,068,520	Δ1 bp	7.7%	coding (686/948 nt)	PP_0926 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,068,520	0	C	.	7.7%	164.2 / 7.4	52	coding (686/948 nt)	PP_0926	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base . (2/2);  total (28/24)
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.15e-01

CCATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGGCGC  >  NC_002947/1068379‑1068666                                                                                                                                                 |                                                                                                                                                    ccATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAt                                                                                                                                              >  3:348249/1‑149 (MQ=255)       cTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGG‑CATCGATGGCCAGc                                                                                                                                       <  1:104225/149‑1 (MQ=255)             gCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTg                                                                                                                                 <  5:190190/149‑1 (MQ=255)                           cccAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCt                                                                                                                   <  8:230171/149‑1 (MQ=255)                                cTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCgc                                                                                                              <  7:302179/149‑1 (MQ=255)                                     tACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTc                                                                                                         >  4:228797/1‑149 (MQ=255)                                          aCCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGc                                                                                                      <  8:164687/147‑1 (MQ=255)                                            cTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGG‑CATCGATGGCCAGc                                                                                                                                       >  6:215065/1‑112 (MQ=255)                                            cTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGG‑CATCGATGGCCAGc                                                                                                                                       <  5:215065/112‑1 (MQ=255)                                                ggCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTgg                                                                                              <  4:348249/149‑1 (MQ=255)                                                     cTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGTGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCt                                                                                         >  2:254285/1‑149 (MQ=255)                                                      tACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTg                                                                                        >  1:219937/1‑149 (MQ=255)                                                      tACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTg                                                                                        >  1:258298/1‑149 (MQ=255)                                                         cTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGtt                                                                                     >  6:215896/1‑149 (MQ=255)                                                               caGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCgg                                                                               <  7:105773/149‑1 (MQ=255)                                                                 gATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGt                                                                             >  1:137673/1‑149 (MQ=255)                                                                      cccTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgc                                                                        >  4:38047/1‑149 (MQ=255)                                                                       ccTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgcg                                                                       <  2:137673/149‑1 (MQ=255)                                                                         tGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgcg                                                                       >  4:77485/1‑147 (MQ=255)                                                                         tGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgcg                                                                       <  3:77485/147‑1 (MQ=255)                                                                                ggCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGc                                                              >  8:41396/1‑149 (MQ=255)                                                                                   gcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGTGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGtt                                                                                     >  2:366599/1‑123 (MQ=255)                                                                                   gcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGTGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGtt                                                                                     <  1:366599/123‑1 (MQ=255)                                                                                   gcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTgg                                                           <  2:258298/149‑1 (MQ=255)                                                                                   gcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTgg                                                           >  1:293092/1‑149 (MQ=255)                                                                                     ggATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGTT‑GCCTGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGc                                                         <  3:274986/148‑1 (MQ=255)                                                                                       aTCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGCTACCCGCTGGCCATCGAT‑GCCACCTGCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTgg                                                           <  1:120895/148‑1 (MQ=255)                                                                                         cTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATcc                                                     >  4:202105/1‑149 (MQ=255)                                                                                         cTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATcc                                                     >  5:169070/1‑149 (MQ=255)                                                                                           gTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGCGCTGGGCATCCtt                                                   >  3:176238/1‑149 (MQ=255)                                                                                            tATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCttt                                                  >  2:42944/1‑149 (MQ=255)                                                                                            tATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCGGGGc                                                         >  4:274986/1‑142 (MQ=255)                                                                                             aTCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCttt                                                  <  1:42944/148‑1 (MQ=255)                                                                                              tCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTc                                                <  1:320325/149‑1 (MQ=255)                                                                                              tCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTc                                                >  1:352560/1‑149 (MQ=255)                                                                                              tCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTc                                                >  5:347630/1‑149 (MQ=255)                                                                                               cggcggCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTcc                                               >  7:31800/1‑149 (MQ=255)                                                                                                        gCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAg                                      <  3:38047/149‑1 (MQ=255)                                                                                                        gCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAg                                      <  3:202105/149‑1 (MQ=255)                                                                                                              gCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGcc                                <  3:350420/149‑1 (MQ=255)                                                                                                                      cgcTGTGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCt                        <  1:254285/149‑1 (MQ=255)                                                                                                                       gcTGCGTAACTGTGG‑T‑‑CGGCTGGCCACCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTg                       <  1:178459/149‑1 (MQ=255)                                                                                                                          gCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCgg                    >  3:270849/1‑149 (MQ=255)                                                                                                                            gTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGcc                  <  6:347630/149‑1 (MQ=255)                                                                                                                             tAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGc                   >  7:172428/1‑147 (MQ=255)                                                                                                                             tAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCa                 <  8:31800/149‑1 (MQ=255)                                                                                                                                 tgtgG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTa             >  4:280150/1‑149 (MQ=255)                                                                                                                                   tgG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACg           >  4:342981/1‑149 (MQ=255)                                                                                                                                     g‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGt          >  3:301241/1‑148 (MQ=255)                                                                                                                                        cgaagCTGG‑CATCGATGGCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGt          >  2:128684/5‑148 (MQ=255)                                                                                                                                            gCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATgg     >  2:206228/1‑149 (MQ=255)                                                                                                                                                gCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGgcgc  <  2:293092/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CCATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGGCGC  >  NC_002947/1068379‑1068666

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