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.96.8% 127.2 / ‑2.1 31intergenic (+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 . (14/16);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

CCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAA  >  NC_000913/321425‑321659
                                                                                                                                     |                                                                                                     
ccGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTgggggg                                                                                                  <  1:422269/139‑1 (MQ=255)
 cGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTgggggga                                                                                                 <  2:157699/139‑2 (MQ=255)
                 aGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAAt                                                                                >  2:799119/1‑139 (MQ=255)
                              gACGTTTGGATGAAAGATTATCATTTGTCATAGAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                   <  2:316278/139‑1 (MQ=255)
                              gACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGCGGTATGTTGCTAGCCAATTAAAAAAGAACGCAATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                   >  2:70011/1‑139 (MQ=255)
                                   ttGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTg                                                                                                       >  2:580689/1‑99 (MQ=255)
                                   ttGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTg                                                                                                       <  1:580689/99‑1 (MQ=255)
                                                       tGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAAt                                          >  2:400681/1‑139 (MQ=255)
                                                        gTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                           <  2:602786/105‑1 (MQ=255)
                                                        gTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                           >  1:602786/1‑105 (MQ=255)
                                                           aTACAAATGAGGGGTGGTATGTTGCTAGCCAAGTAAAAAAGAACGCCATATGTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTACTTGTGATTGTCAATTCCt                                      <  2:487098/139‑1 (MQ=255)
                                                              cAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATg                                                                         >  2:89284/1‑101 (MQ=255)
                                                              cAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATg                                                                         <  1:89284/101‑1 (MQ=255)
                                                                  tGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                               <  1:400681/139‑1 (MQ=255)
                                                                      gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                               >  1:423950/1‑135 (MQ=255)
                                                                      gggTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                               <  2:423950/135‑1 (MQ=255)
                                                                         tggtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAAt                                                                                <  1:658033/83‑1 (MQ=255)
                                                                         tggtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAAt                                                                                >  2:658033/1‑83 (MQ=255)
                                                                         tggtATGTTGCTAGCCAATTAAAAAAGAACACCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                   <  1:337191/96‑1 (MQ=255)
                                                                         tggtATGTTGCTAGCCAATTAAAAAAGAACACCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                   >  2:337191/1‑96 (MQ=255)
                                                                           gtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTggggggaaa                                                                                               <  2:535462/67‑4 (MQ=255)
                                                                           gtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGaaa                                                                                               >  1:535462/1‑67 (MQ=255)
                                                                           gtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAcc                              >  2:309732/1‑131 (MQ=255)
                                                                           gtATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAcc                              <  1:309732/131‑1 (MQ=255)
                                                                                     aGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCa                                <  2:304581/119‑1 (MQ=255)
                                                                                     aGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCa                                >  1:304581/1‑119 (MQ=255)
                                                                                     aGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTg            <  1:325159/139‑1 (MQ=255)
                                                                                               aaaaGAACACCATTTTTATTGATGATTGATCACCCGGTGGGGGGAAACAGCGAGCATAAATGGGCATGAAGTATTGGCGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTCACGCCCTGAAATCAAGaa  <  2:754016/139‑1 (MQ=255)
                                                                                                                     tgattgatCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                           >  2:443291/1‑44 (MQ=255)
                                                                                                                     tgattgatCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCa                                                                           <  1:443291/44‑1 (MQ=255)
                                                                                                                           atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                   <  1:671786/46‑1 (MQ=255)
                                                                                                                           atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTaa                                                                   >  2:671786/1‑46 (MQ=255)
                                                                                                                           atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTAt                                                           >  1:257116/1‑54 (MQ=255)
                                                                                                                           atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTAt                                                           <  2:257116/54‑1 (MQ=255)
                                                                                                                           atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                               <  2:328987/82‑1 (MQ=255)
                                                                                                                           atCCCCCGGTGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCAc                               >  1:328987/1‑82 (MQ=255)
                                                                                                                                     |                                                                                                     
CCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGAATGTCAATTTCTTTGTCACCTGTATTGGTGACGCCCTGAAATCAAGAA  >  NC_000913/321425‑321659

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