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.2% 225.1 / ‑2.2 56intergenic (+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 (0/0);  major base . (28/27);  minor base T (1/0);  total (29/27)
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.84e-01

ATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGAAAA  >  NC_000913/321422‑321690
                                                                                                                                         |                                                                                                                                    
aTGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTggg                                                                                                                                    >  1:1139509/1‑139 (MQ=255)
           aaTAAAATTAGTTATTGCATTTGATGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCg                                                                                                                        >  2:531905/1‑139 (MQ=255)
                         ttGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                                                          <  2:1139509/139‑1 (MQ=255)
                                tGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTa                                                                                                   >  2:900226/1‑139 (MQ=255)
                                  aCGGTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                                                                                                 <  2:903248/139‑1 (MQ=255)
                                    gTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATgg                                                                                               <  2:565041/139‑1 (MQ=255)
                                      ttGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACAt                                                                                                                          >  1:703846/1‑110 (MQ=255)
                                      ttGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACAt                                                                                                                          <  2:703846/110‑1 (MQ=255)
                                       tGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATaaa                                                                                                                >  2:1052966/1‑119 (MQ=255)
                                       tGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATaaa                                                                                                                <  1:1052966/119‑1 (MQ=255)
                                                         ttGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCaa                                                                          <  1:900226/139‑1 (MQ=255)
                                                           gTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACa                                                                                                                           >  1:1153209/1‑88 (MQ=255)
                                                           gTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACa                                                                                                                           <  2:1153209/88‑1 (MQ=255)
                                                                   aaTGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTc                                                                <  2:512307/139‑1 (MQ=255)
                                                                   aaTGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTc                                                                >  1:721758/1‑139 (MQ=255)
                                                                   aaTGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTc                                                                <  1:622885/139‑1 (MQ=255)
                                                                      gAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAcc                                                             <  1:531905/139‑1 (MQ=255)
                                                                      gAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAcc                                                             <  2:721758/139‑1 (MQ=255)
                                                                       aGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAcc                                                             >  1:1090533/1‑139 (MQ=255)
                                                                         gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                                                  >  2:317071/1‑99 (MQ=255)
                                                                         gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                                                  <  1:317071/99‑1 (MQ=255)
                                                                         gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGt                                                          >  1:1209539/1‑139 (MQ=255)
                                                                           gtggtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATg                                                                                                        <  1:1031415/91‑1 (MQ=255)
                                                                           gtggtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATg                                                                                                        >  2:1031415/1‑91 (MQ=255)
                                                                            tggtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                                                                                                 >  2:1003043/1‑97 (MQ=255)
                                                                            tggtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                                                                                                 <  1:1003043/97‑1 (MQ=255)
                                                                                    tGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGc                                               >  1:370353/1‑139 (MQ=255)
                                                                                     gCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGcc                                              >  2:882183/1‑139 (MQ=255)
                                                                                        aGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                                                          >  2:550748/1‑76 (MQ=255)
                                                                                        aGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                                                          <  1:550748/76‑1 (MQ=255)
                                                                                         gCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAcc                                                             >  1:1115200/1‑120 (MQ=255)
                                                                                         gCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAcc                                                             <  2:1115200/120‑1 (MQ=255)
                                                                                          ccAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGaa                                         >  2:6946/1‑139 (MQ=255)
                                                                                             aTTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                                                          <  2:65750/71‑1 (MQ=255)
                                                                                             aTTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                                                          >  1:65750/1‑71 (MQ=255)
                                                                                              ttAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCa                                     >  2:701330/1‑139 (MQ=255)
                                                                                              ttAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCa                                     >  2:1168044/1‑139 (MQ=255)
                                                                                                aaaaaaGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAg                                   <  1:6946/139‑1 (MQ=255)
                                                                                                 aaaaaGAACGCCATATTTATTGATGATTGATCCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAg                                   <  2:1090533/139‑1 (MQ=255)
                                                                                                       aaCGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCa                            <  2:1209539/139‑1 (MQ=255)
                                                                                                         cGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                              <  2:1193608/103‑1 (MQ=255)
                                                                                                         cGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                                                              >  1:1193608/1‑103 (MQ=255)
                                                                                                            cATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTc                                                                >  2:1012501/1‑98 (MQ=255)
                                                                                                            cATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTc                                                                <  1:1012501/98‑1 (MQ=255)
                                                                                                                    ttGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAg                                                                                                     <  2:148379/53‑1 (MQ=255)
                                                                                                                    ttGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAg                                                                                                     >  1:148379/1‑53 (MQ=255)
                                                                                                                     tgatgaTTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCa                                     >  1:804347/1‑116 (MQ=255)
                                                                                                                     tgatgaTTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCa                                     <  2:804347/116‑1 (MQ=255)
                                                                                                                      gatgatTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                                                                                                 >  2:1099757/1‑55 (MQ=255)
                                                                                                                      gatgatTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                                                                                                 <  1:1099757/55‑1 (MQ=255)
                                                                                                                       atgatTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGtgctgc            <  1:579207/139‑1 (MQ=255)
                                                                                                                          attgatCGCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTGGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTAc         >  1:797575/1‑139 (MQ=255)
                                                                                                                          attgatCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTAc         >  2:482331/1‑139 (MQ=255)
                                                                                                                              atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACg                                                >  1:744849/1‑96 (MQ=255)
                                                                                                                              atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACg                                                <  2:744849/96‑1 (MQ=255)
                                                                                                                              atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGa     >  2:250682/1‑139 (MQ=255)
                                                                                                                                  ccccGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGaaaa  >  1:1217918/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                    
ATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAATGGCACGAGACTCCGTGCTGCTACTGGAAAA  >  NC_000913/321422‑321690

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