Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	166,856	A→T	10.7%	A17A (GCT→GCA)	pntAA ←	pyridine nucleotide transhydrogenase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	166,856	0	A	T	10.7%	143.7 / 10.5	56	A17A (GCT→GCA)	pntAA	pyridine nucleotide transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base T (4/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-01

CCGTAGGCATCGGCGGCGCTCCCAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGGGCGGCTGTGGAGCGCCAACCCGCCGGATCCCG  >  NC_002947/166713‑167004                                                                                                                                                |                                                                                                                                                      ccGTAGGCATCGGCGGCGCTCCCAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGAcc                                                                                                                                                 >  2:937990/1‑150 (MQ=255)    tAGGCATCGGCGGCGCTCCCAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCtt                                                                                                                                              <  1:417750/150‑1 (MQ=255)    tAGGCATCGGCGGCGCTCCCAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCtt                                                                                                                                              <  1:423835/150‑1 (MQ=255)       gCATCGGCGGCGCTCCCAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGtc                                                                                                                                           >  2:924485/1‑150 (MQ=255)         aTCGGCGGCGCTCAGCAATTAAGCACCCACGGCTTCAGAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGtctc                                                                                                                                         >  1:587732/1‑150 (MQ=255)             gcggcgCTCCCAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCg                                                                                                                                     <  1:688876/150‑1 (MQ=255)                     cccAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGt                                                                                                                             <  2:527629/150‑1 (MQ=255)                      ccAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCTGCGAc                                                                                                                                                  <  2:759003/128‑1 (MQ=255)                      ccAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGtc                                                                                                                            >  1:573899/1‑150 (MQ=255)                      ccAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCACCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCATTGTATCACCGGTCTGCGta                                                                                                                            >  2:416105/1‑149 (MQ=255)                      ccAAGGGAAGCACCCACGGCTTCATAAGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCTGCGAc                                                                                                                                                  >  1:759003/1‑128 (MQ=255)                        aaGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGtctc                                                                                                                          <  1:937990/150‑1 (MQ=255)                        aaGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGtctc                                                                                                                          <  2:573899/150‑1 (MQ=255)                        aaGGGAAGCACCCACGGCTTCATAGGCACTGTACGGAATGGCGGCGTTGAACCCGGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGtctc                                                                                                                          >  1:139512/1‑150 (MQ=255)                           ggAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCgag                                                                                                                       >  1:484666/1‑150 (MQ=255)                           ggAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCgag                                                                                                                       >  2:972491/1‑150 (MQ=255)                            gAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATAAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCgaga                                                                                                                      <  2:143574/150‑1 (MQ=255)                            gAAGCACCCACGGCTTCATAGGCACGGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCgaga                                                                                                                      <  2:802627/150‑1 (MQ=255)                             aaGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCgaga                                                                                                                      >  2:774437/1‑149 (MQ=255)                             aaGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCgaga                                                                                                                      <  1:774437/149‑1 (MQ=255)                               gCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCa                                                                                                                                        >  2:844368/1‑129 (MQ=255)                               gCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCa                                                                                                                                        <  1:844368/129‑1 (MQ=255)                                cACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGaa                                                                                                                  >  1:1096011/1‑150 (MQ=255)                                 aCCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAac                                                                                                                 <  1:924485/150‑1 (MQ=255)                                             aTAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCAc                                                                                                    >  2:634470/1‑150 (MQ=255)                                               aGGCACTGGCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGa                                                                                                  <  2:139512/150‑1 (MQ=255)                                                     tGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTAttt                                                                                            >  1:928263/1‑150 (MQ=255)                                                         cGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTcc                                                                                        >  2:241354/1‑150 (MQ=255)                                                          ggAATGGCGCCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCt                                                                                       <  2:51351/150‑1 (MQ=255)                                                               ggcggcGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGt                                                                                  <  2:1048994/150‑1 (MQ=255)                                                                  ggcgTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGtctc                                                                                                                          <  2:927931/108‑1 (MQ=255)                                                                  ggcgTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGtctc                                                                                                                          >  1:927931/1‑108 (MQ=255)                                                                       tGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCtttt                                                                          <  1:972491/150‑1 (MQ=255)                                                                             cTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCttt                                                                    <  1:119237/150‑1 (MQ=255)                                                                             cTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTATCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACTGTTTGTCttt                                                                    >  1:465644/1‑150 (MQ=255)                                                                                                   cTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGCGAGGAACTCCCGACCTG‑TCAACACTGGCACCGGCATTGGTCTTTTGTCTTTGTAGAAGCCAACACACTTCGGa                                              <  1:300962/150‑1 (MQ=255)                                                                                                          cccTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTgcg                                       >  1:741518/1‑150 (MQ=255)                                                                                                             tGGCCAATCAGTTTCTTGATGGTTTCGGAGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTg                                    >  2:143654/1‑150 (MQ=255)                                                                                                                 cAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGggc                                >  2:440181/1‑150 (MQ=255)                                                                                                                  aaTCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGggcg                               >  1:113991/1‑150 (MQ=255)                                                                                                                   aTCAGTTTCTTGATGGTCTCGGGGGTCGCAGCGACACTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGggcgg                              >  2:470336/1‑150 (MQ=255)                                                                                                                    tCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGggcggc                             <  2:716789/150‑1 (MQ=255)                                                                                                                       gTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGggc                                >  2:283287/1‑144 (MQ=255)                                                                                                                       gTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGggc                                <  1:283287/144‑1 (MQ=255)                                                                                                                            ttGATGGTTTCGGGGGTCGAAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGTGGAAC‑TCCAAGGTGCACGGCTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGGGCGGCTGTGGAgc                     <  1:407980/150‑1 (MQ=255)                                                                                                                            ttGATGGATTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCGGGGGCGGCTGTGGAgc                     <  2:569985/150‑1 (MQ=255)                                                                                                                              gATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACAGGTCTGCGTCGCGAGGGGAAC‑ACCAATGTGAACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGGGCGGCTGTGGAgcgc                   >  2:157508/1‑150 (MQ=255)                                                                                                                                 ggTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGGGCGGCTGTGGAGCGCCaa                <  2:1096011/150‑1 (MQ=255)                                                                                                                                   tttCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCGACACACTTCGGATGGTGCGCTgg                                   <  2:797146/129‑1 (MQ=255)                                                                                                                                   tttCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCGACACACTTCGGATGGTGCGCTgg                                   >  1:797146/1‑129 (MQ=255)                                                                                                                                   tttCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGGGCGGCTGTGGAGCGCCAAcc              >  2:150464/1‑150 (MQ=255)                                                                                                                                     acaagggTCGCTGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAg                                                             <  1:878522/97‑1 (MQ=255)                                                                                                                                     acaagggTCGCTGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAg                                                             >  2:878522/5‑101 (MQ=255)                                                                                                                                     acaagggTCGCTGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGGGCGGCTGTGGAGCGCCAACccg            >  2:446810/5‑150 (MQ=255)                                                                                                                                     acaagggTCGCTGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGGGCGGCTGTGGAGCGCCAACccg            >  1:161112/5‑150 (MQ=255)                                                                                                                                               cAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGGGCGGCTGTGGAGCGCCAACCCGCCGGATCCcg  <  1:634470/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      CCGTAGGCATCGGCGGCGCTCCCAAGGGAAGCACCCACGGCTTCATAGGCACTGTCCGGAATGGCGGCGTTGAGCCCTGCCCCCCGTTGAACGGTGACCTGATGGCCCTGGCCAATCAGTTTCTTGATGGTTTCGGGGGTCGCAGCGACCCTTGTCTCACCGGTCTGCGTCTCGAGAGGAAC‑ACCAATGTGCACGACTATTTCTCCTGCGGTGACCTTTTGTCTTTGTAGAAGCCAACACACTTCGGATGGTGCGCTGGGGCGGCTGTGGAGCGCCAACCCGCCGGATCCCG  >  NC_002947/166713‑167004

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