Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	23,833	T→A	7.1%	F220I (TTT→ATT)	PP_0018 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	23,833	0	T	A	7.1%	168.7 / 3.2	56	F220I (TTT→ATT)	PP_0018	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/27);  new base A (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TCTAGATCCGGAGTATCTGCACATCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGACTGCTAACTTTGCTTTTAGGACAGTCGCGA  >  NC_002947/23688‑23980                                                                                                                                                  |                                                                                                                                                    tCTAGATCCGGAGTATCTGCACATCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTtg                                                                                                                                                  >  5:144714/1‑149 (MQ=255)  cTAGATCCGGAGTATCTGCACATCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTtgt                                                                                                                                                 >  8:36614/1‑149 (MQ=255)     gATCCGGAGTATCTGCACATCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGc                                                                                                                                               <  7:120477/148‑1 (MQ=255)         cgctgattCTGCACATCAGTTGCTGGGGGCTGACT‑GTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAATTTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGtt                                                                                                                                          <  1:316507/141‑1 (MQ=255)            aGTATCTGCACATCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGc                                                                                                                                       <  7:226176/149‑1 (MQ=255)                     acaTCAGATGCTGGCGGCAGACTGGTATGTGGATACCTGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGt                                                                                                                              <  5:312608/149‑1 (MQ=255)                     acaTCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGt                                                                                                                              <  7:4285/149‑1 (MQ=255)                       aTCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGc                                                                                                                            >  6:6035/1‑149 (MQ=255)                       aTCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCAGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTg                                                                                                                             >  8:262191/1‑148 (MQ=255)                        tCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGaa                                                                                                                                    >  8:52380/1‑140 (MQ=255)                        tCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGaa                                                                                                                                    <  7:52380/140‑1 (MQ=255)                                 ggcggcAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGa                                                                                                                   <  1:119801/148‑1 (MQ=255)                                       aGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGAt                                                                                                            <  6:243340/149‑1 (MQ=255)                                             gTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCAt                                                                                                      >  5:326938/1‑149 (MQ=255)                                             gTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCAt                                                                                                      >  5:259547/1‑149 (MQ=255)                                                ctatcttACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTg                                                                                                    <  5:201018/142‑1 (MQ=255)                                                  tgGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGAt                                                                                                            <  1:272991/138‑1 (MQ=255)                                                  tgGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGAt                                                                                                            >  2:272991/1‑138 (MQ=255)                                                    gATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGc                                                                                               <  6:144714/149‑1 (MQ=255)                                                       aCCGGAAGATTTGGTCGAGCAGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTg                                                                                                    >  6:201018/1‑141 (MQ=255)                                                        ccGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCa                                                                                           <  6:259547/149‑1 (MQ=255)                                                           gAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGACTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGt                                                                                         >  4:66998/1‑148 (MQ=255)                                                            aaGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTg                                                                                                                             <  4:109869/111‑1 (MQ=255)                                                            aaGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTg                                                                                                                             >  3:109869/1‑111 (MQ=255)                                                                tttGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGcc                                                                                    >  6:138658/1‑148 (MQ=255)                                                                         gCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTgg                                                                          <  7:262191/149‑1 (MQ=255)                                                                         gCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTgg                                                                          >  6:204758/1‑149 (MQ=255)                                                                               cATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAAcc                                                                     >  1:329796/1‑148 (MQ=255)                                                                                 tGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCgag                                                                  <  5:204758/149‑1 (MQ=255)                                                                                 tGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCgag                                                                  <  7:253576/149‑1 (MQ=255)                                                                                  ggTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCgag                                                                  >  4:163299/1‑148 (MQ=255)                                                                                     tttCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCgag                                                                  >  1:7585/1‑145 (MQ=255)                                                                                       tCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGAt                                                             <  5:138658/148‑1 (MQ=255)                                                                                         gAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTg                                                          >  1:141815/1‑149 (MQ=255)                                                                                               gCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGGAGTTAGCTGAACGACTGGGAACCGAGGAGATCTGCACga                                                     <  6:326938/148‑1 (MQ=255)                                                                                               gCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACga                                                     <  1:80419/148‑1 (MQ=255)                                                                                                cTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACgaga                                                   <  3:66998/149‑1 (MQ=255)                                                                                                 tAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCAc                                                       <  3:120823/144‑1 (MQ=255)                                                                                                 tAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCAc                                                       >  4:120823/1‑144 (MQ=255)                                                                                                 tAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGCAGATCTGCACgaga                                                   >  4:157058/1‑148 (MQ=255)                                                                                                   aGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTc                                                <  3:163299/149‑1 (MQ=255)                                                                                                        tCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCAATTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACgg                                           >  3:128016/1‑149 (MQ=255)                                                                                                        tCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCAATTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACgg                                           >  3:132941/1‑149 (MQ=255)                                                                                                                ttCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCAt                                    >  6:326540/1‑148 (MQ=255)                                                                                                                   aTCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGa                                >  4:85921/1‑149 (MQ=255)                                                                                                                       aaGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGACt                              >  3:41859/1‑147 (MQ=255)                                                                                                                       aaGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGACt                              <  4:41859/147‑1 (MQ=255)                                                                                                                       aaGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGACTGc                            <  5:326540/149‑1 (MQ=255)                                                                                                                         gCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTc                                                <  5:258803/127‑1 (MQ=255)                                                                                                                         gCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTc                                                >  6:258803/1‑127 (MQ=255)                                                                                                                                  tcctcTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGACTGCTAACTTTGCtt                 >  5:243527/1‑149 (MQ=255)                                                                                                                                    ctctGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCgag                                                                  >  2:317029/1‑98 (MQ=255)                                                                                                                                    ctctGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCgag                                                                  <  1:317029/98‑1 (MQ=255)                                                                                                                                       tGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCATGAGACTCGACGGTTTGCATAAGACTGCTAACTTTGCTTTTa              <  7:36614/147‑1 (MQ=255)                                                                                                                                       tGCAAGGTGCAATTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGACTGCTAACTTTGCTTTTa              <  4:128016/147‑1 (MQ=255)                                                                                                                                       tGCAAGGTGCAATTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGACTGCTAACTTTGCTTTTa              <  4:132941/147‑1 (MQ=255)                                                                                                                                                  tttGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGACTGCTAACTTTGCTTTTAGGACAGTCGCGa  >  7:26109/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TCTAGATCCGGAGTATCTGCACATCAGATGCTGGCGGCAGACTGGTATGTGGATACCGGAAGATTTGGTCGAGCCGGCCATGGTTTTCGAAATTGCTAAGCAATCGAGGAATTCATCCAAGCATGGTGCTCCTCTGCAAGGTGCATTTGTGCCTGTTAGCGAAAATGGTGCTCTTCTCGATACGGATGCGCATTGGCCGCTGCAGTTAGCCGAACGACTGGGAACCGAGGAGATCTGCACGAGACTCGACGGTTTGCATAAGACTGCTAACTTTGCTTTTAGGACAGTCGCGA  >  NC_002947/23688‑23980

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