Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	752,034	C→G	7.3%	F18L (TTC→TTG)	PP_5451 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	752,034	0	C	G	7.3%	170.9 / 3.1	55	F18L (TTC→TTG)	PP_5451	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base G (2/2);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

CGGCACTGGTCGTACGTTCCTTATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTCCGATGCAGAGCTTTACAAAGT  >  NC_002947/751890‑752180                                                                                                                                                 |                                                                                                                                                     cGGCACTGGTCGTACGTTCCTTATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATc                                                                                                                                                  <  6:261478/149‑1 (MQ=255)      cTGGTCGTACGTTCCTTATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAg                                                                                                                                             <  7:13569/149‑1 (MQ=255)             tACGTTCCTTATGTACTTGGCTTACTAGAAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAAt                                                                                                                                      >  6:195591/1‑149 (MQ=255)               cGTTCCTTATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATg                                                                                                                                     <  5:96182/148‑1 (MQ=255)                gTTCCTTATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGt                                                                                                                                   >  2:7957/1‑149 (MQ=255)                    cTTATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGa                                                                                                                               >  4:244936/1‑149 (MQ=255)                     ttATGTACTTTGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCt                                                                                                                                 <  2:35041/146‑1 (MQ=255)                     ttATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCATCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTc                                                                                                                                  <  4:252009/145‑1 (MQ=255)                     ttATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCt                                                                                                                                 >  1:35041/1‑146 (MQ=255)                     ttATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGa                                                                                                                               >  5:91994/1‑148 (MQ=255)                     ttATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGAAAGAATGGTc                                                                                                                                  >  3:252009/1‑145 (MQ=255)                               ggCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAg                                                                                                                    >  7:17945/1‑148 (MQ=255)                               ggCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAg                                                                                                                    >  3:135559/1‑148 (MQ=255)                                 cTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAt                                                                                                                 >  3:258197/1‑149 (MQ=255)                                  ttACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑c                                                                                                               >  2:220359/1‑149 (MQ=255)                                        ggAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCt                                                                                                         <  8:17945/149‑1 (MQ=255)                                          tgTGCTACAGCTTGGTGCCAGGCGCGTTTGGCCTTGGTAGGCGCGTGGAAATGTCGTTTCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCtt                                                                                                        <  8:167226/147‑1 (MQ=255)                                               aaCAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAATAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGa                                                                                                  >  4:39938/1‑149 (MQ=255)                                               aaCAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGa                                                                                                  >  7:148209/1‑149 (MQ=255)                                                         tCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGAt                                                                                                 <  3:111782/140‑1 (MQ=255)                                                         tCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGAt                                                                                                 >  4:111782/1‑140 (MQ=255)                                                              gCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGgcagc                                                                                   <  6:91994/149‑1 (MQ=255)                                                               ccGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTGGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGgcagca                                                                                  >  1:100381/1‑149 (MQ=255)                                                               ccGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGATAACTCACGCGATCAGCAGGGCAGATTTGGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGATgcagca                                                                                  >  1:104263/1‑149 (MQ=255)                                                                 gcgATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATg                                                                                >  1:163502/1‑149 (MQ=255)                                                                   gATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCATATTTCGATCAGAAGCAAGAATGGTCTGATACCGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGc                                                                               <  1:55762/149‑1 (MQ=255)                                                                   gATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATg                                                                                >  5:141242/1‑147 (MQ=255)                                                                      tCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATAAGCAGGGCAGATTTCGATCAGAAGCAAGAATGCTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGAAGCAAGCAGa                                                                            >  7:216751/1‑148 (MQ=255)                                                                           aTGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCa                                                                                                                  >  2:263182/1‑106 (MQ=255)                                                                                aGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGtcatc                                                                    <  6:210579/146‑1 (MQ=255)                                                                                aGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGtcatc                                                                    >  5:210579/1‑146 (MQ=255)                                                                                  gagaGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAATAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCt                                                               <  3:39938/149‑1 (MQ=255)                                                                                       gggAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGcc                                                          >  3:235691/1‑149 (MQ=255)                                                                                       gggAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGcc                                                          <  4:135559/149‑1 (MQ=255)                                                                                       gggAAATGTCGATGCTATCCTTGGAAGGACTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGcc                                                          >  1:33324/1‑149 (MQ=255)                                                                                           aaTGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTa                                                       >  3:126455/1‑148 (MQ=255)                                                                                           aaTGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTa                                                       <  4:126455/148‑1 (MQ=255)                                                                                             tGTCGATGCTATCCTTGGACGCGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTacac                                                    <  8:121815/149‑1 (MQ=255)                                                                                                 gATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAAAAAGTCATCAAGCTGTGCCGTACACGCAt                                                >  1:146147/1‑149 (MQ=255)                                                                                                        tCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGc                                         <  4:235691/149‑1 (MQ=255)                                                                                                               aaGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAgg                                  <  3:244936/149‑1 (MQ=255)                                                                                                                       aCTCACGCGATCAGCAGGGCAGATTTGGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGt                            <  2:104263/147‑1 (MQ=255)                                                                                                                       aCTCACGCGATCAGCAGGGCAGATTTGGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGt                            <  2:100381/147‑1 (MQ=255)                                                                                                                       aCTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAAAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGt                            <  2:146147/147‑1 (MQ=255)                                                                                                                           aCGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTcc                      >  7:201575/1‑149 (MQ=255)                                                                                                                            cgcgATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCATAGGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTcc                      >  1:77592/1‑149 (MQ=255)                                                                                                                            cgcgATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTCCg                     <  5:195591/149‑1 (MQ=255)                                                                                                                              cgATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTCCGAt                   <  4:258197/149‑1 (MQ=255)                                                                                                                               gATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTCCGATg                  <  8:148209/149‑1 (MQ=255)                                                                                                                                       gggCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGt                                                        <  6:238610/103‑1 (MQ=255)                                                                                                                                       gggCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGt                                                        >  5:238610/1‑103 (MQ=255)                                                                                                                                         gCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTCCGATGCAGAGCTTTa        <  1:7957/149‑1 (MQ=255)                                                                                                                                         gCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTCCGATGCAGAGCTTTa        <  2:33324/149‑1 (MQ=255)                                                                                                                                          cAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTCCGATGCAGAGCTTTa        >  7:79097/1‑148 (MQ=255)                                                                                                                                               ttCGATCAGAAGCAAGAATCGTCTGATA‑CGCCTCACCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGAATCAAGAAAGTCATCAAACTGTGCCGTACACGCATTGAAGGCTGACACGGCAGGTTGCTCCGATGCAGAGCATTACAAAGt  >  5:70550/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CGGCACTGGTCGTACGTTCCTTATGTACTTGGCTTACTAGGAGTGCAACAGCTGGCTCCCAGCCGCGATTCGCCATGGAAGGAGAGGGGAAATGTCGATGCTATCCTTGGAAGGGCTAACTCACGCGATCAGCAGGGCAGATTTCGATCAGAAGCAAGAATGGTCTGATA‑CGCCTCAGCAT‑CGGCGCTTGCTGGATGCTTACGAGGCAGCATGCAGAAAGTCATCAAGCTGTGCCGTACACGCATTGAAGGCTGAGAGGGCAGGTTGGTCCGATGCAGAGCTTTACAAAGT  >  NC_002947/751890‑752180

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