Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,068,152	A→T	16.1%	intergenic (‑97/+24)	plsC ← / ← PP_1845	1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,068,152	0	A	T	16.1%	124.7 / 12.7	54	intergenic (‑97/+24)	plsC/PP_1845	1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (26/19);  new base T (3/6);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCG  >  NC_002947/2068005‑2068294                                                                                                                                                     |                                                                                                                                                aaCGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAg                                                                                                                                                >  6:302660/1‑149 (MQ=255)       aaGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGATGGGTCAAAACTTGGTGCAGGCCTGTTACGGCCC‑TTTCTCGGGCAAGCAcgc                                                                                                                                           >  7:176413/1‑148 (MQ=255)        aGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTCCACTGGCTGGGGCAAACCTAGGTGAAGGCCTATTACGGACA‑TTTCGCGGGCAACCCagcg                                                                                                                                          >  8:54867/1‑148 (MQ=255)            aaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTTTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCTTGCCCGCGaaagg                                                                                                                                     <  2:245328/149‑3 (MQ=255)            aaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGaaag                                                                                                                                      <  4:318619/149‑2 (MQ=255)                cAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGAAAgggc                                                                                                                                   <  8:51130/148‑2 (MQ=255)                cAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGAAAgggc                                                                                                                                   <  8:156494/148‑2 (MQ=255)                      cagcagGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCG‑AAAGGGccgtaaca                                                                                                                           <  8:135953/149‑7 (MQ=38)                           ggAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAg                                                                                                                     >  7:309718/1‑149 (MQ=255)                            gAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGAGGGCAAGCCCGCGAAAACGGCCCGTGCGGC‑TAGc                                                                                                                    >  2:103262/1‑149 (MQ=255)                                  aCGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAatat                                                                                                              >  4:46146/1‑149 (MQ=255)                                  aCGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAatat                                                                                                              >  4:27619/1‑149 (MQ=255)                                      aaTGAATAATGCATGGTAGTCCTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTGGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATAtgctg                                                                                                          <  2:126161/149‑1 (MQ=255)                                      aaTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCcgc                                                                                                                                           >  5:78709/1‑117 (MQ=255)                                      aaTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCcgc                                                                                                                                           <  6:78709/117‑1 (MQ=255)                                      aaTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATAtgctg                                                                                                          <  5:302660/149‑1 (MQ=255)                                       aTGAAAAAAGCATGGTACGACTCAAACAAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTCGCTGGGGCAAACCTTGGTGAAGGCCTGTTACGGCCC‑TTTCGCGGGAAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATAtgctgc                                                                                                         >  4:320707/1‑149 (MQ=255)                                           aTAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGAt                                                                                                     <  5:193275/149‑1 (MQ=255)                                                 cATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCt                                                                                               >  3:345663/1‑149 (MQ=255)                                                     ggaCGACTCACACCAGACGGTCGCGGGTCAGTGTACGGGTGTGCACTGGCTTGGGCAAACCTTGGTGCAGGCCTTTTACGGCCA‑TTTCGCGGGCAAGCCCGCTAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGAt                                                                                           <  8:176413/147‑1 (MQ=255)                                                           cTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGAGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCtt                                                                                     >  7:150842/1‑149 (MQ=255)                                                                acCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGaa                                                                                                                                        >  4:120496/1‑94 (MQ=255)                                                                acCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGaa                                                                                                                                        <  3:120496/94‑1 (MQ=255)                                                                        gTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTc                                                                         >  5:339843/1‑148 (MQ=255)                                                                                  attttACGGGTGTGGAATTGCTGGGGCAAACCTTGGTGGATTCCTGTTACGGCCC‑TTTCGGCGGCAATGCCGCGAAAATGGCCCGTGCGGCTTTGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCt                                                               <  1:103262/145‑1 (MQ=255)                                                                                      tACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgc                                                           <  4:345663/148‑1 (MQ=255)                                                                                      tACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGACCTTTgag                                                          >  4:100069/1‑147 (MQ=255)                                                                                             gtgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTc                                                   <  3:27619/149‑1 (MQ=255)                                                                                                     gctggGGCAAACCTTGGAGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                             <  8:150842/147‑1 (MQ=255)                                                                                                      ctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGgtg                                          <  3:46146/149‑1 (MQ=255)                                                                                                       tggGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                             <  4:51553/145‑1 (MQ=255)                                                                                                       tggGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                             >  3:51553/1‑145 (MQ=255)                                                                                                         gggCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGtt                                        <  8:239243/148‑1 (MQ=255)                                                                                                           gCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCtt                                     >  5:296867/1‑149 (MQ=255)                                                                                                               aCCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGcc                                                                                       >  7:352480/1‑95 (MQ=255)                                                                                                               aCCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGcc                                                                                       <  8:352480/95‑1 (MQ=255)                                                                                                                 cTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTgcgc                               >  6:148059/1‑149 (MQ=255)                                                                                                                    ggTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCa                             <  8:309718/148‑1 (MQ=255)                                                                                                                     gTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAg                            >  3:278546/1‑148 (MQ=255)                                                                                                                      tGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCa                          >  2:44340/1‑149 (MQ=255)                                                                                                                        cAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCAt                         <  5:148059/148‑1 (MQ=255)                                                                                                                        cAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATg                        >  6:53752/1‑149 (MQ=255)                                                                                                                         aGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGt                       >  6:234841/1‑149 (MQ=255)                                                                                                                         aGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGt                       >  8:348829/1‑149 (MQ=255)                                                                                                                                   gggCCCTTTTCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGcttctt             >  3:75119/2‑149 (MQ=255)                                                                                                                                   cGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGcttctt             <  2:131755/149‑1 (MQ=255)                                                                                                                                     gcccttttCGCGGGCTTGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGcttctt             >  8:336221/6‑148 (MQ=255)                                                                                                                                     gCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa           <  1:40524/149‑1 (MQ=255)                                                                                                                                      cccttttCGCGGGCTTGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa           <  7:336221/145‑1 (MQ=255)                                                                                                                                       ccttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGc         >  1:151771/4‑149 (MQ=255)                                                                                                                                          tttCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc       >  6:294160/1‑149 (MQ=255)                                                                                                                                               cgGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCcgcg  >  5:328203/1‑149 (MQ=255)                                                                                                                                                   cAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGtt                                        <  1:57923/107‑1 (MQ=255)                                                                                                                                                   cAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGtt                                        >  2:57923/1‑107 (MQ=255)                                                                                                                                                     |                                                                                                                                                AACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGC‑TAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCG  >  NC_002947/2068005‑2068294

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