Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,450,049	G→A	11.3%	R7C (CGC→TGC)	PP_1268 ←	HtrA suppressor protein SohA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,450,049	0	G	A	11.3%	139.3 / 7.0	53	R7C (CGC→TGC)	PP_1268	HtrA suppressor protein SohA
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base A (3/3);  total (29/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 = 2.42e-01

GTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGC  >  NC_002947/1449906‑1450197                                                                                                                                                |                                                                                                                                                      gTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATTAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGatt                                                                                                                                                    <  1:100617/147‑1 (MQ=255) gTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGat                                                                                                                                                     <  1:432314/146‑1 (MQ=255)  tCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAATAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGattgca                                                                                                                                                 >  3:174561/1‑149 (MQ=255)    aCGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGattgcatt                                                                                                                                               >  7:345338/1‑149 (MQ=255)      gTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGattgcattgc                                                                                                                                             <  2:286136/149‑1 (MQ=255)                   cAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAg                                                                                                                                >  2:58661/1‑149 (MQ=255)                      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:388367/148‑1 (MQ=255)                                                                                                         aTGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGAc                                          <  6:223434/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                <  6:171265/143‑1 (MQ=25)                                                                                                                                    ccgtggTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGa                >  5:171265/5‑147 (MQ=25)                                                                                                                                     cgtggTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAg             >  5:363043/4‑149 (MQ=37)                                                                                                                                      aaGGCAATGCAATCGCATTACCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAgg            >  5:310740/1‑149 (MQ=255)                                                                                                                                              gCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg    >  7:122817/1‑149 (MQ=255)                                                                                                                                              gCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg    >  3:281601/1‑149 (MQ=255)                                                                                                                                              gCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg    >  2:318044/1‑149 (MQ=255)                                                                                                                                                gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtggc  >  5:395624/1‑149 (MQ=255)                                                                                                                                                gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtggc  >  4:16380/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      GTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGC  >  NC_002947/1449906‑1450197

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