Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	445,115	G→A	12.4%	intergenic (+88/‑200)	PP_0367 → / → PP_0368	hypothetical protein/acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	445,115	0	G	A	12.4%	124.7 / 7.8	49	intergenic (+88/‑200)	PP_0367/PP_0368	hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (25/17);  new base A (3/3);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

GGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCG  >  NC_002947/444973‑445262                                                                                                                                               |                                                                                                                                                    ggTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                               <  5:194308/149‑1 (MQ=255)           gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                     <  6:83789/149‑1 (MQ=18)                cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGTATGCGTCCGATCAGACAACAATTTTGTCTGGt                                                                                                                                                    >  8:25324/1‑129 (MQ=255)                cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGTATGCGTCCGATCAGACAACAATTTTGTCTGGt                                                                                                                                                    <  7:25324/129‑1 (MQ=255)                  acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGATTCCCGCGAATGCGGCCGAAGAGCCAAAAAATTTTTATGGTGGGGGCTCGTAGCGGGGCAGcc                                                                                                                              >  5:226943/1‑149 (MQ=255)                          cTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                     <  1:283880/134‑1 (MQ=17)                          cTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                     >  2:283880/1‑134 (MQ=17)                                  gATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGCATGCGTCCGATCAGACAACAATTTTGTGTGGTCGGGCCTCTTCGCGGGCAAGCCCGCGCCTAAAAggggt                                                                                                              >  3:322540/1‑147 (MQ=255)                                    tcatcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATAAGACAACAATTTTTTCTGGTCGGGCCTCTTAGCGGGCAATCCCGCTCCGACAACGGTTtg                                                                                                            >  6:101027/1‑149 (MQ=255)                                          acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGTTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcat                                                                                                      >  5:88123/1‑149 (MQ=255)                                                        aTACTTGCGTCATCCCTTATTTCTGCGCCTGCCTTGTTGGTGGGGGCTTGCCGGCTATTGCGTCCTATCACAAAAAATTTTTGTCTGGTCGGGCCTCTTCGCGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                        <  2:341892/149‑1 (MQ=255)                                                            ttGCGTCATCCACTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATGCCCGCGAATGAGTCCCATAACACAACAAAATTTTCTGGTGGGGCCTCGTACCGGGCAAGCCAGGGCCTACAAGCGTTTGAGGCATTGCATGCCAGCTTTggc                                                                                     >  2:263719/1‑146 (MQ=255)                                                                   aTCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGTTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCg                                                                             <  6:88123/149‑1 (MQ=255)                                                                     ccccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATTCGTCCCATAAGACAACAATTTTTTATGGTCGGGCCACTTAGCGGGCCAGCCCCCTCCGAAAAGGGTTTTTGGCATGGCAAGTCTGATTTTGCATTTTTAGTg                                                                           >  2:32186/1‑149 (MQ=255)                                                                             tCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcg                                                                                                                                        >  4:278652/1‑80 (MQ=255)                                                                             tCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcg                                                                                                                                        <  3:278652/80‑1 (MQ=255)                                                                                  ccctgcctTGTTGCTTCGGGCTTGCCCGCGAATGCGTCCGATCATCAAACAATTTTGGCTGGTCGGGACTCTTCGCGGGCATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                               <  1:263719/147‑1 (MQ=255)                                                                                                cgCGGGCTGCCCCGCTAAGGCGTCCGATCCCACAAAAATTTTGTCTGGCGGGGCCTCTCCGCGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                                <  6:226943/148‑1 (MQ=255)                                                                                                  cGGGTTTGCCCGCTAATGCGGCCGATCACAAAACAATTTTGTCGGGGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                              <  8:27770/149‑1 (MQ=255)                                                                                                         gCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                                        >  6:347114/1‑132 (MQ=255)                                                                                                         gCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                                        <  5:347114/132‑1 (MQ=255)                                                                                                          cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCTTAAATGTCACCATAGTCCATCCTTGACAGAAAGc                                      >  2:174557/1‑149 (MQ=255)                                                                                                          cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAg                                       <  1:47634/148‑1 (MQ=255)                                                                                                          cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                                      >  5:39090/1‑149 (MQ=18)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCt                                                 >  6:327631/1‑137 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCt                                                 <  5:327631/137‑1 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                              <  7:350400/140‑1 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                              >  8:350400/1‑140 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                     >  4:21296/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                     >  4:178987/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                     >  5:275895/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                     >  5:134588/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                     >  1:229380/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTACTACAAGGGTTTGTGGCATGGCATGACTGATTTGCAATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAgaa                                     >  4:230645/1‑147 (MQ=255)                                                                                                             gctaaTGCGTCCGATCAGACAACAATTTTGTCTGGGCGGGCCTATTCGCGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                   <  1:221988/146‑1 (MQ=255)                                                                                                             gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCa                                                     <  5:423140/131‑1 (MQ=255)                                                                                                             gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCa                                                     >  6:423140/1‑131 (MQ=255)                                                                                                             gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                     >  2:106113/1‑147 (MQ=255)                                                                                                                            cacacaaaaTTTTTGTTTGGCCGGGCCTCTCCGGGGGAATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGttt                    <  2:365036/146‑1 (MQ=255)                                                                                                                              gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc                   >  1:403752/1‑148 (MQ=255)                                                                                                                                caacaaTTTTTTCTGTTCGGGCCTCTTCGCGTGCAAGGCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                <  3:230645/149‑1 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                >  2:33581/1‑148 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                >  2:27288/1‑148 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGAATTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                >  1:101088/1‑148 (MQ=255)                                                                                                                                   caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAggg                                                                                                                <  5:379266/50‑1 (MQ=255)                                                                                                                                   caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAggg                                                                                                                >  6:379266/1‑50 (MQ=255)                                                                                                                                   caaTTTTGTCTGGTCGGGCCGCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa             <  5:101027/149‑1 (MQ=255)                                                                                                                                   aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                              <  4:179480/111‑1 (MQ=25)                                                                                                                                   aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                              >  3:179480/6‑116 (MQ=25)                                                                                                                                        ttgtcttgtcGGGCATCTTCGCTGTCAATCCCTCTCCTACAAGTGTTTGTGGCATTTCATGACTGATTTGCCATTTTTCGTGGCGTAAATTGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt         <  1:174557/148‑1 (MQ=255)                                                                                                                                              tgTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg  <  1:296435/148‑1 (MQ=255)                                                                                                                                               gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGNTAAACGTAAGTTTCAAACAACTGTTTGATc   <  1:32186/147‑1 (MQ=255)                                                                                                                                               gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAGTTGTTTGaaa   <  6:134588/147‑3 (MQ=255)                                                                                                                                               gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt    <  4:322540/146‑1 (MQ=255)                                                                                                                                               gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATc   <  2:229380/147‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    GGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCG  >  NC_002947/444973‑445262

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑