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	91.4%	101.5 / 4.5	35	intergenic (+139/+76)	PP_4061/PP_4063	hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (12/20);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG  >  NC_002947/4585917‑4586198                                                                                                                                               |                                                                                                                                                 aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGcc                                                                                                                                           <  1:422371/148‑1 (MQ=39) aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCg                                                                                                                                          <  5:35619/149‑1 (MQ=39)        aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGCC‑C‑CTGCCGGCCTTg                                                                                                                                    >  3:180130/1‑148 (MQ=38)          cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAAt                                                                                                                                 >  1:19864/1‑149 (MQ=39)              gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATg‑cgg                                                                                                                            <  5:249008/149‑1 (MQ=39)                    gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGc                                                                                                                      <  7:381637/149‑1 (MQ=39)                       tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc                                                                                                                   <  2:398923/149‑1 (MQ=39)                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCt                                                                                                                  <  5:61565/148‑1 (MQ=39)                                      cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGc                                                                                                         <  8:98290/144‑1 (MQ=38)                                      cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       <  4:377528/146‑1 (MQ=38)                                      cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       <  7:124042/146‑1 (MQ=38)                                           ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTAcc                                                                                                <  2:19864/148‑1 (MQ=38)                                              cATTGAAAGGGTGAGGATGACCGGGGGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcg                                                                                              <  3:292539/147‑1 (MQ=37)                                                 tGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑a                                                                                                                                                   <  6:18495/93‑1 (MQ=37)                                                 tGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑a                                                                                                                                                   >  5:18495/1‑93 (MQ=37)                                                   aaaGGTTGAGGATGCCCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCAACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAaga                                                                                        <  6:186284/148‑1 (MQ=37)                                                   aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       <  6:241750/133‑1 (MQ=38)                                                   aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       >  5:241750/1‑133 (MQ=38)                                                       gTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       >  1:127766/1‑129 (MQ=38)                                                       gTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCgg                                                                                                       <  2:127766/129‑1 (MQ=38)                                                        ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGa                                                                                                           >  4:397859/1‑124 (MQ=38)                                                        ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGa                                                                                                           <  3:397859/124‑1 (MQ=38)                                                                gACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCggg                                                                                                      >  5:245503/1‑121 (MQ=37)                                                                gACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCggg                                                                                                      <  6:245503/121‑1 (MQ=37)                                                                  ccGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATg‑cg                                                                                                                             >  1:224086/1‑96 (MQ=25)                                                                  ccGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATg‑cg                                                                                                                             <  2:224086/96‑1 (MQ=25)                                                                             gCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGttt                                                                                                   <  1:1034/111‑1 (MQ=25)                                                                             gCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGttt                                                                                                   >  2:1034/1‑111 (MQ=25)                                                                              cAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAAt                                                                                                                                 <  4:104149/81‑1 (MQ=21)                                                                              cAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAAt                                                                                                                                 >  3:104149/1‑81 (MQ=21)                                                                                                                  tcgcgggTAAACCCGCTCCTACCAGGG‑G‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGc                                       <  2:273131/129‑1 (MQ=14)                                                                                                                  tcgcgggTAAACCCGCTCCTACCAGGG‑G‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGc                                       >  1:273131/6‑134 (MQ=14)                                                                                                                  tcgcgggTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACg                        >  4:251305/6‑149 (MQ=255)                                                                                                                     cgggTAAACCCGCTCCTACCAGGG‑ACC‑GGGCCGGCCTTGAATG‑CGGCGCTGCCCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCAt                     <  3:251305/147‑1 (MQ=255)                                                                                                                        gTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg                  >  3:180645/1‑149 (MQ=255)                                                                                                                        gTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg                  >  7:74834/1‑149 (MQ=255)                                                                                                                                        cAGGG‑ACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTg  >  4:186570/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                 AGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG  >  NC_002947/4585917‑4586198

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