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.97.7% 171.9 / ‑3.2 43intergenic (+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 . (20/22);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGAAAAAC  >  NC_000913/321422‑321692
                                                                                                                                        |                                                                                                                                      
aTGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTggg                                                                                                                                      <  2:941786/139‑1 (MQ=255)
        agaaaTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACa                                                                                                                             <  2:84708/139‑1 (MQ=255)
         gaaaTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACAt                                                                                                                            <  2:835150/139‑1 (MQ=255)
                              tttGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAg                                                                                                       >  2:244441/1‑139 (MQ=255)
                              tttGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAg                                                                                                       <  1:611645/139‑1 (MQ=255)
                               ttGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGt                                                                                                      >  2:1018510/1‑139 (MQ=255)
                                tGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTa                                                                                                     <  1:924094/139‑1 (MQ=255)
                                 gACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                                                    <  1:244441/139‑1 (MQ=255)
                                       tGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGt                                                                                              >  2:816096/1‑139 (MQ=255)
                                                               tACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCtt                                                                      >  1:440927/1‑139 (MQ=255)
                                                                    aTGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCa                                                                 <  1:816096/139‑1 (MQ=255)
                                                                     tGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                                <  1:1018510/139‑1 (MQ=255)
                                                                      gAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTg                                                                                                                                        >  2:673652/1‑67 (MQ=255)
                                                                      gAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTg                                                                                                                                        <  1:673652/67‑1 (MQ=255)
                                                                      gAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATc                                                                                                                           <  1:681735/79‑1 (MQ=255)
                                                                      gAGGGGTGGTATGTTGCTAGCAAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATc                                                                                                                           >  2:681735/1‑79 (MQ=255)
                                                                         gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGTTCCCCCGGTg                                                                                                                                        >  2:798087/1‑64 (MQ=255)
                                                                         gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTg                                                                                                                                        <  1:798087/64‑1 (MQ=255)
                                                                              gtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTgg                                                       >  2:1210311/1‑139 (MQ=255)
                                                                                   ttGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGt                                                                   >  1:236400/1‑122 (MQ=255)
                                                                                   ttGCTAGCCAATTAAAAAAGAACGCCATATTTAGTGATGATTGATCCCCCGGTGGGGGGAAACATCGTGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGt                                                                   <  2:236400/122‑1 (MQ=255)
                                                                                     gCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGTATAAATGGGCATg                                                                                                          <  1:628008/81‑1 (MQ=255)
                                                                                     gCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGTATAAATGGGCATg                                                                                                          >  2:628008/1‑81 (MQ=255)
                                                                                        aGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                                >  2:251322/1‑120 (MQ=255)
                                                                                        aGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                                <  1:251322/120‑1 (MQ=255)
                                                                                        aGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTg                                             <  1:1210311/139‑1 (MQ=255)
                                                                                            aaTTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                                <  1:794307/116‑1 (MQ=255)
                                                                                            aaTTAAAAAAGAACGCCACATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                                >  2:794307/1‑116 (MQ=255)
                                                                                              ttAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                                                            >  1:811470/1‑70 (MQ=255)
                                                                                              ttAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                                                            <  2:811470/70‑1 (MQ=255)
                                                                                              ttAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                                                                                                   <  2:108280/79‑1 (MQ=255)
                                                                                              ttAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                                                                                                   >  1:108280/1‑79 (MQ=255)
                                                                                                  aaaaGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGaa                                   <  2:1328057/139‑1 (MQ=255)
                                                                                                       aaCGCCATATTTATTGATGATTGATCCCCCGGTgggggga                                                                                                                                  <  2:739437/40‑2 (MQ=255)
                                                                                                       aaCGCCATATTTATTGATGATTGATCCCCCGGTgggggga                                                                                                                                  >  1:739437/1‑39 (MQ=255)
                                                                                                            cATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACg                            <  2:197281/136‑1 (MQ=255)
                                                                                                            cATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACg                            >  1:197281/1‑136 (MQ=255)
                                                                                                                   attGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATgg                                                                                                 >  1:1115535/1‑60 (MQ=255)
                                                                                                                   attGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATgg                                                                                                 <  2:1115535/60‑1 (MQ=255)
                                                                                                                      gatgatTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATg                                                                                                          <  2:532182/48‑1 (MQ=255)
                                                                                                                      gatgatTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATg                                                                                                          >  1:532182/1‑48 (MQ=255)
                                                                                                                          attgatCCCCCGGTGGGGGGAAACATCGGGCATaaa                                                                                                                  <  2:1211102/36‑1 (MQ=255)
                                                                                                                          attgatCCCCCGGTGGGGGGAAACATCGGGCATaaa                                                                                                                  >  1:1211102/1‑36 (MQ=255)
                                                                                                                          attgatCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                                                            <  2:67855/42‑1 (MQ=255)
                                                                                                                          attgatCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                                                            >  1:67855/1‑42 (MQ=255)
                                                                                                                          attgatCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTAc           >  2:707689/1‑139 (MQ=255)
                                                                                                                              atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTc                                                                  >  2:1175162/1‑80 (MQ=255)
                                                                                                                              atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTc                                                                  <  1:1175162/80‑1 (MQ=255)
                                                                                                                                 ccccGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGTGACTCCGTGCTGCTACTGGaaaa    >  1:1022851/1‑139 (MQ=255)
                                                                                                                                   ccGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGAAAAAc  >  1:84956/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGAAAAAC  >  NC_000913/321422‑321692

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