Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,706,189	C→T	7.2%	L101L (CTG→CTA)	PP_4162 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,706,189	0	C	T	7.2%	159.6 / 2.1	55	L101L (CTG→CTA)	PP_4162	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base T (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

GCGCCAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGACCGCCCA  >  NC_002947/4706044‑4706337                                                                                                                                                  |                                                                                                                                                      gcgcCAGGAACCACCTTCCCACCTGCCCAAAGAAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTGAGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAACGAAAGGCACCCAGCCAGGCCGGTGCATCGGGCAgg                                                                                                                                                    <  1:260585/149‑1 (MQ=255) gcgcCAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAgg                                                                                                                                                    >  6:58886/1‑149 (MQ=255)   gcCAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGc‑c                                                                                                                                                 >  3:213785/1‑149 (MQ=255)    ccAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑Cg                                                                                                                                                >  3:83425/1‑149 (MQ=255)      aGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTc                                                                                                                                              <  2:254619/149‑1 (MQ=255)           ccAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAg                                                                                                                                                     >  4:110617/1‑138 (MQ=255)           ccAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAg                                                                                                                                                     <  3:110617/138‑1 (MQ=255)            cAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTg                                                                                                                                        <  3:283929/149‑1 (MQ=255)                 gCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGt                                                                                                                                   >  7:203488/1‑149 (MQ=255)                  cccACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGt                                                                                                                                   <  4:213785/148‑1 (MQ=255)                  cccACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTc                                                                                                                                  >  8:276309/1‑149 (MQ=255)                   ccACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTca                                                                                                                                 <  6:42510/149‑1 (MQ=255)                          cccAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTcacgcac                                                                                                                            <  6:306259/147‑1 (MQ=255)                                   aGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTa                                                                                                                  <  7:276309/148‑1 (MQ=255)                                          aTCGAGGTAAGCAAGGGTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTcgc                                                                                                           <  7:135059/148‑1 (MQ=255)                                             gAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAAt                                                                                                       >  3:114485/1‑149 (MQ=255)                                              aGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAAtg                                                                                                      >  7:231310/1‑149 (MQ=255)                                              aGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAAtg                                                                                                      <  5:58886/149‑1 (MQ=255)                                                 tAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGa                                                                                                   <  8:227409/149‑1 (MQ=255)                                                   aGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATc                                                                                                 <  6:137857/149‑1 (MQ=255)                                                           tGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGc                                                                                                                      <  6:104366/120‑1 (MQ=255)                                                            gcgcCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGtt                                                                                        >  7:253645/1‑149 (MQ=255)                                                                 tttCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAAcc                                                                                   >  3:180932/1‑149 (MQ=255)                                                                                   ccAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCgg                                                                   >  6:60774/1‑147 (MQ=255)                                                                                     aaCCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGt                                                               >  4:159950/1‑149 (MQ=255)                                                                                       ccGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGc                                                             >  4:188798/1‑149 (MQ=255)                                                                                             aCCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAAccc                                                                                  <  1:35066/122‑1 (MQ=255)                                                                                             aCCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAAccc                                                                                  >  2:35066/1‑122 (MQ=255)                                                                                               cAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCTGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATc                                                                                                 <  3:274125/105‑1 (MQ=255)                                                                                               cAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATc                                                                                                 >  4:274125/1‑105 (MQ=255)                                                                                                         ccAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGc                                           <  3:188798/149‑1 (MQ=255)                                                                                                           aaCGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCa                                                     >  1:57302/1‑137 (MQ=255)                                                                                                           aaCGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCa                                                     <  2:57302/137‑1 (MQ=255)                                                                                                                  ggAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcag                                  <  4:180932/149‑1 (MQ=255)                                                                                                                   gAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcagc                                 <  5:60774/149‑1 (MQ=255)                                                                                                                     aTGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTcagcag                               <  5:136256/149‑1 (MQ=255)                                                                                                                        ccACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGCTCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCCCGATCAGGATCCAGGCCATGGTCAGCAGGcc                            >  7:73953/1‑149 (MQ=255)                                                                                                                        ccACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGc                             <  8:203488/148‑1 (MQ=255)                                                                                                                         cACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCATAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCcg                           >  3:256146/1‑149 (MQ=255)                                                                                                                          aCCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCcac                                                                           <  8:282223/100‑1 (MQ=255)                                                                                                                          aCCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCcac                                                                           >  7:282223/1‑100 (MQ=255)                                                                                                                            ccagccagGCCGGTGCATCGGGTAGGCGCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCcggc                         >  1:8560/1‑149 (MQ=255)                                                                                                                            ccagccagGCCGGTGCATCGGGTAGGCGCGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCcggc                         >  1:13272/1‑149 (MQ=255)                                                                                                                            ccagccagGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCcggcg                        >  5:41281/1‑149 (MQ=255)                                                                                                                              agccagGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGt                                     <  4:229077/134‑1 (MQ=255)                                                                                                                              agccagGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGt                                     >  3:229077/1‑134 (MQ=255)                                                                                                                               gccagGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCggcggc                      >  2:184970/1‑148 (MQ=255)                                                                                                                                 cagGCCGGTGCATCGGGTAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAg                   <  2:8560/149‑1 (MQ=255)                                                                                                                                 cagGCCGGTGCATCGGGTAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAg                   <  2:13272/149‑1 (MQ=255)                                                                                                                                    gCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCcggcgg                       <  2:164977/142‑1 (MQ=255)                                                                                                                                    gCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCcggcgg                       >  1:164977/1‑142 (MQ=255)                                                                                                                                         tGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTcgc           <  4:83425/149‑1 (MQ=255)                                                                                                                                           cATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGa         <  6:41281/149‑1 (MQ=255)                                                                                                                                             tCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGAcc       >  8:202743/1‑149 (MQ=255)                                                                                                                                                 gCAGGT‑CGACGGCCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGACCGcc    >  1:296047/1‑148 (MQ=255)                                                                                                                                                  cAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGACCGCCca  <  3:159950/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GCGCCAGCAACCAGTTGCCCACCTGCCCAAAGTAAGGCTTGATCGAGGTAAGCAAGGCTGCGCCTTTCTGGTCCAGCGATTGCCAACCGTTGACCAGACGCTCGCCAACGAAAGGAATGCCACCCAGCCAGGCCGGTGCATCGGGCAGGC‑CGTCGACCTGCACGTCACGCACGAACGCCGTAGCGTCGCGAATGTGATCGGCCAGGTTGAACCCCAGCCACACCAGCGGCAGTGCCACGATCAGGATCCAGGCCATGGTCAGCAGGCCGGCGGCCAGCGTCTCGCGACCGCCCA  >  NC_002947/4706044‑4706337

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