Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	445,107	T→G	13.8%	intergenic (+80/‑208)	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,107	0	T	G	13.8%	129.7 / 14.0	51	intergenic (+80/‑208)	PP_0367/PP_0368	hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (25/19);  new base G (4/3);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

GCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAAC  >  NC_002947/444961‑445252                                                                                                                                                   |                                                                                                                                                   gCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                  <  6:151498/149‑1 (MQ=255) gCTCGGCGCAAAGGTCGTCAACGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                  <  6:15610/149‑1 (MQ=255)     ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                              <  8:207338/149‑1 (MQ=255)     ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCt                                                                                                                                              >  6:302220/1‑149 (MQ=37)       cgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGGCCCCGAATGAGTCCGCTCAGACAACAAATTTGTGTgg                                                                                                                                            >  2:279571/1‑149 (MQ=255)       cgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCAGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTgg                                                                                                                                            <  5:183617/149‑1 (MQ=255)        gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTgat                                                                                                                                           <  7:60047/149‑3 (MQ=25)           aaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                        <  7:292475/149‑1 (MQ=255)                cgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctc                                                                                                                                   <  5:131458/149‑1 (MQ=255)                 gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctct                                                                                                                                  <  1:334800/149‑1 (MQ=255)                      ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGGCCGCGAATGCGTCCCATAAGAAAACAATTTTTTATGGTGCGGCCTCTTagcg                                                                                                                             >  3:301257/1‑149 (MQ=255)                       gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAATATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                            <  1:168365/149‑1 (MQ=18)                       gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                            <  5:302220/149‑1 (MQ=18)                         aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                                             <  8:121305/146‑1 (MQ=255)                         aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                                             >  7:121305/1‑146 (MQ=255)                           ccGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATCAGACAACAATTTTGTGTGGTCGGGCCTCTTAGCGGG‑Caa                                                                                                                       >  5:375033/1‑149 (MQ=255)                                gaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCg                                                                                                                  <  2:332510/149‑1 (MQ=255)                                  aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATTCCCCCGGAATCGTGCCCTCACACAACAAATTTTTAAGGTGGGGCCACTTAGCCGG‑CCAGcccccg                                                                                                                >  7:2557/1‑146 (MQ=255)                                     gCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCCCGAATTCGTCCGATAAGACAAAAATTTTTTCTGGGCGGGCCGATTAGCGGGCCAAGccgcg                                                                                                                 >  8:152349/1‑143 (MQ=255)                                           ccGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGCCTGCCTTGTCGGCGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAgg                                                                                                       <  1:258129/149‑1 (MQ=255)                                                                  gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTACCCCGAAAACGGACCAAAAGACAAAAATTTTTTGTGGTAGGGGCGATAAGCGGG‑CAAAGCCCCGCCAACAAAGGTTTGTGgcatagcatgcaagat                                                                                >  4:86640/1‑143 (MQ=255)                                                                      aCTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                 >  2:144749/1‑97 (MQ=255)                                                                      aCTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                 <  1:144749/97‑1 (MQ=255)                                                                               aTCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATCAGACAACAATTTTTTCTGGTCGGGCCTCTTAGCGGG‑CAAGCCCGCTCCGACAACGGTTTGTGGCATGGCATGACTGATTTGCCAtttttag                                                                   >  7:93497/1‑147 (MQ=255)                                                                                tCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGCGTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CTAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt                                                                  <  5:278511/149‑1 (MQ=255)                                                                                  cccTATTTCCGCCCCGGCCTTGTAGGAGCGGGCTTGCCCGCTAATGCGTCCTATCACACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg                                                                 <  4:50832/148‑1 (MQ=38)                                                                                              gcctgcctTGTAGGAGCGGGCTGGCCCGCTAATGCGTCCTATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCAc                                                    <  1:279571/149‑1 (MQ=255)                                                                                                          ttcgCGGGCAAGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                        >  7:246133/4‑149 (MQ=255)                                                                                                              cGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                    <  7:251288/149‑1 (MQ=255)                                                                                                                 gCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                    >  7:89959/1‑146 (MQ=255)                                                                                                                      cccGCTAATGCGTCCGTTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                            <  3:86640/149‑1 (MQ=255)                                                                                                                      cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                            >  5:212092/1‑149 (MQ=18)                                                                                                                       ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGAc                                                                                    >  2:4866/1‑92 (MQ=255)                                                                                                                       ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                           >  1:278994/1‑149 (MQ=255)                                                                                                                       ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                           >  2:243089/1‑149 (MQ=255)                                                                                                                        cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACaa                                                                                                         >  8:232935/1‑70 (MQ=255)                                                                                                                        cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACaa                                                                                                         <  7:232935/70‑1 (MQ=255)                                                                                                                        cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                           >  2:13328/1‑148 (MQ=255)                                                                                                                         gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                           >  4:330305/1‑147 (MQ=255)                                                                                                                            aaTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg                        >  2:229941/1‑147 (MQ=255)                                                                                                                            aaTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg                        <  1:229941/147‑1 (MQ=255)                                                                                                                                 gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACg                 >  5:159226/1‑149 (MQ=255)                                                                                                                                 gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACg                 >  2:172231/1‑149 (MQ=255)                                                                                                                                  tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACg                 >  3:320823/1‑148 (MQ=255)                                                                                                                                   ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                >  5:316564/1‑148 (MQ=38)                                                                                                                                        cAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGttt          >  2:61283/1‑149 (MQ=37)                                                                                                                                           acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa        >  3:152137/1‑148 (MQ=255)                                                                                                                                            caacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa      <  6:159226/149‑1 (MQ=255)                                                                                                                                             aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCGCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACACAAAGCTGAGATAAACTTAACTTTCaaa      >  8:288766/1‑148 (MQ=255)                                                                                                                                             aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc         >  5:21619/1‑145 (MQ=255)                                                                                                                                             aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc         <  6:21619/145‑1 (MQ=255)                                                                                                                                             aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa      >  3:134284/1‑148 (MQ=255)                                                                                                                                               catttttGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa   <  6:280021/146‑1 (MQ=255)                                                                                                                                                aaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac  <  8:93497/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAAC  >  NC_002947/444961‑445252

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