Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 321,485 T→G 9.4% intergenic (+404/‑123) 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,4850TG 9.4% 174.3 / 10.8 53intergenic (+404/‑123)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 T (24/24);  new base G (3/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GCCATTTTTAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGA  >  NC_000913/321359‑321611
                                                                                                                              |                                                                                                                              
gCCATTTTTAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGtt                                                                                                           >  4:236779/1‑148 (MQ=255)
 ccATTTTTAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAACTTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGggt                                                                                                                   >  4:911402/1‑139 (MQ=255)
 ccATTTTTAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGggt                                                                                                                   <  3:911402/139‑1 (MQ=255)
                 ttCTTATGCCAAAAATGTTCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATa                                                                                                                               >  3:796567/1‑111 (MQ=255)
                 ttCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATa                                                                                                                               <  4:796567/111‑1 (MQ=255)
                     tATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGCCATACAAATg                                                                                                                         <  4:358118/113‑1 (MQ=255)
                     tATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGCCATACAAATg                                                                                                                         >  3:358118/1‑113 (MQ=255)
                          cAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATg                                                                                                                         <  2:953550/108‑1 (MQ=255)
                          cAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATg                                                                                                                         >  1:953550/1‑108 (MQ=255)
                               aTGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCAtattta                                                                            <  3:236779/148‑1 (MQ=255)
                                          tttCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGggtg                                                                                                                  <  2:857392/99‑1 (MQ=255)
                                          tttCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGggtg                                                                                                                  >  1:857392/1‑99 (MQ=255)
                                                             cGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAAt                                                                                                 <  4:165299/97‑1 (MQ=255)
                                                             cGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAAt                                                                                                 >  3:165299/1‑97 (MQ=255)
                                                                          aaTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGAtgattgat                                                                <  2:174326/117‑1 (MQ=255)
                                                                          aaTAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGAtgattgat                                                                >  1:174326/1‑117 (MQ=255)
                                                                            tAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGc                                                                                    >  4:120250/1‑95 (MQ=255)
                                                                            tAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGc                                                                                    <  3:120250/95‑1 (MQ=255)
                                                                            tAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGAtgattg                                                                  >  3:912324/1‑113 (MQ=255)
                                                                            tAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGAtgattg                                                                  <  4:912324/113‑1 (MQ=255)
                                                                             aaaaTTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATgg                              >  3:1009146/1‑148 (MQ=255)
                                                                                aTTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTggggg                                                   >  2:946837/1‑124 (MQ=255)
                                                                                aTTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTggggg                                                   <  1:946837/124‑1 (MQ=255)
                                                                                  tAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTggg                                                     <  4:808303/120‑1 (MQ=255)
                                                                                  tAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTggg                                                     >  3:808303/1‑120 (MQ=255)
                                                                                      tATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCAtattta                                                                            >  4:322034/1‑93 (MQ=255)
                                                                                      tATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCAtattta                                                                            <  3:322034/93‑1 (MQ=255)
                                                                                         tGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATccc                                                             <  1:418828/105‑1 (MQ=255)
                                                                                         tGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATccc                                                             >  2:418828/1‑105 (MQ=255)
                                                                                         tGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGa                        <  2:56876/142‑1 (MQ=255)
                                                                                         tGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGa                        >  1:56876/1‑142 (MQ=255)
                                                                                         tGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAAt                  <  4:1009146/148‑1 (MQ=255)
                                                                                                     ttGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGt             >  2:677353/1‑141 (MQ=255)
                                                                                                     ttGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGt             <  1:677353/141‑1 (MQ=255)
                                                                                                           gAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATg                               <  3:659663/117‑1 (MQ=255)
                                                                                                           gAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATg                               >  4:659663/1‑117 (MQ=255)
                                                                                                            aaaGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAg        <  3:1107737/139‑1 (MQ=255)
                                                                                                            aaaGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAg        >  4:1107737/1‑139 (MQ=255)
                                                                                                                 ttATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATtg                                                                         <  2:475563/69‑1 (MQ=255)
                                                                                                                 ttATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATtg                                                                         >  1:475563/1‑69 (MQ=255)
                                                                                                                  cctCATTTGTATGACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCg                                         >  1:777999/3‑100 (MQ=255)
                                                                                                                  cctCATTTGTATGACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCg                                         <  2:777999/98‑1 (MQ=255)
                                                                                                                    tCATTTGTATGACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACg                                                                                     >  1:649488/1‑54 (MQ=37)
                                                                                                                    tCATTTGTATGACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACg                                                                                     <  2:649488/54‑1 (MQ=37)
                                                                                                                    tCATTTGTATGACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGa  >  2:334244/1‑137 (MQ=255)
                                                                                                                     cATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAAt                                                                                                 <  3:798072/41‑1 (MQ=255)
                                                                                                                     cATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAAt                                                                                                 >  4:798072/1‑41 (MQ=255)
                                                                                                                     cATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATcc                                                              >  3:330958/1‑76 (MQ=255)
                                                                                                                     cATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATcc                                                              <  4:330958/76‑1 (MQ=255)
                                                                                                                     cATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATcc                                                              <  4:416379/76‑1 (MQ=255)
                                                                                                                     cATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATcc                                                              >  3:416379/1‑76 (MQ=255)
                                                                                                                       tttGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATg                               >  1:796484/1‑105 (MQ=255)
                                                                                                                       tttGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATg                               <  2:796484/105‑1 (MQ=255)
                                                                                                                              tACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGt                     >  4:616380/1‑108 (MQ=255)
                                                                                                                              tACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGt                     <  3:616380/108‑1 (MQ=255)
                                                                                                                              |                                                                                                                              
GCCATTTTTAAATGGATTTCTTATGCCAAAAATGATCGTCGATTTCCCCATAAAATGTGAGCGATGCCGAAAGAAATAAAATTAGTTATTGCATTTGACGTTTGGATGAAAGATTATCATTTGTCATACAAATGAGGGGTGGTATGTTGCTAGCCAATTAAAAAAGAACGCCATATTTATTGATGATTGATCCCCCGGTGGGGGGGAAACATCGGGCATAAATGGGCATGAAGTAATGGAGTATTAGTTGTGA  >  NC_000913/321359‑321611

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: