Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,178,752	G→T	11.5%	intergenic (+57/+37)	guaA → / ← PP_1033	glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,178,752	0	G	T	11.5%	143.9 / 9.8	52	intergenic (+57/+37)	guaA/PP_1033	glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base T (3/3);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

GTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTT  >  NC_002947/1178613‑1178897                                                                                                                                              |                                                                                                                                                  gTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATAGCGGGc                                                                                                                                          >  3:77385/1‑149 (MQ=255)      atcaACGAACTCGACGGTATCTCGGGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGGGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCAATCATGCGGGCTTTCt                                                                                                                                     <  7:287601/149‑1 (MQ=255)      atcaACGAAATCGACGGTATCTCGTGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCAATCATGCGGGCTTTCt                                                                                                                                     <  7:132087/149‑1 (MQ=255)      atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCt                                                                                                                                     <  3:209296/149‑1 (MQ=255)       tcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                    <  8:90921/149‑1 (MQ=255)       tcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                    >  6:153399/1‑149 (MQ=255)         aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAATCCCGCATGATTGCGGGCTTTCTGGt                                                                                                                                  >  1:313252/1‑149 (MQ=255)         aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGt                                                                                                                                  >  4:141322/1‑149 (MQ=255)         aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGCAC‑GGATA‑CACCAGAAAGCCCGCATGAGTGCGGGCTTTCTGGt                                                                                                                                  >  8:216379/1‑149 (MQ=255)          aCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGtt                                                                                                                                 <  3:91572/149‑1 (MQ=255)                  gACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATACCACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttc                                                                                                                          <  1:229409/149‑1 (MQ=255)                           tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCAATCATGCGGGCTTTCTGGtgtat                                                                                                                              >  8:221498/1‑131 (MQ=38)                           tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCAATCATGCGGGCTTTCTGGtgtat                                                                                                                              <  7:221498/135‑5 (MQ=38)                             gcgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTg                                                                                                                <  2:180513/147‑1 (MQ=255)                               gcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCAGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCATTCTGGTTCTTGTTCGTTAAGATTGCTTa                                                                                                            >  3:44631/1‑149 (MQ=255)                                     ccTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAttcttc                                                                                                      >  8:247057/1‑149 (MQ=255)                                        acgacgTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGAACAGTCTGGCACCCGCAGCAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATAATTCTgg                                                                                                   >  7:139292/1‑149 (MQ=255)                                               gtCGAGCAAGCCGCCTGCCACTATCTAGTTGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCt                                                                                             <  2:313252/148‑1 (MQ=255)                                                      aaGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCaaa                                                                                       >  7:5828/1‑147 (MQ=255)                                                         ccgccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACTGGAAA‑AACCACAAACCCCGCATGATTGCGGGCTTACTGGTTGTTTTACGTTAAAATTGCTTATTATTATGCCTTTATGGCAAAAAgc                                                                                   >  1:259982/1‑149 (MQ=255)                                                              tGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTc                                                                             >  5:212399/1‑149 (MQ=255)                                                                        gAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGt                                                                   <  6:212399/149‑1 (MQ=255)                                                                           tGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCTGTTTGCTTATTCTTCTGGCTTTCTGGCATATAGCGCTTTCTGGTAGGCGTta                                                                 <  4:77385/148‑1 (MQ=255)                                                                             ggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtata                                                               <  7:247057/148‑1 (MQ=255)                                                                              gAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGtt                                                                  <  1:226674/144‑1 (MQ=255)                                                                              gAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGtt                                                                  >  2:226674/1‑144 (MQ=255)                                                                              gAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAg                                                              >  1:138650/1‑148 (MQ=255)                                                                                gTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGa                                                            >  3:193298/1‑148 (MQ=255)                                                                                   aTCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAag                                                         >  4:26710/1‑148 (MQ=255)                                                                                    tCGGGTAGCCAGCCCAGGCTGGCACCCGGAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAagca                                                       <  7:169641/149‑1 (MQ=255)                                                                                        gTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCt                                                   <  3:141322/149‑1 (MQ=255)                                                                                            gcagacCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCGCTAATGATTGCGGGCTTTCTGGTTCTAGTTCGTTCAGATTTCTTATTCTTCTGGCTTTCTGTAAAATAGCTCTTTGTGGTAGGCGTTATAGGTAAGCAGACTGctg                                               <  8:139292/148‑2 (MQ=255)                                                                                             cagccCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCt                                                                                             >  4:96908/1‑102 (MQ=255)                                                                                             cagccCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCt                                                                                             <  3:96908/102‑1 (MQ=255)                                                                                             cagccCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctcc                                              >  6:291225/1‑149 (MQ=255)                                                                                                      cTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGa                                     <  5:153399/149‑1 (MQ=255)                                                                                                            cccGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACt                               >  4:41520/1‑149 (MQ=255)                                                                                                                 aGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTg                                        <  6:308381/135‑1 (MQ=255)                                                                                                                 aGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTg                                        >  5:308381/1‑135 (MQ=255)                                                                                                                   gAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCAcc                        >  4:215405/1‑149 (MQ=255)                                                                                                                   gAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATAGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCAcc                        >  8:288220/1‑149 (MQ=255)                                                                                                                           agaacCAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGAAAATAGAGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCAc                  >  7:37395/4‑149 (MQ=38)                                                                                                                                cAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCATCCACGCGtt             >  2:31797/1‑149 (MQ=255)                                                                                                                                  gAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACTCGTTgg           >  5:74225/1‑149 (MQ=255)                                                                                                                                   aaaGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggg          <  7:216379/149‑1 (MQ=255)                                                                                                                                    aaGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgggg         <  7:288220/149‑1 (MQ=255)                                                                                                                                        ccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGc     <  5:291225/149‑1 (MQ=255)                                                                                                                                        ccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGc     >  1:212273/1‑149 (MQ=255)                                                                                                                                         cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTc                           <  4:147786/126‑1 (MQ=255)                                                                                                                                         cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTc                           >  3:147786/1‑126 (MQ=255)                                                                                                                                         cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTgac    <  8:5828/149‑3 (MQ=255)                                                                                                                                           cATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCttt  <  3:41520/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                  GTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAC‑GGATA‑CACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTT  >  NC_002947/1178613‑1178897

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