Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,884,963	C→A	8.1%	E139* (GAA→TAA)	yedI ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,884,963	0	C	A	8.1%	174.0 / 3.1	57	E139* (GAA→TAA)	yedI	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base A (3/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATC  >  NC_002947/4884817‑4885105                                                                                                                                                   |                                                                                                                                                acCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCa                                                                                                                                                >  6:81364/1‑148 (MQ=255)  cccagtgTCCTCACCACCGCAATGCCCGACAGCACGATGATCTGCTGGGTCAGCGGCTAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACATGATGCAGTCGGTACGCACCGCCCCCTTGATCTGGGCGTTTTCaaa                                                                                                                                              <  6:114393/143‑1 (MQ=255)   cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGACGACAGGATGAAGTCGGTACGCACCGCCACCTTGATCTTGGCTTTTTCaa                                                                                                                                               >  5:261848/1‑147 (MQ=255)         tcatcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGccacc                                                                                                                                       <  2:163375/149‑1 (MQ=255)         tcatcaCCACCGCAAAGCCGGAAAGCACGATGATCTCCGGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGGCGGTACGCACCGCCCCCTTGATCTTGGCTTTTGAAAAAGccaag                                                                                                                                       <  7:267587/149‑3 (MQ=255)              accaccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGt                                                                                                                                   >  4:42911/1‑148 (MQ=255)               ccaccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCa                                                                                                                                 <  8:296309/149‑1 (MQ=255)                  ccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGcc                                                                                                                                          >  5:82359/1‑137 (MQ=255)                  ccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGcc                                                                                                                                          <  6:82359/137‑1 (MQ=255)                    gCAATGCCGCAAAGCACGATGATCTGCTGTGTCAGCGGCGATTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACt                                                                                                                            <  6:261848/149‑1 (MQ=255)                     cAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTg                                                                                                                           <  6:4285/149‑1 (MQ=255)                         gCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCAt                                                                                                                       >  3:329156/1‑149 (MQ=255)                                   aCGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGc                                                                                                             <  5:10972/149‑1 (MQ=255)                                     gatgatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcct                                                                                                           <  1:175019/149‑1 (MQ=255)                                      atgatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcctc                                                                                                          <  6:280575/149‑1 (MQ=255)                                          tCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCATAAGccac                                                                                                                                        >  5:312069/1‑115 (MQ=255)                                          tCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGccac                                                                                                                                        <  6:312069/115‑1 (MQ=255)                                              ctgGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTgcg‑c                                                                                                  <  5:324090/148‑1 (MQ=255)                                                              tCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCt                                                                                  >  2:337845/1‑148 (MQ=255)                                                                aGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCTCCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCAt                                                                               >  1:43952/1‑149 (MQ=255)                                                                          cccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGGCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGGCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCtt                                                                     <  1:264770/149‑1 (MQ=255)                                                                            cAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGc                                                                   >  5:58380/1‑149 (MQ=255)                                                                                gTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCtgtg                                                               >  7:156666/1‑149 (MQ=255)                                                                                    tGGCGATGATTTCTGCCGAGAGGATGAAGTAGGTACGCACCGCCCCCTTGATCTTGGCTTTTTAAAAAGCCACCAGGTCAACACTGGCATCGTCAAAGGCCTCCTGGCT‑CGCCTCGTGCTGTGCCTAATCCTATTACTTGCTGTGc                                                              >  8:317484/1‑146 (MQ=255)                                                                                        gatgatTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcctc                                                                                                          >  1:144858/1‑99 (MQ=255)                                                                                        gatgatTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcctc                                                                                                          <  2:144858/99‑1 (MQ=255)                                                                                        gatgatTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACt                                                       >  3:91312/1‑149 (MQ=255)                                                                                             tttCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTtgtg                                                   <  3:42911/148‑1 (MQ=255)                                                                                             tttCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTgc                                                  <  6:58380/149‑1 (MQ=255)                                                                                                cTGCCGACAGGATGAATTCGGTACGCACCGCCCCCTTGCTCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTgcgc                                                <  7:317484/148‑1 (MQ=255)                                                                                                     gACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCtt                                          <  5:81364/149‑1 (MQ=255)                                                                                                      aCAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCg                                                                                            <  8:30790/98‑1 (MQ=255)                                                                                                      aCAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCg                                                                                            >  7:30790/1‑98 (MQ=255)                                                                                                        aGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCtg                                                                                       >  3:112942/1‑101 (MQ=255)                                                                                                        aGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCtg                                                                                       <  4:112942/101‑1 (MQ=255)                                                                                                          gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAAc                                                        >  5:104937/1‑130 (MQ=255)                                                                                                          gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAAc                                                        <  6:104937/130‑1 (MQ=255)                                                                                                          gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCg                                      >  6:239863/1‑148 (MQ=255)                                                                                                          gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCg                                      >  5:86243/1‑148 (MQ=255)                                                                                                          gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCg                                      >  7:278417/1‑148 (MQ=255)                                                                                                           aTGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTgc                                                  >  4:22536/1‑135 (MQ=255)                                                                                                           aTGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTgc                                                  <  3:22536/135‑1 (MQ=255)                                                                                                            tGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGaaa                                   <  1:337845/149‑1 (MQ=255)                                                                                                             gAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtc                                         <  6:233751/142‑1 (MQ=255)                                                                                                             gAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtc                                         >  5:233751/1‑142 (MQ=255)                                                                                                                 tCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCGTCGTGCTGGGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGACACCCTc                              >  4:314412/1‑149 (MQ=255)                                                                                                                        gCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGaaa                                   >  4:328869/1‑137 (MQ=255)                                                                                                                        gCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGaaa                                   <  3:328869/137‑1 (MQ=255)                                                                                                                         cACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGcaca                      <  8:156666/149‑1 (MQ=255)                                                                                                                                     gATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCa           <  5:239863/148‑1 (MQ=255)                                                                                                                                     cctggtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCAc          >  1:200430/6‑149 (MQ=255)                                                                                                                                        ggtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCAcc         >  4:322600/3‑147 (MQ=255)                                                                                                                                         ttGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGTCCCCTCGTGCTGTGCATCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGACCAGGAACGCGCCACCCa       >  1:277383/1‑149 (MQ=255)                                                                                                                                         ttGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAg      >  5:122094/1‑149 (MQ=255)                                                                                                                                         taggCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAg      >  3:249536/3‑149 (MQ=255)                                                                                                                                           ggCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGc     <  4:249536/148‑1 (MQ=255)                                                                                                                                             cTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATc  <  2:200430/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                ACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCG‑CGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATC  >  NC_002947/4884817‑4885105

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