Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,111,320	G→A	7.7%	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	7.7%	147.2 / 3.1	52	intergenic (+41/‑69)	PP_0970/PP_0971	cyanate MFS transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base A (2/2);  total (25/27)
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.23e-01

GACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTGGC  >  NC_002947/1111174‑1111465                                                                                                                                                   |                                                                                                                                                      gACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCg                                                                                                                                                     <  3:90684/149‑1 (MQ=255) gACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGGTGGCGCCGTGGTGTTCGGTCTGGGGGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGc                                                                                                                                                      <  2:187388/148‑1 (MQ=255)  aCGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCgg                                                                                                                                                    <  5:147007/149‑1 (MQ=255)    ggCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGTCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCggg                                                                                                                                                   >  6:225214/1‑148 (MQ=255)           ggTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGGTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTgg                                                                                                                                              <  2:349336/146‑1 (MQ=255)            gTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCGTGCCAAGAGGGCGGGACTGGCAt                                                                                                                                           >  6:192443/1‑148 (MQ=255)            gTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCAt                                                                                                                                           >  1:21566/1‑148 (MQ=255)            gTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGCCCTCTTGGCAGGCACGAGGACGGGa                                                                                                                                                  <  3:136988/141‑1 (MQ=39)            gTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGCCCTCTTGGCAGGCACGAGGACGGGa                                                                                                                                                  >  4:136988/1‑141 (MQ=39)                gATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGGGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCGTGCCAAGAGGGCGGGACTGGCATAAGcc                                                                                                                                      <  5:192443/149‑1 (MQ=255)                       gCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCggg                                                                                                                                <  4:301156/148‑1 (MQ=255)                        cGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGc                                                                                                                              <  6:389358/149‑1 (MQ=255)                           tGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGTATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCgggg                                                                                                                               >  7:105267/1‑145 (MQ=255)                           tGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGCCCTCTTGGCAGGCACGAGGACgg                                                                                                                                                    >  6:75824/1‑124 (MQ=38)                           tGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGCCCTCTTGGCAGGCACGAGGACgg                                                                                                                                                    <  5:75824/124‑1 (MQ=38)                             cTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGAt                                                                                                                          >  6:88099/1‑148 (MQ=255)                                ggCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATc                                                                                                                      <  3:152695/149‑1 (MQ=255)                                ggCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATc                                                                                                                      <  4:257296/149‑1 (MQ=255)                                  cGTTGGCGCCGTGGTGTTCGGGCT‑GGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑g                                                                                                                 <  5:20076/149‑1 (MQ=255)                                       gcgcCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGt                                                                                                             >  2:5644/1‑149 (MQ=255)                                           cGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTc                                                                                                         <  5:88099/149‑1 (MQ=255)                                               gtgtTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAAc                                                                                                     <  5:38687/149‑1 (MQ=255)                                                         tGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑ggacgga                                                                                         <  4:321154/149‑1 (MQ=255)                                                           ggTGCCGGGCGGGCGCTGCATGTGCAGTTGACCTGCGAAAAAGTCTGAGGCGTGTTTTAGCAGTCCCGTCCTCGTGCCTGACAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑ggacgga                                                                                         <  4:147634/147‑1 (MQ=255)                                                                    cgggcgCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCaa                                                                              <  5:225214/149‑1 (MQ=255)                                                                     tgtgtcTGCATGTGCATGTTACCTGTGAAAAAGTCTGATGCCTGTTTTTGCAGTCCCGTCCTCGTGCCTGGCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCaa                                                                              <  4:248643/143‑1 (MQ=255)                                                                     gggcgCTGGATGTGCTGGTGACGTGCGAAAAAGTCTGAGGCGTGTTTTAGCAGTGCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATTAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAAt                                                                             <  4:272006/149‑1 (MQ=255)                                                                        cgcTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAg                                                                          <  3:301487/149‑1 (MQ=255)                                                                         gcTGCATGTGCAGGTGACCTGTGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑A‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCaa                                                                              >  3:248643/1‑144 (MQ=255)                                                                             cATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATg                                                                     >  4:52294/1‑149 (MQ=255)                                                                                   cAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAAcc                                                               >  8:168010/1‑149 (MQ=255)                                                                                         aCCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGcc                                                          >  1:294990/1‑148 (MQ=255)                                                                                          ccTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCt                                                        <  8:105267/149‑1 (MQ=255)                                                                                            tGCGAACAAGTCTGTGGCGTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGccc                                                         <  6:57395/146‑1 (MQ=255)                                                                                            tGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGccc                                                         >  5:57395/1‑146 (MQ=255)                                                                                              cGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGAt                                                     <  3:52294/148‑1 (MQ=255)                                                                                                aaaaaGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCAc                                                  >  1:348524/1‑149 (MQ=255)                                                                                                 aaaaGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCgg                                                                                                                                 >  7:285640/1‑73 (MQ=255)                                                                                                             ccTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAg                                     >  6:360096/1‑149 (MQ=255)                                                                                                             ccTGTTTTAGCAGACCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTTTCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAg                                     >  6:9215/1‑149 (MQ=255)                                                                                                               tGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGc                                   >  2:366429/1‑149 (MQ=255)                                                                                                                gTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCGATGGCGGTTTTCGAACCTAAGGGACGGACCAGCCCCCAATGAGTGATGCCCACCGCGCCCTGATCACCCGTTTCTAc                                        >  1:315583/1‑147 (MQ=255)                                                                                                                          tCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGcc                        <  2:21566/149‑1 (MQ=255)                                                                                                                            ccGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTg                      >  4:336856/1‑149 (MQ=255)                                                                                                                            ccGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTg                      >  3:270995/1‑149 (MQ=255)                                                                                                                             cGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTg                      >  1:370821/1‑148 (MQ=255)                                                                                                                                        ctgccAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGc           >  8:227027/1‑148 (MQ=255)                                                                                                                                           ccAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCCACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATg       <  1:5644/149‑1 (MQ=255)                                                                                                                                             aaGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGACGACGCCATGCt     >  4:274450/1‑149 (MQ=255)                                                                                                                                              agagGGCGGGACTGGCATAATCCACCGGGGGTGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACAAGGGACAAATGAGTTATGCCAACCGCGCCGTGATAAACCTTTTATAACAGGCTTTGCAGCGCCTGGATGCGGACGCCAttgtg    >  4:146944/1‑145 (MQ=255)                                                                                                                                                agGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGTCGCCATGCttg   <  2:348524/148‑3 (MQ=255)                                                                                                                                                agGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTGGc  <  7:168010/149‑1 (MQ=255)                                                                                                                                                 gggCTGGACTTGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTTATCACCAGTTTCTACCAGTCTTTCCAGCGCCTGGATGCCGACGCCATGCTGGc  <  3:146944/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GACGGCACTGGGTTGGATTTTTGCGGTGCTGGGCGTTGGCGCCGTGGTGTTCGGGCTGGGTGCCGGGCGGGCGCTGCATGTGCAGGTGACCTGCGAAAAAGTCTGAGGCCTGTTTTAGCAGTCCCGTCCTCGTGCCTGCCAAGAGGGCGGGACTGGCATAAGCCACCGGGGCAGATTATCG‑T‑GCAGTTTTCGAACCTA‑‑GGACGGACCAGCCCCCAATGAGTGATGCCAACCGCGCCCTGATCACCCGTTTCTACCAGGCTTTCCAGCGCCTGGATGCCGACGCCATGCTGGC  >  NC_002947/1111174‑1111465

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