Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	445,086	T→G	10.7%	intergenic (+59/‑229)	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,086	0	T	G	10.7%	118.8 / 8.0	47	intergenic (+59/‑229)	PP_0367/PP_0368	hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base G (2/3);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GCGGTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAG  >  NC_002947/444938‑445230                                                                                                                                                     |                                                                                                                                                  gcgggtggCAAGACCGAACTCGAGCGCCGCGCAAACGGCGTCAACGCATCCGCCGAAATGCGTGGCCGGATCATCGCCCCCTACGCCTGAATACTTGCGGCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAt                                                                                                                                                   <  6:379514/143‑1 (MQ=255)     tGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGGTCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTCTTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGt                                                                                                                                               <  1:76965/149‑1 (MQ=255)           aGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                         <  7:269339/149‑1 (MQ=255)               cGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTACGAACAGa                                                                                                                                      >  6:278511/1‑148 (MQ=255)                      gAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGCCACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATGCCCGCGGATGCGTCCGCTCAGCcaacaaat                                                                                                                              >  5:280021/1‑147 (MQ=255)                        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                                                                                                       >  7:121305/1‑146 (MQ=255)                                                aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                       <  8:121305/146‑1 (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)                                                                                       cTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                         >  3:20451/1‑73 (MQ=255)                                                                                       cTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                         <  4:20451/73‑1 (MQ=255)                                                                                         gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTACCCCGAAAACGGACCAAAAGACAAAAATTTTTTGTGGTAGGGGCGATAAGCGGG‑CAAAGCCCCGCCAACAAAGGTTTGTGgcatagcatgcaagat                                                          >  4:86640/1‑143 (MQ=255)                                                                                             aCTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                           <  1:144749/97‑1 (MQ=255)                                                                                             aCTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                           >  2:144749/1‑97 (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     >  2:243089/1‑149 (MQ=255)                                                                                                                                              ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt     >  1:278994/1‑149 (MQ=255)                                                                                                                                               cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACaa                                                                                   <  7:232935/70‑1 (MQ=255)                                                                                                                                               cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACaa                                                                                   >  8:232935/1‑70 (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)                                                                                                                                                     |                                                                                                                                                  GCGGTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAG  >  NC_002947/444938‑445230

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