Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,068,152	A→T	18.5%	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	18.5%	86.9 / 13.8	43	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 (19/16);  new base T (6/2);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

AAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCG  >  NC_002947/2068011‑2068299                                                                                                                                               |                                                                                                                                                    aaGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCcgc                                                                                                                                               >  7:140806/1‑148 (MQ=255)   gAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCcgcg                                                                                                                                              >  2:222807/1‑148 (MQ=255)   gAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCAGGCAAGCCCGCGa                                                                                                                                             >  3:35361/1‑148 (MQ=255)     aaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                          <  7:96916/149‑1 (MQ=255)       aTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGaaagg                                                                                                                                         <  1:222807/149‑3 (MQ=255)             ccccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGt                                                                                                                                  >  7:236673/1‑149 (MQ=255)             ccccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGt                                                                                                                                  >  1:34667/1‑149 (MQ=255)                     ggAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAg                                                                                                                          <  2:34667/149‑1 (MQ=255)                            aCGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGAGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAatat                                                                                                                   >  4:109720/1‑149 (MQ=255)                                aaTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCTGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                               <  4:35361/149‑1 (MQ=255)                                 aTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctgc                                                                                                              <  3:101741/149‑1 (MQ=255)                                  tGACTAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATgctgct                                                                                                             <  1:316473/149‑1 (MQ=255)                                   gAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATgctgct                                                                                                             >  4:305072/1‑148 (MQ=255)                                    aaTAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGctgctg                                                                                                            >  7:205993/1‑148 (MQ=255)                                    aaTAAAGCATGGTACGACTCACACAAGACGGACGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGAAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGctgctg                                                                                                            >  7:210103/1‑148 (MQ=255)                                                     cTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCTTGCCCGCGaaagg                                                                                                                                         <  8:26878/102‑3 (MQ=255)                                                     cTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCTTGCCCGCGaaagg                                                                                                                                         >  7:26878/1‑100 (MQ=255)                                                     cTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCtt                                                                                          >  2:209999/1‑149 (MQ=255)                                                      tCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTg                                                                                         >  2:3247/1‑149 (MQ=255)                                                             aGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGAGCGGCTAGCCAata                                                                                                                    >  6:157999/1‑115 (MQ=255)                                                                 ggTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                               >  6:296839/1‑116 (MQ=255)                                                                 gcgcgcgGGGCAGTGTACGGGTGGGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTc                                                                              <  4:236528/146‑1 (MQ=255)                                                                   tCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGcc                                                                                >  5:114717/1‑145 (MQ=255)                                                                   tCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGcc                                                                                <  6:114717/145‑1 (MQ=255)                                                                                               gctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                                  <  3:305072/147‑1 (MQ=255)                                                                                                     gCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCtt                                          <  3:109720/149‑1 (MQ=255)                                                                                                       aaaCCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTtgt                                        <  1:3247/149‑1 (MQ=255)                                                                                                           cTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgc                                                                <  3:112296/121‑1 (MQ=255)                                                                                                           cTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgc                                                                >  4:112296/1‑121 (MQ=255)                                                                                                                     gCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGtt                           >  3:54688/1‑148 (MQ=255)                                                                                                                     gCCTGTTACGGCCC‑TTTCGCGCGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTc                          <  8:74101/149‑1 (MQ=255)                                                                                                                          ttACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGct                      <  8:205993/148‑1 (MQ=255)                                                                                                                          ttACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGct                      <  8:210103/148‑1 (MQ=255)                                                                                                                          ttACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGc                       <  8:236673/147‑1 (MQ=255)                                                                                                                           tACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGctt                     <  1:209999/148‑1 (MQ=255)                                                                                                                            aCGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTTCGCACTTTCAGCTTGGTGTTATTGTGCGCCAGCATGTTCAGcttct                   >  5:58366/1‑149 (MQ=255)                                                                                                                              ggCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGcttctt                  >  7:161095/1‑148 (MQ=255)                                                                                                                              ggCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGcttctt                  >  3:310947/1‑148 (MQ=255)                                                                                                                                ccc‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAg               >  5:32277/1‑149 (MQ=255)                                                                                                                                 ccttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGc              >  5:42929/4‑149 (MQ=255)                                                                                                                                  cttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCt             >  2:201126/3‑149 (MQ=255)                                                                                                                                   ttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc            >  7:229248/2‑149 (MQ=255)                                                                                                                                   ttttCGCGGGCTTGCCCGCG‑AAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc            >  3:130973/2‑149 (MQ=255)                                                                                                                                              aaGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCg  <  8:161095/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    AAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCC‑TTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCG  >  NC_002947/2068011‑2068299

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