Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,199,645	G→C	11.3%	intergenic (+20/+53)	PP_4579 → / ← PP_4580	LysR family transcriptional regulator/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,199,645	0	G	C	11.3%	118.9 / 5.1	44	intergenic (+20/+53)	PP_4579/PP_4580	LysR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base C (2/3);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCCTTGAGCTGGTGGTGGCGTGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAGGCCAGC  >  NC_002947/5199499‑5199788                                                                                                                                                   |                                                                                                                                                gAGCCTTGAGCTGGTGGTGGCGTGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCGGTAGCg                                                                                                                                               >  4:215741/1‑149 (MQ=255)    ccTTGAGCTGGTGGTGGCGTGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGcc                                                                                                                                            <  6:117181/149‑1 (MQ=255)    ccTTGAGCTGGTGGTGGCGTGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGc                                                                                                                                             <  4:178640/148‑1 (MQ=255)         aGCTGGTGGTGGCGTGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTTGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTc                                                                                                                                       <  6:314168/149‑1 (MQ=255)         aGCTGGTGGTGGCGTGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTc                                                                                                                                       >  3:62616/1‑149 (MQ=255)                ggtggCGTGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGAAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTACCGGCCCTTTCGCGGGCa                                                                                                                                >  6:295278/1‑149 (MQ=255)                 gtggcgtgtcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCACCTGTACCGGCCCTTTCGCGGGCAt                                                                                                                               >  1:356981/1‑149 (MQ=255)                      gtggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCcg                                                                                                                          >  4:145519/1‑149 (MQ=255)                           aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTACCGGCCCTTTCGCGGGCATGCCCGCGaaa                                                                                                                     <  5:295278/149‑1 (MQ=255)                            gACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCACCTGTACCGGCCCTTTCGCGGGCATGCCCGCGAAAg                                                                                                                    <  2:356981/149‑1 (MQ=255)                               cGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGaaa                                                                                                                     >  2:73741/1‑145 (MQ=255)                               cGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGaaa                                                                                                                     <  1:73741/145‑1 (MQ=255)                                 gggTGGCGGTGGGCGATGAGGTGGGTGGGGGTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGAGGTGGGGGGGCAGCGGATCGTGCTGGGCTGACCAACAGGCTCCGCCGGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCg                                                                                                               <  3:215741/149‑1 (MQ=38)                                      gCGTTGGTCGATGCGGGGGTTGGGGTTGCCCGGCTGGGGTGGGCGCGGTCTGCGCGGGTGGTGGGTGGGCAGCGGATCGGGCTGGCCTTACCAACAGTCTCCGCCTGTCGCGCCCCTTTCGCGGGCATGCCCGGGAAAGGGCCggtaca                                                                                                          <  2:304919/149‑1 (MQ=255)                                        gTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCggtacagg                                                                                                        >  1:2383/1‑149 (MQ=255)                                            gTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTTGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTcgcg                                                                                                                                    <  5:173508/117‑1 (MQ=255)                                            gTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTcgcg                                                                                                                                    >  6:173508/1‑117 (MQ=255)                                              cGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCATTCGCGGGCATGCCCGCGAAATGGCCCgtacaggtacag                                                                                                   >  4:330377/1‑148 (MQ=255)                                                                 tGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGctc                                                                                >  7:62172/1‑148 (MQ=255)                                                                    cGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCggtacaggt                                                                                                       <  5:95401/122‑1 (MQ=255)                                                                    cGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCggtacaggt                                                                                                       >  6:95401/1‑122 (MQ=255)                                                                      gCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTc                                                                          <  8:62172/149‑1 (MQ=255)                                                                        aGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGccc                                                                                                                           <  2:231062/98‑1 (MQ=255)                                                                        aGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGccc                                                                                                                           >  1:231062/1‑98 (MQ=255)                                                                        aGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTcac                                                                        >  4:244907/1‑149 (MQ=255)                                                                                 cgcgGTAGGCGCGTGATGTGGGTGGGCATCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTgcg                                                               <  3:330377/149‑1 (MQ=255)                                                                                        tGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGCCCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGCAAGGGCCGGGACAGGTACAGCAGAACTTAGGGCCAGCTGTCTTTCACCACCCAGCGccgcccc                                                        >  6:346420/1‑148 (MQ=39)                                                                                        tGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgcccg                                                        >  8:264633/1‑149 (MQ=255)                                                                                        tGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgcccg                                                        >  5:147872/1‑149 (MQ=255)                                                                                          cgcgGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCttt                                                                                                                                        >  7:249802/1‑67 (MQ=255)                                                                                          cgcgGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCttt                                                                                                                                        <  8:249802/67‑1 (MQ=255)                                                                                          cgcgGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCATTGGCGGGCATGCCCGCGAAAGGGCCGGGACAGGTACCGCAGAACTTAGCGCCAGCGCACGTTCACCACCCGGCGCcgcccgcc                                                      >  2:316691/1‑149 (MQ=255)                                                                                             gggATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAg                                                   >  5:205726/1‑149 (MQ=255)                                                                                                     ggtgggCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCCTGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagc                                            >  5:330841/1‑148 (MQ=255)                                                                                                         ggCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagccagc                                        <  3:145519/148‑1 (MQ=255)                                                                                                             gCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATc                                   <  6:147872/149‑1 (MQ=255)                                                                                                                      gCTGGTGTGACCAACAGTCCCCGCCTGTTCCGGCCCTTTCGCGGGCAAGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagccagc                                        <  5:13414/135‑1 (MQ=39)                                                                                                                      gCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCaa                           >  3:286106/1‑148 (MQ=255)                                                                                                                               aCCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTc                  <  3:244907/148‑1 (MQ=255)                                                                                                                                       tctcCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGATGGTAAAGGTAGAGCAGAACTTAGCGCGAGCTCTCTTTCTCAACCCTGCGCCGCACGACCAGGATGAGCCAGGCCATCCCCAGCAACAGGCTTTCGACcaca           >  5:351530/1‑147 (MQ=37)                                                                                                                                           cGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAg        <  7:264633/146‑1 (MQ=255)                                                                                                                                           cGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAGGc      <  4:286106/148‑1 (MQ=255)                                                                                                                                              cTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAGGCCAGc  <  6:205726/149‑1 (MQ=255)                                                                                                                                              cTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAGGCCAGc  <  1:316691/149‑1 (MQ=255)                                                                                                                                                   gCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACaaa         >  3:190599/1‑137 (MQ=255)                                                                                                                                                   gCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACaaa         <  4:190599/137‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GAGCCTTGAGCTGGTGGTGGCGTGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAGGCCAGC  >  NC_002947/5199499‑5199788

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