Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,499,508	+C	intergenic (+85/‑79)	trpS → / → zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,499,506	1	.	C	95.9%	137.3 / ‑0.2	49	intergenic (+83/‑81)	trpS/zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (24/23);  total (25/24)
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

GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCCA‑C‑GGCTTTGCCGGTGTT‑TCGGG‑CACG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC  >  NC_002947/1499361‑1499648                                                                                                                                                        |                                                                                                                                                gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTAGCGGG‑CAcg‑cc                                                                                                                                               >  8:4744/1‑149 (MQ=37)            tGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCTCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCcaca                                                                                                                                     <  1:299499/149‑1 (MQ=21)                 cGGCCATTCAGAGCCTCAACGAAGCCTTGGCACCACAGCAGGACTGTGGCTGCCTAACGGTGGGAGGGGGCGTGCCCGCGAAACAGGCGG‑CGCGGTGAATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGa                                                                                                                               <  1:274379/149‑1 (MQ=21)                   gCCATTCAGCGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGAtt                                                                                                                             <  3:31486/149‑1 (MQ=21)                          agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑ACTCCCACAGACTGATTGCCAtt                                                                                                                       <  3:330634/148‑1 (MQ=21)                                 gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGCCCGCGAAACAGCCGC‑CCCGGTGCATGGCAC‑C‑GGCTTTGACGGTGTTCGCGGG‑CACGCCC‑GCTCCCACACACAGATTGCCAATAACCGGGc                                                                                                               >  5:336191/1‑149 (MQ=18)                                 gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACAGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCGTTGCCCGTGTTAGCGGG‑CACGCCC‑GCTCCCACCGACTGATTGCCCTTTAGCGGGc                                                                                                               >  1:273026/1‑149 (MQ=21)                                      aaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                                                          <  1:310546/149‑1 (MQ=21)                                             ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTAGCGGG‑GACGCCC‑GCTCCCCCAGACTGATTGCCATTAAGCGGGCCTATCGCTACAg                                                                                                   >  8:285146/1‑149 (MQ=25)                                                accacAGCAGGACTTAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                                                <  3:40447/149‑1 (MQ=21)                                                accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAACGCGC‑CCCGGTGCCTTGCAC‑C‑CGCTTTTGCGGTGTTCGCGGGCCCCC‑CC‑GCTCCCCCAGACTGATTGGCCTTAATCGGGCCAATCGGTACAGTgg                                                                                                >  8:232790/1‑148 (MQ=16)                                                          gccTGAGGCTCCCAAAGGTTGGGAGCGGGCGTGCCCGCTAAACAGGCGC‑CGCGGTGCATGGAAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGa                                                                                                                               <  1:120108/106‑1 (MQ=14)                                                          gACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGGATGGCAC‑C‑CGCATTGCCGGTGTTGGCGGG‑CCCGCCC‑CCGCCCACAGACTGATTGCCAATAAGCGGGGCTATAGCGACAGTTACGGCCCGgtt                                                                                      >  4:84180/1‑149 (MQ=18)                                                            cTGCGGCTGCAAAAGGGTGGGAGGGGGCGTGCCCGCTAAACCGGCGC‑CGCGG‑GGATGGTCACC‑GGTTTTGCGGGTGTTCGCGGG‑AACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                                                    <  8:241208/149‑1 (MQ=255)                                                                      cAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCCC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CCCGCCC‑GCGCCCCCAGACTGCTTGCCATTAAGCGGGCCTATGGCGACAGTGACGGCCCGGTTTTTGTTGCcgcg                                                                          >  3:61986/1‑149 (MQ=18)                                                                      cAAACGGTGGGAGCGGGCGTGACCGCGAAACAGGCGC‑CGCCGTGCATGGCAC‑C‑CGCGTTGCCCGTGTTACCGCG‑CACGCCC‑GCTCCCACAGACTGAATGCCATTAAGCGGGGCTAAGGCTACAGTGAAGGCCCGTTTTTTGTTGCCACg                                                                          >  8:106976/1‑149 (MQ=255)                                                                         aCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                                                <  3:284964/124‑1 (MQ=18)                                                                         aCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCCGTGCATGGCAC‑C‑CGCATTGCCCGTGTTAGCCGG‑CACGCCC‑CCGCCCCCAGACAGATTGCCATTAAGCGGGCCGATAGCGACAGTgg                                                                                                >  4:284964/1‑123 (MQ=12)                                                                                 aGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                               >  5:125995/1‑149 (MQ=25)                                                                                 aGCGGGCGTCCCGGCTAACCCGGCGC‑CGCGGTGCA‑GGCCA‑CGGGCTTTGCCGGTGTT‑CGGGGCCACG‑CCGGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                               <  6:220837/149‑1 (MQ=255)                                                                                  gCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCGGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                              <  3:84180/149‑1 (MQ=25)                                                                                  gCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                              <  7:4744/149‑1 (MQ=25)                                                                                         tGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGCCGGCCCGTTTTTTGTTGCCTCGCTAACGACTCATGACTccc                                                       >  4:209523/1‑149 (MQ=25)                                                                                          gCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                      >  1:260380/1‑149 (MQ=37)                                                                                           cccGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                      >  8:104788/1‑148 (MQ=37)                                                                                              gcgAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                                                   >  8:297421/1‑148 (MQ=37)                                                                                                gAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                       >  5:223380/1‑142 (MQ=37)                                                                                                gAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                       <  6:223380/142‑1 (MQ=37)                                                                                                 aaaCAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                                               <  6:125995/149‑1 (MQ=37)                                                                                                     aGGTGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                                >  6:307118/1‑112 (MQ=21)                                                                                                     aGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                                <  5:307118/112‑1 (MQ=21)                                                                                                       gcgc‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                          >  7:287077/1‑148 (MQ=37)                                                                                                         gc‑cgcGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                       >  8:73905/1‑149 (MQ=37)                                                                                                              cgGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                       >  7:186512/1‑145 (MQ=37)                                                                                                               ggTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                                  <  7:104788/149‑1 (MQ=38)                                                                                                                          acc‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc                        <  7:232790/147‑1 (MQ=38)                                                                                                                          acc‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc                        >  1:154898/3‑149 (MQ=38)                                                                                                                            c‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc                      <  7:73905/149‑1 (MQ=39)                                                                                                                                ctatcccGGTGTTCGCGCG‑CACGCCC‑GCTCCCACAGACTGAATGCCATTAAGCGGGCCAATCGCTACAGTGACGGCCCGGATTATGATGCCTCGCTACCGCAACAAGACACCCAGAGAACGCTATCACGCAGATCTTCACCGTcccca                    >  8:168460/6‑146 (MQ=255)                                                                                                                                cTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc                   <  7:297421/149‑1 (MQ=39)                                                                                                                                cTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc                   <  8:314896/149‑1 (MQ=39)                                                                                                                                 tttGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                                    >  4:136164/1‑67 (MQ=18)                                                                                                                                 tttGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                                    <  3:136164/67‑1 (MQ=18)                                                                                                                                     ccGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt              >  5:132424/1‑149 (MQ=255)                                                                                                                                     ccGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt              <  1:359268/149‑1 (MQ=255)                                                                                                                                     ccGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct               <  7:285146/148‑1 (MQ=255)                                                                                                                                             cgcGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCa           >  6:336690/3‑144 (MQ=255)                                                                                                                                             cgcGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCa           <  5:336690/142‑1 (MQ=255)                                                                                                                                             cgcGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAACTGAAACGACCCGACTTCTTCCATGAct      >  2:149924/3‑148 (MQ=255)                                                                                                                                                 gg‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcggc  >  1:352723/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCCA‑C‑GGCTTTGCCGGTGTT‑TCGGG‑CACG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC  >  NC_002947/1499361‑1499648

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