Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	445,115	G→A	11.9%	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	11.9%	106.0 / 5.5	42	intergenic (+88/‑200)	PP_0367/PP_0368	hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (22/15);  new base A (3/2);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

GGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATC  >  NC_002947/444973‑445261                                                                                                                                               |                                                                                                                                                    gggcgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                               <  7:247930/146‑1 (MQ=255) gctcggcAAGGCATCCGACGAAATGCTTCGCCGCATCATCGACACCTACGCCTGAATACTTGCGGCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                               <  8:298974/147‑1 (MQ=255)             aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGc                                                                                                                                   <  6:333242/149‑1 (MQ=255)                 gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc                                                                                                                                                   >  4:37213/1‑129 (MQ=255)                 gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc                                                                                                                                                   <  3:37213/129‑1 (MQ=255)                                    tcatcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCCCGCATTCGTACGCTAAAACAACAATTTTTTAAGGTCGGGCCTCTTAGCGGGGC‑AGCCCGCGCCAAACAGGGGTtg                                                                                                           >  7:56890/1‑149 (MQ=255)                                                                      cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGCATGCGTCCGCTCAGACAACAATTTTTTGTGGTGGGGCCTAATAGCCGGCA‑AGCCCGCGCCGACAAGGGGTTGTGGGATGGCATGAGTGATTTTGCATTTTTAGTgg                                                                         >  6:155680/1‑149 (MQ=255)                                                                      cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTc                                                                                                                                         <  1:70690/86‑1 (MQ=255)                                                                      cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTc                                                                                                                                         >  2:70690/1‑86 (MQ=255)                                                                         accttCAGCGCCTGCCTTGTCGGCGCGGGCTTGCCCGCGAATGCGTCCTACCACAAAACAATTTTGTCTGGCCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGt                                                                      <  7:334441/146‑1 (MQ=38)                                                                                               gAGCGGGCTTGCCCGCTAATGCGTCCGATCAGAAAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                                       <  1:228687/142‑1 (MQ=255)                                                                                                     gCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCGATCCTGGa                                             >  7:221795/1‑146 (MQ=255)                                                                                                          cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                                     >  7:150484/1‑149 (MQ=255)                                                                                                           gcgcgAATGCGTACGATCAGACAACAATTTTGTATGGTCGGGCCTCTTCGCGTGCA‑AGCGCGCTCCTTCAAGGGTTTGTTGCTTGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGTCTCCTTAGTCCATCCTTGa                                             <  4:1761/139‑1 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAGAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                                     >  1:98631/1‑149 (MQ=255)                                                                                                           ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGTCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGACCATCATTGACAGATAgg                                     >  3:1761/1‑147 (MQ=255)                                                                                                            cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                    >  5:28086/1‑148 (MQ=255)                                                                                                                    cGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAc                           >  8:38742/1‑149 (MQ=255)                                                                                                                     gcccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                       <  5:252000/86‑1 (MQ=18)                                                                                                                     gcacGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                       >  6:252000/4‑88 (MQ=18)                                                                                                                     gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACg                          >  8:71125/1‑149 (MQ=255)                                                                                                                       ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg                                 >  2:46324/1‑140 (MQ=255)                                                                                                                       ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg                                 <  1:46324/140‑1 (MQ=255)                                                                                                                       ccGACCAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGTCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                         >  6:106845/1‑148 (MQ=38)                                                                                                                         gATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAg                      >  4:275588/1‑149 (MQ=255)                                                                                                                         gATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATCGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAg                      >  6:364900/1‑149 (MQ=255)                                                                                                                               acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAt                                                          >  2:295711/1‑107 (MQ=255)                                                                                                                               acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa                 >  1:257373/1‑148 (MQ=255)                                                                                                                               acaacaATTTTCTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAt                                                          <  1:295711/107‑1 (MQ=255)                                                                                                                               acaacaATTTTCCCGTGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGCCATGGCATGACTGATTTCCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                <  8:221795/149‑1 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                  >  7:97998/1‑129 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                  <  8:97998/129‑1 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                       <  6:354344/140‑1 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa                       >  5:354344/1‑140 (MQ=255)                                                                                                                                 aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa               >  3:2099/1‑148 (MQ=255)                                                                                                                                    aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac           >  6:129744/1‑149 (MQ=255)                                                                                                                                    aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac           <  5:155680/149‑1 (MQ=255)                                                                                                                                     ttgttGTCTGATCGGACGCATTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAATCAACt          >  7:338777/4‑149 (MQ=38)                                                                                                                                     aTTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGCTGAGTGGCTTGGCATGACTGATTTCCCATTTTTCGTGGCGTAAATGGCACCAAAGTCCATCCTTGACAGAAACGTGAGCTAAAAGTAAGTTTCAAACAACt          >  4:314428/1‑149 (MQ=255)                                                                                                                                       gttGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg         <  5:106845/147‑1 (MQ=255)                                                                                                                                         tGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt       <  2:257373/148‑1 (MQ=255)                                                                                                                                           tCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGa    >  8:103607/1‑149 (MQ=255)                                                                                                                                               gTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg                                                                          <  7:231469/75‑1 (MQ=255)                                                                                                                                               gTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg                                                                          >  8:231469/1‑75 (MQ=255)                                                                                                                                               gTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATc  <  4:2099/147‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    GGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCA‑AGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATC  >  NC_002947/444973‑445261

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