Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,858,519	A→C	7.8%	intergenic (‑25/+200)	PP_mr50 ← / ← ligA	Bacteria_large_SRP ncRNA/DNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,858,519	0	A	C	7.8%	153.4 / 4.2	51	intergenic (‑25/+200)	PP_mr50/ligA	Bacteria_large_SRP ncRNA/DNA ligase
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base C (2/2);  total (26/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 = 9.87e-01

CGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCA  >  NC_002947/4858377‑4858666                                                                                                                                                |                                                                                                                                                    cGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAAAGCGCGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGagcagc                                                                                                                                                >  4:56813/1‑148 (MQ=255)  gACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGCCCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGagcagc                                                                                                                                                >  8:116426/1‑147 (MQ=255)  gACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGagcagc                                                                                                                                                >  8:325072/1‑147 (MQ=255)  gACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGagcagc                                                                                                                                                >  6:342144/1‑147 (MQ=255)     cccAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCCATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTg                                                                                                                                           >  1:358578/1‑149 (MQ=255)        aGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGgcg                                                                                                                                        >  5:368755/1‑149 (MQ=255)             tGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGTGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCAtt                                                                                                                                   <  8:161210/149‑1 (MQ=255)                  gcggcCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATCGTCTCGCTGAGCAGCGAGTGGCGCCATTGTAc                                                                                                                               <  1:45765/149‑1 (MQ=255)                  gcggcCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGagc                                                                                                                                                   >  6:196659/1‑128 (MQ=255)                    ggcCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCgg                                                                                                                            >  6:295908/1‑149 (MQ=255)                            caGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTg                                                                                                                      >  7:182132/1‑147 (MQ=255)                              gCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGc                                                                                                                  <  6:368755/149‑1 (MQ=255)                               ccACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCa                                                                                                                 <  8:182132/149‑1 (MQ=255)                                 aaacCCCGGACTTCTATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGTGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGa                                                                                                               <  1:320633/147‑1 (MQ=255)                                     cccGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAg                                                                                                                   <  2:370370/141‑1 (MQ=255)                                     cccGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAg                                                                                                                   >  1:370370/1‑141 (MQ=255)                                             tCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGc                                                                                                   >  8:6504/1‑149 (MQ=255)                                                       cGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCa                                                                                               <  7:104261/143‑1 (MQ=255)                                                       cGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCa                                                                                               >  8:104261/1‑143 (MQ=255)                                                             ggCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGtcat                                                                                    <  2:358578/148‑1 (MQ=255)                                                              gctgctCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGtcatca                                                                                  <  7:211588/149‑1 (MQ=255)                                                                   tCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTtct                                                                             >  3:101787/1‑149 (MQ=255)                                                                    cccTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACAGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTtct                                                                             >  8:173274/1‑148 (MQ=255)                                                                     ccTTCCGGGCCTGACCAGGTTAACGGGTAATCTATGCGGGGGGACCGATGGGGCCACCACTACTT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTtctc                                                                            <  3:56813/148‑1 (MQ=255)                                                                      cTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCaa                                                                          <  5:342144/149‑1 (MQ=255)                                                                      cTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCaa                                                                          <  7:325072/149‑1 (MQ=255)                                                                            ggCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACtt                                                                    <  5:295908/149‑1 (MQ=255)                                                                             gCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTg                                                                   >  5:221391/1‑149 (MQ=255)                                                                                  aCCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAgg                                                              >  2:90278/1‑149 (MQ=255)                                                                                   ccAGGTTAACGGGTAATCAATGCGGGTGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAgg                                                              <  7:173274/148‑1 (MQ=255)                                                                                    cAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGtcat                                                                                    >  8:21817/1‑125 (MQ=255)                                                                                    cAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGtcat                                                                                    <  7:21817/125‑1 (MQ=255)                                                                                      ggTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCg                                                                                                                             >  7:90308/1‑82 (MQ=255)                                                                                      ggTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCg                                                                                                                             <  8:90308/82‑1 (MQ=255)                                                                                       gTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCg                                                          <  4:101787/148‑1 (MQ=255)                                                                                              ggTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGa                                                                                                                    >  2:333492/1‑83 (MQ=255)                                                                                                   tCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGa                                                                                                                    <  7:158056/78‑1 (MQ=255)                                                                                                   tCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGa                                                                                                                    >  8:158056/1‑78 (MQ=255)                                                                                                     aaTGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACgg                                           >  5:169421/1‑149 (MQ=255)                                                                                                     aaTGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACgg                                           >  2:73157/1‑149 (MQ=255)                                                                                                     aaTGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACg                                            <  7:6504/148‑1 (MQ=255)                                                                                                           gggggACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCctct                                     <  1:73157/149‑1 (MQ=255)                                                                                                              ggACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTcg                                  >  7:85912/1‑149 (MQ=255)                                                                                                                  cGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCggg                              >  1:96813/1‑149 (MQ=255)                                                                                                                  cGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTAGAGATGGCCTGGACGGGCCTCTTCGCggg                              >  2:255695/1‑149 (MQ=255)                                                                                                                       gggCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAAc                         <  8:84779/149‑1 (MQ=255)                                                                                                                        ggCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAAcc                        >  4:3887/1‑149 (MQ=255)                                                                                                                            accacTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCt                    >  7:293719/1‑149 (MQ=255)                                                                                                                             ccacTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTc                   <  3:3887/149‑1 (MQ=255)                                                                                                                                     gccactcGCTGCTCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCAAGATGGCCt                                                 >  1:341956/5‑112 (MQ=255)                                                                                                                                     gccactcGCTGCTCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCAAGATGGCCt                                                 <  2:341956/108‑1 (MQ=255)                                                                                                                                      ccactcGCTGCTCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTtct                                                                             <  1:89192/80‑1 (MQ=255)                                                                                                                                      ccactcGCTGCTCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTtct                                                                             >  2:89192/4‑83 (MQ=255)                                                                                                                                        tctcGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGa         >  1:318303/1‑149 (MQ=255)                                                                                                                                            gCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACcg      <  6:169421/148‑1 (MQ=255)                                                                                                                                               gAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACcgcg    >  4:36969/1‑147 (MQ=255)                                                                                                                                                agcagcGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCa  <  6:221391/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    CGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACAT‑GTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCA  >  NC_002947/4858377‑4858666

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑