Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,937	A→T	10.4%	I99F (ATC→TTC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,937	0	A	T	10.4%	168.6 / 8.3	57	I99F (ATC→TTC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base T (4/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAG  >  NC_002947/1549792‑1550084                                                                                                                                                  |                                                                                                                                                    cAAGGACGGCGTTTCCGTCGCCAAAGAAATCCAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGt                                                                                                                                                 >  2:877100/1‑150 (MQ=255)         gCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAg                                                                                                                                         <  1:1174532/150‑1 (MQ=255)           gTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGc                                                                                                                                       >  2:444352/1‑150 (MQ=255)                cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGag                                                                                                                                  <  1:444352/150‑2 (MQ=255)                 gTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                                 <  1:1233735/150‑1 (MQ=255)                      cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                   <  2:820465/143‑1 (MQ=255)                      cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACGATGGCCTGa                                                                                                                                   >  1:820465/1‑143 (MQ=255)                       aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                                                                                                           >  2:472524/1‑150 (MQ=255)                       aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTc                                                                                                                           <  1:156522/150‑1 (MQ=255)                        aaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCg                                                                                                                          >  1:933894/1‑150 (MQ=255)                                            gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                      >  1:536789/1‑150 (MQ=255)                                              tttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCACCGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGa                                                                                                    >  2:624353/1‑150 (MQ=255)                                              tttCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGa                                                                                                    >  1:552593/1‑150 (MQ=255)                                                tCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAcc                                                                                                  <  2:252893/150‑1 (MQ=255)                                                  gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTg                                                                                                >  2:881890/1‑150 (MQ=255)                                                   aaaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                                 <  2:111473/116‑1 (MQ=255)                                                   aaaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGaaa                                                                                                                                 >  1:111473/1‑116 (MQ=255)                                                   aaaaCATGGGCGCCCAGCTGGGCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGa                                                                                               >  1:1052279/1‑150 (MQ=255)                                                        aTGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAgcgc                                                                                          >  1:378146/1‑150 (MQ=255)                                                           ggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGTATGAACCCGATGGACCTGAAGCGCGGt                                                                                       >  1:757539/1‑150 (MQ=255)                                                                ccAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTCAAGTCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGa                                                                                  <  1:797610/150‑1 (MQ=255)                                                                   gCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGTATGAACCCGATGGACCTGAAGCGCGGTATCGACaa                                                                               <  2:757539/150‑1 (MQ=255)                                                                      ggTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                            <  1:881890/150‑1 (MQ=255)                                                                        tCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTTTCGACAAGGCCa                                                                          >  1:647255/1‑150 (MQ=255)                                                                        tCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCa                                                                          >  1:1119809/1‑150 (MQ=255)                                                                         cAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCAc                                                                         <  1:1376972/150‑1 (MQ=255)                                                                                gTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGa                                                                                  >  1:794177/1‑134 (MQ=255)                                                                                gTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGa                                                                                  <  2:794177/134‑1 (MQ=255)                                                                                 ttGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCcg                                                                 >  2:1165807/1‑150 (MQ=255)                                                                                        cAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGa                                                                                  <  1:1332784/126‑1 (MQ=255)                                                                                        cAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGa                                                                                  >  2:1332784/1‑126 (MQ=255)                                                                                        cAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCa                                                                          <  1:699671/134‑1 (MQ=255)                                                                                        cAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCa                                                                          >  2:699671/1‑134 (MQ=255)                                                                                                   acgCTGCCGGTGACGGCACCACCCATGCTACTGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaaga                                               <  2:536789/150‑1 (MQ=255)                                                                                                   acgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCGCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaaga                                               >  1:403558/1‑150 (MQ=255)                                                                                                     gCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAAc                                             >  1:608739/1‑150 (MQ=255)                                                                                                     gCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAAc                                             <  1:472524/150‑1 (MQ=255)                                                                                                       gcacGGTGACGGAACCACAACCGCTACCGTTCTGGCTCAGTCTATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGACCCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCt                                           <  2:265637/147‑1 (MQ=255)                                                                                                             tGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                                     <  2:608739/150‑1 (MQ=255)                                                                                                                      caccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTAt                                                                                     <  2:947178/93‑1 (MQ=255)                                                                                                                      caccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTAt                                                                                     >  1:947178/1‑93 (MQ=255)                                                                                                                      caccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATgcgc                            >  1:428516/1‑150 (MQ=255)                                                                                                                       accaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGcc                           >  1:243640/1‑150 (MQ=255)                                                                                                                        ccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATgc                              <  1:641915/146‑1 (MQ=255)                                                                                                                        ccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATgc                              >  2:641915/1‑146 (MQ=255)                                                                                                                            cGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAACCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTc                      >  2:1033933/1‑150 (MQ=255)                                                                                                                               tACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCaa                   >  1:1354119/1‑150 (MQ=255)                                                                                                                                aCCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAg                  <  2:1354119/150‑1 (MQ=255)                                                                                                                                 ccGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAgg                 <  2:933894/150‑1 (MQ=255)                                                                                                                                  cGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGc                >  1:671213/1‑150 (MQ=255)                                                                                                                                  cGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGc                <  1:1165807/150‑1 (MQ=255)                                                                                                                                   gTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGcc               >  1:145079/1‑150 (MQ=255)                                                                                                                                         tttCAGGCCTTCGTTGACGATGACCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGcc         >  1:207686/3‑150 (MQ=255)                                                                                                                                           tCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACGTGTCAAAGCAATGGGCCGACTCGGAAGCCATCGGCCa       >  2:682000/1‑150 (MQ=255)                                                                                                                                              ggCCATCGTCAACGAAGGCCTGAAAGCCGTCGCAGTCCGCATGAACCCGAGGCATCGCAAGCGCGCTATCGACACGTCCACCACTGCCATCGTTGCCAAAATGAAGAACCTCTCCTAGCCATGCGCCGACTCCAAGGCCATCGCCCcgga    <  2:364297/150‑2 (MQ=255)                                                                                                                                               gCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggta   >  1:1425203/1‑150 (MQ=255)                                                                                                                                                ccATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggtag  <  2:378146/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CAAGGACGGCGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAG  >  NC_002947/1549792‑1550084

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑