Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,111,312	C→G	10.3%	intergenic (+33/‑77)	PP_0970 → / → PP_0971	cyanate MFS transporter/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,111,312	0	C	G	10.3%	109.6 / 4.7	39	intergenic (+33/‑77)	PP_0970/PP_0971	cyanate MFS transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/16);  new base G (2/2);  total (21/18)
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.84e-01

CGGCGGCTGGACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATG  >  NC_002947/1111165‑1111460                                                                                                                                                    |                                                                                                                                                     cggcggCTGGACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTAGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCGGTTTTAGCAGTCCCGTCCTCGtgcctgc                                                                                                                                                      >  4:222699/1‑148 (MQ=255)     ggcTGGACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAaga                                                                                                                                                 <  3:6219/149‑1 (MQ=255)             ggCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGCCCTCTTGGCAGGCACGAGGACGGGa                                                                                                                                         <  3:332208/149‑1 (MQ=255)                 cTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGCCCTCTTGGCAGGCACGAGGACGGGACTgc                                                                                                                                     <  4:195582/149‑2 (MQ=38)                  tGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTgg                                                                                                                                     >  4:260766/1‑148 (MQ=255)                   gggTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGCCCTCTTGGCAGGCACGAGGACGGGACTgct                                                                                                                                    >  5:345352/1‑146 (MQ=38)                    ggTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGc                                                                                                                                    <  3:260766/147‑1 (MQ=255)                      ttGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCAt                                                                                                                                  >  3:188699/1‑147 (MQ=255)                        ggATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAg                                                                                                                               <  1:335155/148‑1 (MQ=255)                                 cGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGc                                                                                                                     <  8:18832/149‑1 (MQ=255)                                        gggCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTAt                                                                                                              >  7:194954/1‑149 (MQ=255)                                         ggCGTTGGCGCGGTGGGGGTGGGGCTGGGGGCGGGGCGGGCGCTGCATGGGCAGGTGACCTGCGAAAACGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATc                                                                                                             <  5:247865/149‑1 (MQ=255)                                            gTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTg                                                                                                          >  6:30181/1‑149 (MQ=255)                                                  gcCGTGGTGTTCGGGCTGGGTGCCGAGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGttt                                                                                                    >  5:227281/1‑149 (MQ=255)                                                       ggggtTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCCGCGAAAAAGGCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTTCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGaa                                                                                               <  7:187432/146‑1 (MQ=255)                                                           ttCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAAcc                                                                                             >  3:274225/1‑147 (MQ=255)                                                                     gTGCCGTGCGGGCGCGGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCa                                                                                 <  8:231658/149‑1 (MQ=255)                                                                                           gCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAAc                                                           >  7:226936/1‑149 (MQ=255)                                                                                                  aCCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAAc                                                           <  4:39561/142‑1 (MQ=255)                                                                                                  aCCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAAc                                                           >  3:39561/1‑142 (MQ=255)                                                                                                    cTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCTAGAGGGCGGGACTGGCATAAGCCACCCGGGCAGATTATCGTGCAGTATTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCCACAGCGCCATg                                                  >  4:263043/1‑149 (MQ=255)                                                                                                    cTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCt                                                   >  5:118179/1‑148 (MQ=255)                                                                                                    cTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTCATGCCAACCGCGCCCt                                                   >  4:22464/1‑148 (MQ=255)                                                                                                    cTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGATTTCGAACCTAGGACGGACCAGCCTCCAATGAGTGATGCCAACCGCGCCCt                                                   >  6:107432/1‑148 (MQ=255)                                                                                                       cGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATc                                               <  8:194954/149‑1 (MQ=255)                                                                                                          aaaaGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCAc                                             <  5:30181/148‑1 (MQ=255)                                                                                                            aaGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCg                                          >  6:134315/1‑149 (MQ=255)                                                                                                                 tGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTc                                      >  3:90047/1‑148 (MQ=255)                                                                                                                     gCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTAcc                                  <  2:32989/148‑1 (MQ=255)                                                                                                                     gCCTGTTTTAGCAGTCCCGTCCTCGTGCCTCCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTAcc                                  >  1:32989/1‑148 (MQ=255)                                                                                                                          ttttAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGc                              <  3:326595/147‑1 (MQ=255)                                                                                                                             tAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTc                          <  5:134315/148‑1 (MQ=255)                                                                                                                               gCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAg                       >  3:364938/1‑149 (MQ=255)                                                                                                                                  gTCCCGCCCTCTTGGCAGGCACGAGGACGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAgcgc                    >  3:109203/1‑149 (MQ=39)                                                                                                                                   tCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGcc                   >  1:232803/1‑149 (MQ=255)                                                                                                                                        tCCTCGTGCCTTCCTAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGAt              <  3:263043/149‑1 (MQ=255)                                                                                                                                           tCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCGGCGCCt                  >  5:13978/1‑142 (MQ=255)                                                                                                                                           tCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCGGCGCCt                  <  6:13978/142‑1 (MQ=255)                                                                                                                                           tCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAACCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGttt                                       >  3:88051/1‑121 (MQ=255)                                                                                                                                           tCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAACCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGttt                                       <  4:88051/121‑1 (MQ=255)                                                                                                                                                    ccAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGc            >  7:214433/1‑139 (MQ=255)                                                                                                                                                    ccAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGc            <  8:214433/139‑1 (MQ=255)                                                                                                                                                    ccAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATg  >  6:253188/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGGCGGCTGGACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATG  >  NC_002947/1111165‑1111460

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