Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,288,843	G→A	7.2%	P94S (CCG→TCG)	PP_3760 ←	methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,288,843	0	G	A	7.2%	155.9 / 2.2	55	P94S (CCG→TCG)	PP_3760	methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base A (2/2);  total (30/25)
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.33e-01

TGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCAAGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAA  >  NC_002947/4288701‑4288990                                                                                                                                                 |                                                                                                                                                    tGCGGGTTGATGTCTGTGGCGTAGATGATGGTGCGTTCAAGCAGCCCGTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACCTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGcaccac                                                                                                                                               <  5:296680/149‑1 (MQ=255)       ttGATGTCGGTGGCGTAGATGATGGTGCGTTCAAGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGtt                                                                                                                                         <  3:314020/149‑1 (MQ=255)            gTCGGTGGCGTAGATGATGGTGCGTTCAAGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCa                                                                                                                                    >  4:339335/1‑149 (MQ=255)              cGGTGGCGTAGATGATGGTGCGTTCAAGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACg                                                                                                                                  >  7:259553/1‑149 (MQ=255)                 tGGCGTAGATGATGGTGCGTTCAAGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCTGACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCCTTgccgtgcgc                                                                                                                                <  1:74276/149‑8 (MQ=255)                            atgGTGCGTTCAAGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCCTTGCCGTGCGCAACGAAGTGGt                                                                                                                     <  1:218822/148‑1 (MQ=255)                                 gCGTTCAAGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGt                                                                                                               <  8:259553/149‑1 (MQ=255)                                           agcCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAg                                                                                                                  <  5:19931/136‑1 (MQ=255)                                           agcCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAg                                                                                                                  >  6:19931/1‑136 (MQ=255)                                              ccTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACAtctc                                                                                                  <  3:339335/149‑1 (MQ=255)                                               cTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCg                                                                                                 >  5:233590/1‑149 (MQ=255)                                               cTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCg                                                                                                 <  6:12954/149‑1 (MQ=255)                                                   cTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACAtctc                                                                                                  <  3:60414/144‑1 (MQ=255)                                                   cTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACAtctc                                                                                                  >  4:60414/1‑144 (MQ=255)                                                    tCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTg                                                                                                                       >  7:251620/1‑122 (MQ=255)                                                    tCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTg                                                                                                                       <  8:251620/122‑1 (MQ=255)                                                          aaCAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGCCAGGGTCACGGAACATCTCGCTGACCGGGAt                                                                                      >  8:113333/1‑149 (MQ=255)                                                              ggATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTc                                                                                  <  4:321942/149‑1 (MQ=255)                                                                          aaTACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTgcag                                                                       >  5:158848/1‑148 (MQ=255)                                                                          aaTACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGACACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTgcag                                                                       >  8:182963/1‑148 (MQ=255)                                                                          aaTACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGACACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTgcag                                                                       >  2:318880/1‑148 (MQ=255)                                                                             acacCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCCAAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCaa                                                                    >  1:274900/1‑149 (MQ=255)                                                                               acCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCgg                                                                                         <  7:325434/125‑1 (MQ=255)                                                                               acCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCgg                                                                                         >  8:325434/1‑125 (MQ=255)                                                                                cctcctCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTg                                                                 >  5:289098/1‑148 (MQ=255)                                                                                 ctcctcCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAg                                                                                >  2:123866/1‑132 (MQ=255)                                                                                    ctccCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATc                                                                                     <  6:254458/124‑1 (MQ=255)                                                                                    ctccCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATc                                                                                     >  5:254458/1‑124 (MQ=255)                                                                                                    cccGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTg                                            >  2:225164/1‑149 (MQ=255)                                                                                                     ccGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGACACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGc                                           <  7:182963/149‑1 (MQ=255)                                                                                                     ccGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGACACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGc                                           <  1:318880/149‑1 (MQ=255)                                                                                                     ccGCGATCCAGACTTTGATC‑GACGCCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAATGAAGTGGTCAGGGTCACGCAACATCTCGCTGACCGGGATCGTCACGTACTTCAGCAACTTCAAGAACATGCCCGCGTCGTg                                            >  7:339956/1‑148 (MQ=255)                                                                                                      cgcgATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTg                                            >  3:205719/1‑147 (MQ=255)                                                                                                      cgcgATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCa                                          <  2:340495/149‑1 (MQ=255)                                                                                                      cgcgATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCa                                          <  6:289098/149‑1 (MQ=255)                                                                                                      cgcgATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGCCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCa                                          <  7:113333/149‑1 (MQ=255)                                                                                                          aTCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGaa                                      >  1:309909/1‑148 (MQ=255)                                                                                                             cAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGAcgcg                                   >  8:55799/1‑149 (MQ=255)                                                                                                             cAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGAcgcg                                   >  8:101920/1‑149 (MQ=255)                                                                                                             cAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGAcgcg                                   >  3:316464/1‑149 (MQ=255)                                                                                                              aGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCAt                                                              >  1:230227/1‑121 (MQ=255)                                                                                                              aGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCAt                                                              <  2:230227/121‑1 (MQ=255)                                                                                                              aGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACgcgc                                  >  2:46873/1‑149 (MQ=255)                                                                                                              aGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGAGCAGGACgagc                                  >  6:2806/1‑149 (MQ=255)                                                                                                                  tttGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCt                              >  3:51511/1‑149 (MQ=255)                                                                                                                    tGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTg                             <  4:51511/148‑1 (MQ=255)                                                                                                                    tGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTg                             >  3:139053/1‑148 (MQ=255)                                                                                                                     gATC‑GACGGCG‑AGGTGCGCTGCAGCGGCACGTCTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCAGGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCATCTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGc                            <  8:339956/148‑1 (MQ=255)                                                                                                                           aCGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACg                        >  6:23268/1‑147 (MQ=255)                                                                                                                             ggCC‑AGGGGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGGCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGGGCAGGACGCGCTcc                               <  5:193928/138‑1 (MQ=255)                                                                                                                               cc‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAAc                                                          <  2:273963/109‑1 (MQ=255)                                                                                                                               cc‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAAc                                                          >  1:273963/1‑109 (MQ=255)                                                                                                                                    gTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGt              >  6:48991/1‑149 (MQ=255)                                                                                                                                       cgcAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAg           <  7:55799/149‑1 (MQ=255)                                                                                                                                       cgcAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAg           <  7:101920/149‑1 (MQ=255)                                                                                                                                                ggCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGaa  <  4:205719/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCAAGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATC‑GACGGCC‑AGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAA  >  NC_002947/4288701‑4288990

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