Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	445,120	G→A	10.9%	intergenic (+93/‑195)	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,120	0	G	A	10.9%	148.2 / 7.1	55	intergenic (+93/‑195)	PP_0367/PP_0368	hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base A (3/3);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

GGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTC  >  NC_002947/444973‑445267                                                                                                                                                    |                                                                                                                                                    ggTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                                    <  5:194308/149‑1 (MQ=255)           gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                          <  6:83789/149‑1 (MQ=18)                  acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGATTCCCGCGAATGCGGCCGAAGAGCCAAAAAATTTTTATGGTGGGGGCTCGTAGCGGGGCAGcc                                                                                                                                   >  5:226943/1‑149 (MQ=255)                          cTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                          >  2:283880/1‑134 (MQ=17)                          cTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                          <  1:283880/134‑1 (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                                                             <  5:347114/132‑1 (MQ=255)                                                                                                         gCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                                             >  6:347114/1‑132 (MQ=255)                                                                                                          cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCTTAAATGTCACCATAGTCCATCCTTGACAGAAAGc                                           >  2:174557/1‑149 (MQ=255)                                                                                                          cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAg                                            <  1:47634/148‑1 (MQ=255)                                                                                                          cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                                           >  5:39090/1‑149 (MQ=18)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCt                                                      <  5:327631/137‑1 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCt                                                      >  6:327631/1‑137 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                                   >  8:350400/1‑140 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                                   <  7:350400/140‑1 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                          >  4:178987/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                          >  4:21296/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                          >  1:229380/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                          >  5:134588/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                          >  5:275895/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:27288/1‑148 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                     >  2:33581/1‑148 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGAATTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                     >  1:101088/1‑148 (MQ=255)                                                                                                                                   caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAggg                                                                                                                     >  6:379266/1‑50 (MQ=255)                                                                                                                                   caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAggg                                                                                                                     <  5:379266/50‑1 (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)                                                                                                                                                tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg       >  2:296435/1‑147 (MQ=255)                                                                                                                                                tCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACGGATTTTCCATTTTTAGTGGCGTAAATGGCACCATTGTGCATACATGACAGAAAGCAGAGCTAAAACGAAATTTCAAACAACAGTTTGATAGGc     >  1:318680/1‑149 (MQ=255)                                                                                                                                                   ggCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTc  <  3:21296/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGCTTC  >  NC_002947/444973‑445267

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