Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,450,074	T→C	9.1%	intergenic (‑7/+23)	PP_1268 ← / ← PP_1269	HtrA suppressor protein SohA/universal stress protein family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,450,074	0	T	C	9.1%	144.8 / 5.6	55	intergenic (‑7/+23)	PP_1268/PP_1269	HtrA suppressor protein SohA/universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base T (31/19);  new base C (2/3);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

GGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAGTTCGCAGTT  >  NC_002947/1449927‑1450219                                                                                                                                                    |                                                                                                                                                   ggCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTg                                                                                                                                                   <  3:240447/149‑1 (MQ=255)  gcgcgcCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGc                                                                                                                                                  <  2:212222/149‑1 (MQ=255)          cACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCAc                                                                                                                                          >  5:388367/1‑149 (MQ=255)           aCCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACg                                                                                                                                         <  1:129209/149‑1 (MQ=255)             cGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGt                                                                                                                                       >  1:382373/1‑149 (MQ=255)              gTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGt                                                                                                                                       >  5:359048/1‑148 (MQ=255)                 tCGGGAACATTGGATTTATAGCGATCTGTCATAGTTGACGCGACGTAGATGGTTGAAGCTATGAGTTATGTCCTGTGTGTTAGAGCGATACAACTCAAGGTAATTCGTTTGCTTTGCATTGTCATTTGAGTTGCATTACCACgtttat                                                                                                                                    <  4:125996/148‑6 (MQ=255)                    ggAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGc                                                                                                                               <  4:156394/149‑1 (MQ=255)                         aGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGt                                                                                                                            >  5:223434/1‑147 (MQ=255)                           tGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCg                                                                                                                              >  4:292275/1‑143 (MQ=255)                           tGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCg                                                                                                                              <  3:292275/143‑1 (MQ=255)                                       gATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTAcctt                                                                                                                                          <  2:15597/120‑3 (MQ=21)                                       gATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTAcctt                                                                                                                                          <  2:8379/120‑3 (MQ=21)                                        aTCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGc                                                                                                                               >  8:224459/1‑129 (MQ=255)                                        aTCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCaa                                                                                                           >  2:314314/1‑149 (MQ=255)                                          cGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAaca                                                                                                         >  4:400426/1‑149 (MQ=255)                                            gTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGAAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCACACCAACACt                                                                                                       >  7:13718/1‑149 (MQ=255)                                              cAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTgg                                                                                                     >  3:23706/1‑149 (MQ=255)                                                 cGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCAt                                                                                                  <  7:78136/149‑1 (MQ=255)                                                   tGGACTCGACGTCGAAGGTTGCAGCCATGAGTAATGTGATGTGGGTGACAGCGATACAACTCATGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATT‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGt                                                                                                <  4:23706/149‑1 (MQ=255)                                                    ggACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGt                                                                                                >  3:232469/1‑148 (MQ=255)                                                     gACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTc                                                                                               >  2:93639/1‑148 (MQ=255)                                                        tCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTAcc                                                                                           >  2:165935/1‑149 (MQ=255)                                                          gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTg                                                                                                      <  7:102723/136‑1 (MQ=255)                                                          gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTg                                                                                                      >  8:102723/1‑136 (MQ=255)                                                                gAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGCAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCaa                                                                                    <  1:314314/149‑1 (MQ=255)                                                                 aaGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAAt                                                                                   >  1:174865/1‑148 (MQ=255)                                                                       gCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTc                                                                            <  1:93639/149‑1 (MQ=255)                                                                                  tAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGa                                                                 >  5:220332/1‑149 (MQ=255)                                                                                    aTGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGAc                                                                <  6:223434/148‑1 (MQ=255)                                                                                    aTGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGAc                                                                <  6:388367/148‑1 (MQ=255)                                                                                    aTGTCCTGAGGGTGACAGCGATACAACTAAAGGAAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGAAG‑AGCGGTCAGGCAGAGCACACCAACACTGGCATGTCTACCCGGTAAATCACCTTAGCTCTGATTGACg                                                               >  1:85687/1‑149 (MQ=255)                                                                                       tCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACggg                                                             <  4:232469/148‑1 (MQ=255)                                                                                        ccTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTc                                                           >  7:113977/1‑149 (MQ=255)                                                                                            tgGGTGTCAGCGATACAACTGATGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGAATTACCACGGTGATG‑AGCGTTCAGTCAGAGCCCACCAACACTTGGATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAg                                                        <  2:85687/148‑1 (MQ=255)                                                                                            tgGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCa                                                         <  4:174561/147‑1 (MQ=255)                                                                                              ggTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCaa                                                     <  6:220332/149‑1 (MQ=255)                                                                                              ggTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGACAATCACCTTCGCGCTGATTGACGGGTCCAGCaa                                                     <  1:58661/149‑1 (MQ=255)                                                                                                       gATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGccacct                                            >  8:264952/1‑145 (MQ=255)                                                                                                         tACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAAcc                                          >  5:123236/1‑149 (MQ=255)                                                                                                               ccgtggTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGa                                      >  5:171265/5‑147 (MQ=25)                                                                                                               ccgtggTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGa                                      <  6:171265/143‑1 (MQ=25)                                                                                                                cgtggTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAg                                   >  5:363043/4‑149 (MQ=37)                                                                                                                 aaGGCAATGCAATCGCATTACCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAgg                                  >  5:310740/1‑149 (MQ=255)                                                                                                                         gCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg                          >  3:281601/1‑149 (MQ=255)                                                                                                                         gCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg                          >  7:122817/1‑149 (MQ=255)                                                                                                                         gCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg                          >  2:318044/1‑149 (MQ=255)                                                                                                                           gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtggc                        >  4:16380/1‑149 (MQ=255)                                                                                                                           gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtggc                        >  5:395624/1‑149 (MQ=255)                                                                                                                                cattgcCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCat                    >  3:158580/1‑148 (MQ=255)                                                                                                                                 attgcCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatga                  >  3:195330/1‑149 (MQ=255)                                                                                                                                 attgcCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatga                  <  1:318044/149‑1 (MQ=255)                                                                                                                                            aaTGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCTACCGCGACAGGTGTCGGTGCCGCATGAAGATAAGTTCg       >  3:55742/1‑149 (MQ=255)                                                                                                                                            aaTGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAGTTCg       >  4:327199/1‑149 (MQ=255)                                                                                                                                                 aGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAGTTCGCAGtt  >  6:201409/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGATCAGTTCGCAGTT  >  NC_002947/1449927‑1450219

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