Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,111,320	G→A	9.1%	intergenic (+41/‑69)	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,320	0	G	A	9.1%	119.9 / 2.6	44	intergenic (+41/‑69)	PP_0970/PP_0971	cyanate MFS transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/18);  new base A (2/2);  total (24/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 = 9.82e-01

GGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTGGCCT  >  NC_002947/1111177‑1111467                                                                                                                                                |                                                                                                                                                    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                         <  6:13978/142‑1 (MQ=255)                                                                                                                               tCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCGGCGCCt                         >  5:13978/1‑142 (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)                                                                                                                                           agagGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTg      <  8:226936/149‑1 (MQ=255)                                                                                                                                            gagGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGAt                                                                        >  6:184922/1‑82 (MQ=255)                                                                                                                                            gagGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGAt                                                                        <  5:184922/82‑1 (MQ=255)                                                                                                                                               ggCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTGGCCt  >  4:146392/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    GGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGTGCAGTTTTCGAACCTAGGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTGGCCT  >  NC_002947/1111177‑1111467

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