Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,814,541	A→T	10.0%	intergenic (+41/‑72)	frmC → / → ispF	S‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase/2‑C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,814,541	0	A	T	10.0%	130.1 / 3.6	48	intergenic (+41/‑72)	frmC/ispF	S‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase/2‑C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base T (3/2);  total (25/23)
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.00e+00

AATGACCCTGCGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGC  >  NC_002947/1814395‑1814682                                                                                                                                                   |                                                                                                                                               aaTGACCCTGCGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAcag                                                                                                                                               <  5:217895/148‑1 (MQ=255) aaTGACCCTGCGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAcag                                                                                                                                               >  6:217895/1‑148 (MQ=255) aaTGACCCTGCGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAcag                                                                                                                                               >  1:271376/1‑148 (MQ=255)       ccTGCGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGccc                                                                                                                                         >  4:118848/1‑148 (MQ=255)        cTGCGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAGGGGCTGCTGTGCAGccctt                                                                                                                                       >  1:242734/1‑147 (MQ=255)        cTGCGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCAGCGAAGGGCTGCACAGAAGcccc                                                                                                                                        >  5:190552/1‑148 (MQ=255)         tGCGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCaa                                                                                                                                      >  5:77685/1‑149 (MQ=255)           cGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCaa                                                                                                                                      <  4:172412/147‑1 (MQ=255)           cGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATc                                                                                                                                    <  3:244353/149‑1 (MQ=255)              aTGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAGGGGCTGCTGTGCAGcccttcgcgg                                                                                                                                  <  2:242734/148‑8 (MQ=255)               tGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCca                                                                                                                                 <  6:77685/148‑1 (MQ=255)                 cAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGAGAAGGGCTGCACAGCAGCCCCAATCCcaca                                                                                                                               <  4:187326/148‑1 (MQ=255)                                    gCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAACCAGGTa                                                                                                           <  2:271376/149‑1 (MQ=255)                                         tactTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTa                                                                                                           <  3:86710/144‑1 (MQ=255)                                         tactTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTa                                                                                                           >  4:86710/1‑144 (MQ=255)                                               aTTGCCTGATTGATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑cc                                                                                                <  6:190552/148‑1 (MQ=255)                                                ttGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCt                                                                                              >  5:5307/1‑149 (MQ=255)                                                  gCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCt                                                                                                                       <  3:124903/123‑1 (MQ=255)                                                  gCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCt                                                                                                                       >  4:124903/1‑123 (MQ=255)                                                            tCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAgg                                                                                  >  4:53298/1‑149 (MQ=255)                                                            tCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCACTGCAAAGAAGCTAGAATCACG‑CCATGACTCCTTAAgg                                                                                  >  4:163070/1‑149 (MQ=38)                                                                    cACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTc                                                                          <  4:226751/149‑1 (MQ=255)                                                                      ccTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCtt                                                                        <  3:118848/149‑1 (MQ=255)                                                                       cTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTg                                                                                                                      >  2:201734/1‑103 (MQ=255)                                                                       cTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTg                                                                                                                      <  1:201734/103‑1 (MQ=255)                                                                         gcgcCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGc                                                                                >  5:199572/1‑138 (MQ=255)                                                                         gcgcCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGc                                                                                <  6:199572/138‑1 (MQ=255)                                                                          cggcataatGCGGTGGCACTGGGTCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGGACATGAGCCCCAATGCCCCAAATGCCCTGCAAATCCGGTAGAATCACGACCCTGACTCATTCAGGGCGTTTTCTTAt                                                                      <  1:153588/145‑1 (MQ=255)                                                                            ccATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCg                                                                   >  6:149114/1‑148 (MQ=255)                                                                                  tGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGc                                                            >  6:201065/1‑149 (MQ=255)                                                                                               gggCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGa                                                 <  6:5307/147‑1 (MQ=255)                                                                                                   gggTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCa                                           >  5:100189/1‑149 (MQ=255)                                                                                                     gTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCAcc                                         >  5:225503/1‑149 (MQ=255)                                                                                                     gTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCAcc                                         >  7:48031/1‑149 (MQ=255)                                                                                                            gCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTg                                  >  1:135859/1‑149 (MQ=255)                                                                                                            gCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTAGTGGCCACGGCTACGATGTGCACCGTTTCTg                                  >  3:39300/1‑149 (MQ=255)                                                                                                            gCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGAATGAGTCAGGGCGtga                                                                            >  6:59523/1‑105 (MQ=37)                                                                                                            gCCTCCTGTGGGAGCGGCCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGAATGAGTCAGGGCGtga                                                                            <  5:59523/107‑3 (MQ=37)                                                                                                              ctcctGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCt                                   >  8:47307/1‑146 (MQ=255)                                                                                                                    tgGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGt                           <  8:48031/148‑1 (MQ=255)                                                                                                                         gCGGGCGTGCCCGCTAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTc                     <  7:47307/149‑1 (MQ=255)                                                                                                                            ggCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCAtt                  <  3:53298/149‑1 (MQ=255)                                                                                                                             gCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTAAGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTa                 >  1:219771/1‑149 (MQ=255)                                                                                                                                  cccGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCt             <  3:163070/148‑1 (MQ=255)                                                                                                                                    cgcgAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTggg          >  1:16330/1‑149 (MQ=255)                                                                                                                                     gattggggCTGCTGTGCAGCCCCAATCACACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCGGACTCAGTCAGGGCGTTTTCTTATGGGGTTTGGCCACGGGTAAGATGTGCACAGTTTCTGCGACGGTGATTTCATTAGCGTGGGc         >  4:218669/6‑149 (MQ=255)                                                                                                                                      attggggCTGCTGTGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCg        >  8:92900/5‑149 (MQ=255)                                                                                                                                         aGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCgggg     >  7:47675/1‑149 (MQ=255)                                                                                                                                            gctgctCTGCAGCCCCAATCCCACAACTCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGc  <  3:218669/149‑1 (MQ=255)                                                                                                                                            gCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGc  <  4:507/149‑1 (MQ=255)                                                                                                                                            gCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGc  <  2:219771/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               AATGACCCTGCGTATGCAGCCAGGCTATGACCATAGCTACTACTTCATTGCCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACG‑CCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCGTTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGC  >  NC_002947/1814395‑1814682

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