Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,998,856	C→G	7.7%	G44G (GGG→GGC)	PP_2622 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,998,856	0	C	G	7.7%	161.9 / 3.3	52	G44G (GGG→GGC)	PP_2622	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base G (2/2);  total (24/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 = 7.26e-01

CGTTGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCC  >  NC_002947/2998711‑2999004                                                                                                                                                   |                                                                                                                                                       cGTTGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCAcc                                                                                                                                                       <  7:342107/147‑1 (MQ=255)   ttGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTg                                                                                                                                                   >  7:26175/1‑149 (MQ=255)      aGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑Gc                                                                                                                                                <  2:331763/148‑1 (MQ=255)       gcgtgcgtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAAGGGCTGCTCTGGGAGCACCCCTG‑GC‑g                                                                                                                                              <  8:75477/149‑1 (MQ=255)           gcgtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCAt                                                                                                                                         >  1:348275/1‑149 (MQ=255)               aCCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTcag                                                                                                                                     >  1:294228/1‑149 (MQ=255)                 cTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTcag                                                                                                                                     >  2:97995/1‑147 (MQ=255)                 cTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTcag                                                                                                                                     >  2:173994/1‑147 (MQ=255)                         gcgGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTc                                                                                                                            >  1:358338/1‑148 (MQ=255)                          cgGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGa                                                                                                                          <  2:63699/149‑1 (MQ=255)                           ggCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAg                                                                                                                         >  5:119150/1‑149 (MQ=255)                             cTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGt                                                                                                                       <  1:9641/149‑1 (MQ=255)                             cTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGt                                                                                                                       <  1:97995/149‑1 (MQ=255)                                gAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTgg                                                                                                                     <  2:294228/148‑1 (MQ=255)                                  aGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGgcgc                                                                                                                  <  5:3030/149‑1 (MQ=255)                                       tGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGAt                                                                                                             >  3:183190/1‑149 (MQ=255)                                         gAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATgg                                                                                                           >  5:59070/1‑149 (MQ=255)                                           gggTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATggcg                                                                                                         <  1:177446/149‑1 (MQ=255)                                             gTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATggcggc                                                                                                       >  5:143074/1‑149 (MQ=255)                                                cggcggATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGaa                                                                                                    <  6:308936/149‑1 (MQ=255)                                                cggcggATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGaa                                                                                                    <  6:119150/149‑1 (MQ=255)                                                   cggATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGAGGGCGGCGAACtt                                                                                                 <  2:348275/149‑1 (MQ=255)                                                     gATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgt                                                                                               <  6:220718/149‑1 (MQ=255)                                                                tCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGTGCACCCCTG‑GC‑CGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGct                                                                                     >  7:121713/1‑148 (MQ=255)                                                                  aCCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCtctc                                                                                  <  5:365686/149‑1 (MQ=255)                                                                  aCCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCtctc                                                                                  <  6:59070/149‑1 (MQ=255)                                                                   ccTCCACAAACAAGTCGCCGTTGCTCAAGTCATGCCCATTGAAGTCGCGCAAGCCCAGCTC‑AGGGCTCCTCTCCCAGCACCCCTC‑GC‑GGGCATTCAGCAGCGTATCGAGGTGGCTCTAGATGGCGGCGAACTTGTGGCAACCGctctac                                                                                 >  4:274808/1‑147 (MQ=255)                                                                    cTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGACCTTGTGGGAAGCGCTCTccc                                                                                >  6:44227/1‑149 (MQ=255)                                                                        acaAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCg                                                                            <  1:173994/149‑1 (MQ=255)                                                                          aagcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCggg                                                                          <  6:143074/149‑1 (MQ=255)                                                                           agcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGa                                                                         <  5:354342/149‑1 (MQ=255)                                                                            gcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAg                                                                        <  7:251931/149‑1 (MQ=255)                                                                              aagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGc                                                                      <  2:358338/149‑1 (MQ=255)                                                                               agTCGCCGGTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCGGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGACGGGAgagc                                                                     <  7:286870/149‑4 (MQ=255)                                                                                   gCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCAtt                                                                                                                                        >  1:58463/1‑78 (MQ=255)                                                                                   gCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCAtt                                                                                                                                        <  2:58463/78‑1 (MQ=255)                                                                                                  aTGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGt                                                  <  2:302079/149‑1 (MQ=255)                                                                                                    gCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTc                                                 >  7:276932/1‑148 (MQ=255)                                                                                                            aaGTCGCGCAGGCCCAGCGC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCa                                         <  4:93861/148‑1 (MQ=255)                                                                                                             aGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTcca                                               >  4:183592/1‑141 (MQ=255)                                                                                                             aGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTcca                                               <  3:183592/141‑1 (MQ=255)                                                                                                               tCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGc                                      >  2:341147/1‑148 (MQ=255)                                                                                                                 gcgcAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGAt                                   >  2:356682/1‑149 (MQ=255)                                                                                                                       gCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaa                               >  2:175563/1‑147 (MQ=255)                                                                                                                           ccgcc‑AGGGGTGCTCCCAGAGCAGCCCTGAGCTGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGa                                                                         <  1:271675/99‑1 (MQ=21)                                                                                                                           ccgcc‑AGGGGTGCTCCCAGAGCAGCCCTGAGCTGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGa                                                                         >  2:271675/5‑103 (MQ=21)                                                                                                                             cTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCAc                       >  3:89845/1‑149 (MQ=255)                                                                                                                                 aGGGGTGCTCCCAGAGCAGCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaaaa                             >  1:102287/1‑140 (MQ=255)                                                                                                                                 aGGGGTGCTCCCAGAGCAGCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaaaa                             <  2:102287/140‑1 (MQ=255)                                                                                                                                   ggCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgt                                                                                               >  2:234705/1‑72 (MQ=255)                                                                                                                                   ggCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgt                                                                                               <  1:234705/72‑1 (MQ=255)                                                                                                                                             ggAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAg        <  1:341147/149‑1 (MQ=255)                                                                                                                                              gAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAgc       >  5:210516/1‑149 (MQ=255)                                                                                                                                                   ccccTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGccc  >  8:172577/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                       CGTTGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCGAGCCGCCC  >  NC_002947/2998711‑2999004

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