Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,928,379	C→G	10.3%	intergenic (+55/‑85)	peaE → / → PP_3464	phenylacetaldehyde dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,928,379	0	C	G	10.3%	99.8 / 2.9	38	intergenic (+55/‑85)	peaE/PP_3464	phenylacetaldehyde dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (16/18);  new base G (2/2);  total (18/20)
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.82e-01

GGTGGCTACAAGATGTCTGGCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTAC  >  NC_002947/3928232‑3928526                                                                                                                                                     |                                                                                                                                                    ggtggCTACAAGATGTCTGGCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATc                                                                                                                                                     <  1:226908/148‑1 (MQ=255)  gtggCTACAAGATGTCTGGCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATca                                                                                                                                                    >  3:47661/1‑148 (MQ=255)           aGATGTCTGGCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTg                                                                                                                                           >  7:14105/1‑148 (MQ=255)            gATGTCTGGCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCGCCGACACCAGGGATGACATGCCCTggt                                                                                                                                         >  8:243561/1‑147 (MQ=255)            gATGTCTGGCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGACCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTa                                                                                                                                         <  2:274179/149‑1 (MQ=255)             aTGTCTGGCGTTGGTCGTGATATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACAt                                                                                                                                                 <  5:199766/140‑1 (MQ=255)             aTGTCTGGCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACAt                                                                                                                                                 >  6:199766/1‑140 (MQ=255)                    gCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGGGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCgg                                                                                                                                  <  3:252561/148‑1 (MQ=255)                    gCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCgg                                                                                                                                  <  3:114866/148‑1 (MQ=255)                          gTCGTGAGATGGGTTCCGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGGTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTCTTTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCa                                                                                                                             <  1:290549/148‑1 (MQ=255)                          gTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCAc                                                                                                                            <  7:272386/148‑1 (MQ=255)                                      gTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGcc                                                                                                               >  8:157242/1‑149 (MQ=255)                                       ttCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTGAACCCGCTCCTACAGGGCATGTGATACCTGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                                              <  6:81237/149‑1 (MQ=14)                                       ttCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGGGCTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCg                                                                                                              <  7:243561/149‑1 (MQ=255)                                                cATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGc                                                                                                      >  4:255306/1‑148 (MQ=255)                                                                 aGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATa                                                                                    <  7:157242/149‑1 (MQ=255)                                                                         tCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCCGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCAt                                                                            <  8:14105/149‑1 (MQ=255)                                                                            gtgtTGACCCACCCCTCAACCTTTTTCCTGTACCCGCCTC‑AACCCGTGAAAACCCGCTCCAACCCCACGTATCACATGCCCTGTACGAGCGGGATCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGa                                                                         <  4:238632/149‑1 (MQ=25)                                                                                   tCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCCCCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAAccc                                                                  >  2:14975/1‑149 (MQ=255)                                                                                     aaGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACAAGGTATCACATGCCATGTAGGAGAGGGTTCACCCGCGAATAGGCCGGCACTGGCAACAAAAAAAGAAAAATAAAAACCATGGAGCACCCCCa                                                                >  5:286231/1‑149 (MQ=37)                                                                                      aGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACAAGCCCTGTAGGAGCGGGTTCAACCGCGAAGAGGCCGGCACAGGCAACAAAAAAAAAACAAt                                                                                     >  3:119780/1‑127 (MQ=255)                                                                                           tGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTAACCCGCGAAGAGGCCCGCACAGGCAACAAAAAAAAAAAAAAAAAAACCATAGAGCAACCCCATGacg                                                           >  8:143713/1‑146 (MQ=37)                                                                                                   tttCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCacacac                                                    <  3:255306/147‑1 (MQ=255)                                                                                                       cTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGcc                                               <  7:143713/148‑1 (MQ=255)                                                                                                           aCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAgg                                                                                                                 >  3:39269/1‑78 (MQ=255)                                                                                                           aCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAgg                                                                                                                 <  4:39269/78‑1 (MQ=255)                                                                                                             cGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACAt                                        >  8:286242/1‑149 (MQ=255)                                                                                                              ggCCTC‑TTCGCGGGTGAACCCGCTCCTACAGGGCATGTGATACCTGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACAtt                                       >  6:199760/1‑149 (MQ=14)                                                                                                               gcgtc‑tTGGCGGGTAAACCGGCTCCTACACCAGGTATCACATCCCCTGACGGCGCGGTTTAACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTg                                      <  3:48030/146‑1 (MQ=38)                                                                                                                ccTC‑TTCGCGGGTGAACCCGCTCCTACAGGGCATGTGATACCTGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTgg                                     <  5:199760/149‑1 (MQ=14)                                                                                                                     ttCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGAAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATCGAGCAACCCCAAGAGGCACACACACGCCTCCACATTGGCCct                                 >  4:60494/1‑149 (MQ=38)                                                                                                                          gggTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGTAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGt                            <  7:237686/149‑1 (MQ=255)                                                                                                                          gggTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGt                            <  3:60494/149‑1 (MQ=255)                                                                                                                           ggTAAACCCGCTCCTACACCAGGTATCAGATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGt                            >  1:209352/1‑148 (MQ=255)                                                                                                                                  ccGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGcc                    >  4:298808/1‑149 (MQ=255)                                                                                                                                  ccGCTACTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGcc                    >  7:235400/1‑149 (MQ=255)                                                                                                                                           cacCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCgccagc           <  7:286242/149‑1 (MQ=255)                                                                                                                                                 gTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGcc     >  5:83010/1‑149 (MQ=255)                                                                                                                                                   aTCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTa   >  4:363744/1‑149 (MQ=255)                                                                                                                                                    tCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTAc  <  2:209352/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGTGGCTACAAGATGTCTGGCGTTGGTCGTGAGATGGGTTCAGCGGCCATCGAGCATTACACGGAGCTGAAATCGGTGTTGATCAAGCTCTGATCCTTTTTCCTGTACCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACACCAGGTATCACATGCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACAAAAAAACAACAATAAAAACCATGGAGCAACCCCATGAGGCACACACTCGCCTGCACATTGGCCCTGCTGTTAACCGCCCCCGCCAGCCAAGCCTAC  >  NC_002947/3928232‑3928526

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