Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,836,055	C→T	8.0%	A79T (GCC→ACC)	yohC ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,836,055	0	C	T	8.0%	142.8 / 2.4	50	A79T (GCC→ACC)	yohC	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (20/26);  new base T (2/2);  total (22/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 = 4.57e-01

GGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAACGCCGAAACGGCCGGAATCGCCGCGAGGATCAGGGTG  >  NC_002947/1835907‑1836191                                                                                                                                                      |                                                                                                                                         gggTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGc                                                                                                                                          <  4:309494/149‑1 (MQ=255)    tAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAg                                                                                                                                       <  3:13377/149‑1 (MQ=255)      aaGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCTTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAgg                                                                                                                                      <  5:206833/148‑1 (MQ=255)         cgccgcCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAgg                                                                                                                                  <  1:364286/149‑1 (MQ=255)           ccgccAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGa                                                                                                                                 >  6:43203/1‑148 (MQ=255)           ccgccAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGAc                                                                                                                                <  5:340907/149‑1 (MQ=255)             gccAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACat                                                                                                                              <  5:43203/149‑1 (MQ=255)                 ggCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACatgatg                                                                                                                          <  5:154036/149‑1 (MQ=255)                  gccgccGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGc                                                                                                                         <  4:243149/149‑1 (MQ=255)                              cAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACatgatg                                                                                                                          <  8:249527/136‑1 (MQ=255)                              cAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACatgatg                                                                                                                          >  7:249527/1‑136 (MQ=255)                              cAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTg                                                                                                             <  2:98567/149‑1 (MQ=255)                               aaGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCCAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTgc                                                                                                            >  3:50462/1‑149 (MQ=255)                                     gTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTc                                                                                                      >  6:132343/1‑149 (MQ=255)                                      tGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTcc                                                                                                     <  5:85960/149‑1 (MQ=255)                                        gCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCAt                                                                                                   >  7:135647/1‑149 (MQ=255)                                         cGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATg                                                                                                  <  5:290063/149‑1 (MQ=255)                                             gtAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGt                                                                                                >  2:316779/1‑147 (MQ=255)                                                aGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGc                                                                                            <  8:316420/148‑1 (MQ=255)                                                 gTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGTTTGCGCtt                                                                                                        <  8:164508/135‑1 (MQ=255)                                                 gTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGTTTGCGCtt                                                                                                        >  7:164508/1‑135 (MQ=255)                                                  tGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCAt                                                                                          >  1:280223/1‑148 (MQ=255)                                                         aGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCActgctg                                                                                  >  1:116788/1‑149 (MQ=255)                                                            gATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTggtc                                                                               <  2:280223/149‑1 (MQ=255)                                                                 aTTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACGACGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTggtcggtc                                                                           >  1:281054/1‑149 (MQ=255)                                                                       ccATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATc                                                                    <  1:316779/149‑1 (MQ=255)                                                                        cATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCa                                                                   >  2:19629/1‑149 (MQ=255)                                                                               gggTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGcc                                                            >  7:247131/1‑149 (MQ=255)                                                                                     gCATCGTAGCTGCGGCCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCAATGCTGGTCGGTCACCGATCACCCAGCCCACCTg                                                      >  2:38039/1‑149 (MQ=255)                                                                                      cATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCGGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTgg                                                     >  3:86659/1‑149 (MQ=255)                                                                                      cATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTgg                                                     <  6:161786/149‑1 (MQ=255)                                                                                      cATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTgg                                                     <  1:19629/149‑1 (MQ=255)                                                                                              gTGCGGGCCATCCAGTGGATGAACG‑CGCCGATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCGATGGTCAGGATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGcc                                             <  1:38039/149‑1 (MQ=255)                                                                                               tGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGcc                                             >  3:91658/1‑148 (MQ=255)                                                                                                gCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGa                                           >  6:55777/1‑149 (MQ=255)                                                                                                 cGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGa                                           <  5:132343/148‑1 (MQ=255)                                                                                                 cGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGAt                                          <  2:116788/149‑1 (MQ=255)                                                                                                     ccATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAAc                                      >  7:113865/1‑149 (MQ=255)                                                                                                     ccATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAAc                                      <  8:224524/149‑1 (MQ=255)                                                                                                             tGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCGGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAACGCCGAAAc                              <  4:86659/149‑1 (MQ=255)                                                                                                                aTGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAACGTCGAAACGGc                           >  7:73875/1‑149 (MQ=255)                                                                                                                  gAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGTCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTggg                                                    >  5:43178/1‑122 (MQ=255)                                                                                                                  gAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGTCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTggg                                                    <  6:37853/122‑1 (MQ=255)                                                                                                                  gAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGTCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTggg                                                    <  6:43178/122‑1 (MQ=255)                                                                                                                  gAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGTCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTggg                                                    >  5:37853/1‑122 (MQ=255)                                                                                                                    aCG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAACGCCGAAACGGCCgg                        >  4:328123/1‑148 (MQ=255)                                                                                                                      g‑cgcCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAAc                                      >  7:330168/1‑132 (MQ=255)                                                                                                                      g‑cgcCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAAc                                      <  8:330168/132‑1 (MQ=255)                                                                                                                                 ccGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAACGCCGAAACGGCCGGAATCGCCGCGAg            <  5:55777/148‑1 (MQ=255)                                                                                                                                       cTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAACGCCGAAACGGCCGGAATCGCCGCGAGGATCAgg     <  2:281054/149‑1 (MQ=255)                                                                                                                                          cGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAACGCCGAAACGGCCGGAATCGCCGCGAGGATCAGGgtg  <  8:247131/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                         GGGTAAAGCGCCGCCAGGCCGCCGATGAACAAGGGGGTGGCGGTGTAAGTGGCAAAAGCGATGCATTGCGCCATGGACGGGTTGGCATCGTAGGTGCGGGCCATCCAGTGGATGAACG‑CGCCCATTACCGCCACTCCGGCGAGCATGGCCAGGTAGGACATGATGCTCATCCAGATTGCGCTTTCCATGGTCAGCATCACTGCTGGTCGGTCACCGATCACCCAGCCCACCTGGGTAGTGCCGATAAACGCCGAAACGGCCGGAATCGCCGCGAGGATCAGGGTG  >  NC_002947/1835907‑1836191

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