Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	445,088	C→G	10.7%	intergenic (+61/‑227)	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,088	0	C	G	10.7%	118.8 / 7.6	47	intergenic (+61/‑227)	PP_0367/PP_0368	hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (20/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 = 2.16e-01

TGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAG  >  NC_002947/444942‑445230                                                                                                                                                   |                                                                                                                                                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                                                                                                                                         <  4:20451/73‑1 (MQ=255)                                                                                   cTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                         >  3:20451/1‑73 (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     >  2:243089/1‑149 (MQ=255)                                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt     >  1:278994/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  <  1:229941/147‑1 (MQ=255)                                                                                                                                               aaTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg  >  2:229941/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                TGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAG  >  NC_002947/444942‑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: ‑