Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,696,568 C→T 11.9% intergenic (+20/‑192) pyrC' → / → hslV dihydroorotase‑like protein/ATP‑dependent HslVU protease peptidase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,696,5680CT11.9% 128.4 / 7.0 51intergenic (+20/‑192)pyrC'/hslVdihydroorotase‑like protein/ATP‑dependent HslVU protease peptidase subunit
Reads supporting (aligned to +/- strand):  ref base C (19/25);  new base T (2/4);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

CGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATAT  >  NC_002947/5696421‑5696713
                                                                                                                                                   |                                                                                                                                                  
cGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACg                                                                                                                                                   >  7:123891/1‑149 (MQ=255)
 gCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGc                                                                                                                                                  <  8:246918/149‑1 (MQ=255)
           aCCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCca                                                                                                                                         >  1:337061/1‑148 (MQ=255)
            ccGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTc                                                                                                                                            <  4:278302/144‑1 (MQ=255)
            ccGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTc                                                                                                                                            >  3:278302/1‑144 (MQ=255)
             cGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAg                                                                                                                                      >  5:354693/1‑149 (MQ=255)
                  cTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATca                                                                                                                                 <  6:21374/149‑1 (MQ=255)
                   tGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcac                                                                                                                                >  4:69260/1‑149 (MQ=255)
                     gCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcacca                                                                                                                              >  6:94356/1‑149 (MQ=255)
                     gCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcacca                                                                                                                              <  8:166456/149‑1 (MQ=255)
                      cGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTtacca                                                                                                                              >  2:297356/1‑148 (MQ=255)
                      cGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTtacca                                                                                                                              >  2:293240/1‑148 (MQ=255)
                      cGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcaccac                                                                                                                             <  2:286235/149‑1 (MQ=255)
                       gAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcaccac                                                                                                                             >  6:13655/1‑148 (MQ=255)
                          cAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAAtt                                                                                                                         <  2:337061/149‑1 (MQ=255)
                           aaTGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTc                                                                                                                        >  4:82775/1‑149 (MQ=255)
                           aaTGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACAGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTc                                                                                                                        >  7:92870/1‑149 (MQ=255)
                               gctctcCCGGGGCGAGAACTGCCCGTTCCTCGGCCACTGCCTGCCGGGCGCGGTGCGATATACCTTTGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAg                                                                                                                    <  8:288730/147‑1 (MQ=255)
                                  ctcCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGc                                                                                                                  <  3:82775/148‑1 (MQ=255)
                                     ccGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑tgt                                                                                                             >  4:313495/1‑149 (MQ=255)
                                     ccGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑tg                                                                                                              <  3:163901/148‑1 (MQ=255)
                                     ccGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑agt                                                                                                             >  2:208206/1‑149 (MQ=255)
                                       gTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTgg                                                                                                           >  3:181830/1‑149 (MQ=255)
                                          gCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTaa                                                                                                        <  5:94356/149‑1 (MQ=255)
                                                aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCtgt                                                                                                  <  2:274228/149‑1 (MQ=255)
                                                aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCtgt                                                                                                  <  8:123891/149‑1 (MQ=255)
                                                 cTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑tg                                                                                                              <  2:258278/136‑1 (MQ=255)
                                                 cTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑tg                                                                                                              >  1:258278/1‑136 (MQ=255)
                                                   gCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTggg                                                                                               >  7:340203/1‑149 (MQ=255)
                                                    cccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCAC‑ATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                              <  8:36928/149‑1 (MQ=21)
                                                    cccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCAC‑ATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                              <  1:293240/149‑1 (MQ=21)
                                                    cccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGCCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCAC‑ATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                              <  1:297356/149‑1 (MQ=21)
                                                     ccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCAC‑ATGCCTTGGAATTGTGGTGATCCTGTGGGAg                                                                                             <  1:292988/149‑1 (MQ=21)
                                                         tCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCggg                                                                                         >  8:156949/1‑149 (MQ=255)
                                                               gCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGcc                                                                                   >  2:17954/1‑149 (MQ=255)
                                                                ccACTGCCTGCCTTGCGCGGTGAGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGccc                                                                                  <  6:354693/149‑1 (MQ=255)
                                                                      cctgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAAt                                                                            <  7:156949/149‑1 (MQ=255)
                                                                        tgccGGGCGCGGTGCGTTATACCTAGGTCGATGGCCACATCTGCCACCAGGCCTGAGTACCATGCTTAGCGGGCACGCACGCTCCCACACGAACACCACAAATCCACGGCA‑TGTTGTACCCCTGTGGGAGCGGGCATGCCCCCGAAAGc                                                                          >  8:171130/1‑149 (MQ=255)
                                                                           cGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCtt                                                                       >  1:105528/1‑149 (MQ=255)
                                                                              gcgcGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTcac                                                                    <  8:340203/149‑1 (MQ=255)
                                                                                                        ggCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCAt                                          >  1:237204/1‑149 (MQ=255)
                                                                                                          ccACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTc                                        <  2:237204/149‑1 (MQ=255)
                                                                                                             caTCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATc                                     <  3:69260/149‑1 (MQ=255)
                                                                                                                        aGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCAACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTg                          <  5:4801/149‑1 (MQ=255)
                                                                                                                                gTACCTTGCTTAGGGGCCTCGGGCTCTGCCACAGGACCACCACATTTGCAACGCA‑TCTGGTAACCCTGGGGTAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCtt                  <  3:144030/149‑1 (MQ=255)
                                                                                                                                 tACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCtt                  <  2:105528/148‑1 (MQ=255)
                                                                                                                                             cGGGAAGGCCCGCTCCAACAGTAGCACAACAATTCCAAGTCC‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccc      <  1:301778/148‑1 (MQ=255)
                                                                                                                                              gggCACGCCCGCTCCCACACGATCACCACAATTCGCAGGCA‑TGTGGATACCCAGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCTCCCAGCGCGACCATTCATCGCCCCCGGTTGAAATCCTTCCCCCCACCCCCat    <  8:92870/149‑1 (MQ=255)
                                                                                                                                              gggCACGCCAGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCat    <  7:171130/149‑1 (MQ=255)
                                                                                                                                                gCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCatat  <  5:13655/149‑1 (MQ=255)
                                                                                                                                                gCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCatat  <  3:278291/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCAAGCGTCCACCGTTGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCA‑TGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATAT  >  NC_002947/5696421‑5696713

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