Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,586,057	+C	intergenic (+140/+75)	PP_4061 → / ← PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,586,056	1	.	C	88.9%	104.6 / 8.0	36	intergenic (+139/+76)	PP_4061/PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (17/15);  total (19/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 = 1.00e+00

AGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCC‑TACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG  >  NC_002947/4585910‑4586198                                                                                                                                                      |                                                                                                                                                 aGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGa                                                                                                                                                   <  4:234495/148‑1 (MQ=255)     ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGA‑ccgtg                                                                                                                                             >  7:124707/1‑145 (MQ=39)     ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGA‑ccgtg                                                                                                                                             >  7:105173/1‑145 (MQ=39)                ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGaa                                                                                                                                  <  5:22450/149‑1 (MQ=39)                 cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCACCAAACCCGGCC‑C‑GTGCCGGCCATGCAt                                                                                                                                 >  6:278267/1‑149 (MQ=255)                            cGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcc                                                                                                                     <  1:322949/149‑1 (MQ=39)                                  cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc                                                                                                                   >  1:301980/1‑145 (MQ=38)                                  cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc                                                                                                                   <  2:301980/145‑1 (MQ=38)                                   aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTa                                                                                                              <  8:124707/149‑1 (MQ=39)                                           tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTgg                                                                                                                <  1:3752/139‑1 (MQ=38)                                           tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTgg                                                                                                                >  2:3752/1‑139 (MQ=38)                                             cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCgg                                                                                                       <  5:161656/146‑1 (MQ=14)                                             cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       <  8:105173/146‑1 (MQ=38)                                             cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       <  4:255475/146‑1 (MQ=38)                                                  ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTAccc                                                                                               <  3:260050/149‑1 (MQ=14)                                                  ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       >  8:131944/1‑141 (MQ=38)                                                  ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       <  7:131944/141‑1 (MQ=38)                                                   tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTAccc                                                                                               >  4:260050/1‑148 (MQ=14)                                                     cATTGAAAGCTTGAGGATGACCGGTTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TCCCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgc                                                                                             <  5:36951/148‑1 (MQ=37)                                                       ttGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGc                                                                                                                      >  8:140657/1‑121 (MQ=38)                                                       ttGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGc                                                                                                                      >  8:145662/1‑121 (MQ=38)                                                       ttGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGc                                                                                                                      <  7:145662/121‑1 (MQ=38)                                                       ttGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGc                                                                                                                      <  7:140657/121‑1 (MQ=38)                                                         gAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCt                                                                                                                                      >  7:280116/1‑104 (MQ=37)                                                         gAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCt                                                                                                                                      <  8:280116/104‑1 (MQ=37)                                                          aaaGATTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTg                                                                                                                       <  4:176110/117‑1 (MQ=37)                                                          aaaGATTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTg                                                                                                                       >  3:176110/1‑117 (MQ=37)                                                                    gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgc                                                                                             <  8:338691/133‑1 (MQ=37)                                                                    gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgc                                                                                             >  7:338691/1‑133 (MQ=37)                                                                            ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAAt                                                                                                                                 <  6:272343/90‑1 (MQ=25)                                                                            ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAAt                                                                                                                                 >  5:272343/1‑90 (MQ=25)                                                                              tGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTAGCGGGTAAACCCGCTAC‑TACCAGGGACC‑GTGGCGGCCGTTAATG‑CGGCGCTGGCCCTGGTAGGGGCGGGGTTTACCGCGAAGCGGCCGGCCGCGCCAGCACAta                                                                   >  6:367768/1‑148 (MQ=255)                                                                                                                                    cccGCTCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt             >  2:92823/1‑149 (MQ=255)                                                                                                                                     ccGCTCC‑TACCAGCGACC‑GTGCCGCCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCTGGCCGCGCCAGCACATCAACCCCCAGATGCAGCCGCAATCTCCTCCACACTGAACTCACGCATGCGGAAca             >  4:251496/1‑147 (MQ=255)                                                                                                                                         tCC‑TACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCGGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTg         >  6:221440/1‑148 (MQ=255)                                                                                                                                             tACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc     >  4:195450/1‑149 (MQ=255)                                                                                                                                              aCCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCt    >  1:288130/1‑149 (MQ=255)                                                                                                                                              aCCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCACCACACTGATCTCACGCATGCGGAACTTCTGCACCt    >  2:96583/1‑149 (MQ=255)                                                                                                                                                cAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTg  >  2:371811/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                 AGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCC‑TACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG  >  NC_002947/4585910‑4586198

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑