Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,173,244	T→C	10.6%	intergenic (‑285/+47)	fadD‑I ← / ← fadD‑II	long‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,173,244	0	T	C	10.6%	123.0 / 5.8	47	intergenic (‑285/+47)	fadD‑I/fadD‑II	long‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base C (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 = 8.81e-01

AGGCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGT  >  NC_002947/5173098‑5173389                                                                                                                                                    |                                                                                                                                                  agGCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTg                                                                                                                                                 <  2:52837/149‑1 (MQ=255)  ggCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTg                                                                                                                                                 <  7:30010/148‑1 (MQ=255)    ccAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTggg                                                                                                                                               <  7:126463/148‑1 (MQ=255)     cAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTggg                                                                                                                                               >  8:66363/1‑147 (MQ=255)     cAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACACTgccc                                                                                                                                              >  7:54834/1‑145 (MQ=255)      aaaaCAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCCGGGTAAACCCGCTCCTACATTGGGCa                                                                                                                                             <  1:3562/149‑1 (MQ=255)                ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCg                                                                                                                                  >  6:16797/1‑149 (MQ=255)                  aTGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCgg                                                                                                                                 <  4:151986/148‑1 (MQ=255)                   tGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCggg                                                                                                                                <  8:54834/148‑1 (MQ=255)                                 aCTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCgaagaa                                                                                                                 <  7:66363/149‑1 (MQ=255)                                  cTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCGAAAATGTCGCAGATATGGCTGTATAGAGCCTGTTAGAGCTGTTTAAGGCATTCG‑GGGGTAAACCGGGTCCTACATTGGGCAGTTTAGGAGCGGGTTTACCCGCgaagaa                                                                                                                 <  6:111871/148‑1 (MQ=255)                                               aCAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcgc                                                                                                                       <  3:153654/129‑1 (MQ=255)                                               aCAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcgc                                                                                                                       >  4:153654/1‑129 (MQ=255)                                                    tCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGc                                                                                               >  1:251492/1‑148 (MQ=255)                                                                tAATTATCCAAAATGTCGCAGAT‑AGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACt                                                                                 >  2:338525/1‑149 (MQ=255)                                                                     aTCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCccg                                                                             >  4:104956/1‑149 (MQ=255)                                                                         aaaaTGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCgcgc                                                                                     >  8:201861/1‑137 (MQ=255)                                                                         aaaaTGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCgcgc                                                                                     <  7:201861/137‑1 (MQ=255)                                                                          aaaTGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCggg                                                                                                                                <  4:26983/93‑1 (MQ=25)                                                                                 gCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCtt                                                                 >  3:95648/1‑149 (MQ=255)                                                                                            tGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCg                                                      >  7:156291/1‑149 (MQ=255)                                                                                             gTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGt                                                     >  3:163885/1‑149 (MQ=255)                                                                                                  gagCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATgg                                                <  4:162640/149‑1 (MQ=255)                                                                                                  gagCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATg                                                 <  3:104956/148‑1 (MQ=255)                                                                                                   agCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATggg                                               >  5:298142/1‑149 (MQ=255)                                                                                                    gCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGc                                              <  6:377767/149‑1 (MQ=255)                                                                                                       tGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAAc                                           >  3:316847/1‑149 (MQ=255)                                                                                                               cGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATcgc                                                                                       <  4:275749/97‑1 (MQ=255)                                                                                                               cGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATcgc                                                                                       >  3:275749/1‑97 (MQ=255)                                                                                                                 gTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGt                                                                                                      >  6:183334/1‑80 (MQ=255)                                                                                                                 gTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGt                                                                                                      <  5:183334/80‑1 (MQ=255)                                                                                                                        gCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGAt                          <  4:316847/149‑1 (MQ=255)                                                                                                                                    tAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAATCCATCACGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTTCCCACCGGAGTCATCGGCAACTATTCGCGCAGAACGAATTGCTTGGGTa               >  3:13242/1‑149 (MQ=38)                                                                                                                                        cccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACttt           >  6:256134/1‑149 (MQ=255)                                                                                                                                         ccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg          <  4:163885/149‑1 (MQ=255)                                                                                                                                           gCTCCTACATTGGGAAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTa        <  5:16797/149‑1 (MQ=255)                                                                                                                                           gCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt            >  8:344249/1‑145 (MQ=255)                                                                                                                                           gCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt            <  7:344249/145‑1 (MQ=255)                                                                                                                                            cTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTg                                                                <  6:368564/92‑1 (MQ=255)                                                                                                                                            cTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTg                                                                >  5:368564/1‑92 (MQ=255)                                                                                                                                             tCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCtata                                                                                           <  6:220062/64‑1 (MQ=255)                                                                                                                                             tCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCtata                                                                                           >  5:220062/1‑64 (MQ=255)                                                                                                                                               cTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCg    >  7:196356/1‑149 (MQ=255)                                                                                                                                                tACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                    <  5:134345/132‑1 (MQ=255)                                                                                                                                                tACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCtt                    >  6:134345/1‑132 (MQ=255)                                                                                                                                                 aCATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGt  >  1:161810/1‑149 (MQ=255)                                                                                                                                                  cATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCgg   <  4:95648/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  AGGCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCG‑CGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGT  >  NC_002947/5173098‑5173389

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