Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 321,564 (G)7→6 intergenic (+483/‑44) rclR → / → ykgE reactive chlorine species (RCS)‑specific activator of the rcl genes/cysteine‑rich LutA family protein, putative electron transport chain YkgEFG component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913321,5580G.98.4% 250.4 / ‑3.1 62intergenic (+477/‑50)rclR/ykgEreactive chlorine species (RCS)‑specific activator of the rcl genes/cysteine‑rich LutA family protein, putative electron transport chain YkgEFG component
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (31/30);  total (32/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 = 9.74e-01

CGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGAAAA  >  NC_000913/321420‑321690
                                                                                                                                          |                                                                                                                                    
cGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTg                                                                                                                                      <  2:1608935/139‑1 (MQ=255)
cGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTg                                                                                                                                      <  2:69985/139‑1 (MQ=255)
  aTGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTggg                                                                                                                                    >  2:131678/1‑139 (MQ=255)
  aTGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTggg                                                                                                                                    >  1:207284/1‑139 (MQ=255)
        aaagaaaTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGaaa                                                                                                                              >  2:164684/1‑139 (MQ=255)
          agaaaTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACa                                                                                                                           <  2:1259131/139‑1 (MQ=255)
             aaTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCg                                                                                                                        >  1:305224/1‑139 (MQ=255)
               tAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCggg                                                                                                                      >  2:347724/1‑139 (MQ=255)
                aaaaTTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGc                                                                                                                     <  1:131678/139‑1 (MQ=255)
                             gCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATg                                                                                                        <  1:164684/139‑1 (MQ=255)
                                        ttGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAg                                                                                             >  1:1393377/1‑139 (MQ=255)
                                         tGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATg                                                                                                <  2:240951/135‑1 (MQ=255)
                                         tGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATg                                                                                                >  1:240951/1‑135 (MQ=255)
                                           gATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTAt                                                                                          <  2:305224/139‑1 (MQ=255)
                                                            tGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAAt                                                                         >  2:543495/1‑139 (MQ=255)
                                                             gTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATg                                                                                                        <  1:174051/107‑1 (MQ=255)
                                                             gTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATg                                                                                                        >  2:174051/1‑107 (MQ=255)
                                                                 tACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACAt                                                                                                                          <  2:1036853/85‑1 (MQ=255)
                                                                 tACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCAACCGGT‑GGGGGGAAACAt                                                                                                                          >  1:1036853/1‑85 (MQ=255)
                                                                    aaaTGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGt                                                                 <  1:347724/139‑1 (MQ=255)
                                                                           gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACAt                                                                                                                          <  1:1290630/75‑1 (MQ=255)
                                                                           gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACAt                                                                                                                          >  2:1290630/1‑75 (MQ=255)
                                                                           gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTAt                                                                                          <  1:513421/107‑1 (MQ=255)
                                                                           gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTAt                                                                                          >  2:513421/1‑107 (MQ=255)
                                                                           gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGa                                                                                 >  2:163592/1‑116 (MQ=255)
                                                                           gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGa                                                                                 <  1:163592/116‑1 (MQ=255)
                                                                                gtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GTGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTAAATTCCTTTGTCAc                                                              >  1:1457065/1‑130 (MQ=255)
                                                                                gtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTc                                                                >  1:1532025/1‑128 (MQ=255)
                                                                                gtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTc                                                                <  2:1532025/128‑1 (MQ=255)
                                                                                gtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                              <  2:1457065/130‑1 (MQ=255)
                                                                                   tGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATgg                                                                                               <  2:1612678/94‑1 (MQ=255)
                                                                                   tGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATgg                                                                                               >  1:1612678/1‑94 (MQ=255)
                                                                                         tAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCt                                            <  2:207284/139‑1 (MQ=255)
                                                                                           gCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCa                                                                                                          <  1:56711/75‑1 (MQ=255)
                                                                                           gCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCa                                                                                                          >  2:56711/1‑75 (MQ=255)
                                                                                           gCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                                                                                                 >  2:567172/1‑84 (MQ=255)
                                                                                           gCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                                                                                                 <  1:567172/84‑1 (MQ=255)
                                                                                           gCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                              <  1:524914/119‑1 (MQ=255)
                                                                                           gCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                              >  2:524914/1‑119 (MQ=255)
                                                                                           gCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGa                                          >  1:477150/1‑139 (MQ=255)
                                                                                                       aGAACGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATgg                              <  2:477150/139‑1 (MQ=255)
                                                                                                         aaCGCCATATTTATTGATGATTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCa                            <  2:1393377/139‑1 (MQ=255)
                                                                                                                       tgatgaTTGATCCCCCGGT‑GGGGGGAAACATCGGGCATaaa                                                                                                                >  1:788904/1‑41 (MQ=255)
                                                                                                                       tgatgaTTGATCCCCCGGT‑GGGGGGAAACATCGGGCATaaa                                                                                                                <  2:788904/41‑1 (MQ=255)
                                                                                                                        gatgatTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATg                                                                                                        >  1:724327/1‑48 (MQ=255)
                                                                                                                        gatgatTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATg                                                                                                        <  2:724327/48‑1 (MQ=255)
                                                                                                                        gatgatTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGt                                                                 >  1:448950/1‑87 (MQ=255)
                                                                                                                        gatgatTGATCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGt                                                                 <  2:448950/87‑1 (MQ=255)
                                                                                                                          tgattgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTgctgct           >  1:1512676/1‑139 (MQ=255)
                                                                                                                          tgattgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGAAGCCCTGAAATCAAGAATGGCACGAGAc                      <  1:960789/128‑1 (MQ=255)
                                                                                                                          tgattgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGAATTGGTGACGCCCTGAAATCAAGAATGGCACGAGAc                      >  2:960789/1‑128 (MQ=255)
                                                                                                                           gattgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATg                                                                                                >  2:1531624/1‑53 (MQ=255)
                                                                                                                           gattgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATg                                                                                                <  1:1531624/53‑1 (MQ=255)
                                                                                                                            attgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGa                                                                                 >  2:1424735/1‑67 (MQ=255)
                                                                                                                            attgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGa                                                                                 <  1:1424735/67‑1 (MQ=255)
                                                                                                                            attgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACt                     <  2:1427525/127‑1 (MQ=255)
                                                                                                                            attgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACt                     >  1:1427525/1‑127 (MQ=255)
                                                                                                                              tgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTg                                           >  2:1150041/1‑103 (MQ=255)
                                                                                                                              tgatCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTg                                           <  1:1150041/103‑1 (MQ=255)
                                                                                                                                 tCCCCCGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGaa    >  1:1415436/1‑139 (MQ=255)
                                                                                                                                   ccccGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGaaaa  >  2:687426/1‑139 (MQ=255)
                                                                                                                                   ccccGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGaaaa  >  1:1428115/1‑139 (MQ=255)
                                                                                                                                     ccGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                                                  <  1:93329/41‑1 (MQ=255)
                                                                                                                                     ccGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                                                  >  2:93329/1‑41 (MQ=255)
                                                                                                                                     ccGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGAc                                                 <  2:372019/90‑1 (MQ=255)
                                                                                                                                     ccGGT‑GGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGAc                                                 >  1:372019/1‑90 (MQ=255)
                                                                                                                                          |                                                                                                                                    
CGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGAAAA  >  NC_000913/321420‑321690

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