Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,111,315	A→C	9.8%	intergenic (+36/‑74)	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,315	0	A	C	9.8%	115.9 / 5.4	41	intergenic (+36/‑74)	PP_0970/PP_0971	cyanate MFS transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base C (2/2);  total (23/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 = 7.14e-01

GGCTGGACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTG  >  NC_002947/1111169‑1111463                                                                                                                                                   |                                                                                                                                                     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                                                              >  3:39561/1‑142 (MQ=255)                                                                                              aCCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAAc                                                              <  4:39561/142‑1 (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                                          <  4:88051/121‑1 (MQ=255)                                                                                                                                       tCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAACCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGttt                                          >  3:88051/1‑121 (MQ=255)                                                                                                                                                ccAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGc               >  7:214433/1‑139 (MQ=255)                                                                                                                                                ccAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGc               <  8:214433/139‑1 (MQ=255)                                                                                                                                                ccAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATg     >  6:253188/1‑149 (MQ=255)                                                                                                                                                   agagGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTg  <  8:226936/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GGCTGGACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTG  >  NC_002947/1111169‑1111463

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