Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 702,967 A→G intergenic (+32/+184) ECB_00650 → / ← kdpE conserved hypothetical protein/DNA‑binding response regulator in two‑component regulatory system with KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606702,9670AG98.2% 219.3 / NA 54intergenic (+32/+184)ECB_00650/kdpEconserved hypothetical protein/DNA‑binding response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (27/26);  minor base T (0/1);  total (27/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.99e-01

TAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGTTTTTTAC  >  REL606/702896‑703039
                                                                       |                                                                        
tAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATa                                                                      <  2:317420/76‑1 (MQ=255)
     cGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGtt                                                                 <  2:492181/76‑1 (MQ=255)
          gTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACt                                                            <  1:617379/76‑1 (MQ=255)
                   cATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTg                                                   <  2:808724/76‑1 (MQ=255)
                    atttaatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGa                                                  <  2:452994/76‑1 (MQ=255)
                    atttaatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGa                                                  <  2:1115911/76‑1 (MQ=255)
                      ttaatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATc                                                >  1:1059157/1‑76 (MQ=255)
                       taatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCa                                               >  1:373523/1‑76 (MQ=255)
                       taatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCa                                               <  1:44332/76‑1 (MQ=255)
                       taatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCa                                               >  1:1069076/1‑76 (MQ=255)
                        aatttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCaa                                              <  2:1121355/76‑1 (MQ=255)
                         atttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAAt                                             <  2:1003182/76‑1 (MQ=255)
                         atttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAAt                                             >  2:801267/1‑76 (MQ=255)
                          tttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATa                                            >  2:1007999/1‑76 (MQ=255)
                           ttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGTAATACGGTTGAACTGTGAAGGTGATCAATaa                                           <  1:177981/76‑1 (MQ=255)
                           ttaCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATaa                                           <  1:783829/76‑1 (MQ=255)
                            taCCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATaaa                                          <  2:53046/76‑1 (MQ=255)
                              ccTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATaaaaa                                        <  1:832251/76‑1 (MQ=255)
                                 ttCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATg                                     >  1:1121750/1‑76 (MQ=255)
                                 ttCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATg                                     >  2:57661/1‑76 (MQ=255)
                                  tCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGa                                    >  1:361916/1‑76 (MQ=255)
                                    cTGATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATc                                  >  1:697550/1‑76 (MQ=255)
                                      gATGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCaa                                >  1:115643/1‑76 (MQ=255)
                                       aTGTGAATGAAGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAAt                               >  1:915791/1‑76 (MQ=255)
                                          tgaatgaaGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCtt                            <  2:166765/76‑1 (MQ=255)
                                           gaatgaaGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTa                           >  1:657783/1‑76 (MQ=255)
                                           gaatgaaGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTa                           >  1:611514/1‑76 (MQ=255)
                                              tgaaGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTCatt                        >  2:956854/1‑76 (MQ=255)
                                              tgaaGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAatt                        >  2:956856/1‑76 (MQ=255)
                                               gaaGGTGAAGATATTCTATCCGGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAattt                       <  2:690609/76‑1 (MQ=255)
                                               gaaGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAattt                       >  1:775031/1‑76 (MQ=255)
                                               gaaGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAattt                       >  1:1157366/1‑76 (MQ=255)
                                               gaaGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAattt                       <  2:690601/76‑1 (MQ=255)
                                                 aGGTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAatttat                     <  2:543693/76‑1 (MQ=255)
                                                   gTGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAatttattt                   <  1:447140/76‑1 (MQ=255)
                                                    tGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATTTTAAtttattta                  <  1:241427/76‑1 (MQ=255)
                                                    tGAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAAtttattta                  <  2:35462/76‑1 (MQ=255)
                                                     gAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTaa                 <  1:861022/76‑1 (MQ=255)
                                                     gAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTaa                 >  2:485634/1‑76 (MQ=255)
                                                     gAAGATATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTaa                 >  2:485629/1‑76 (MQ=255)
                                                     gAAGAGATTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTaa                 <  1:861014/76‑1 (MQ=255)
                                                          tatTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAatgatg            <  1:795094/76‑1 (MQ=255)
                                                          tatTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAatgatg            <  1:500549/76‑1 (MQ=255)
                                                           atTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAAtgatga           >  1:690944/1‑76 (MQ=255)
                                                           atTCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAAtgatga           >  1:1119596/1‑76 (MQ=255)
                                                            ttCTATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAg          >  1:709417/1‑76 (MQ=255)
                                                               tATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGttt       >  1:497680/1‑76 (MQ=255)
                                                               tATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGttt       >  1:1002114/1‑76 (MQ=255)
                                                               tATCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTAAATTTATTTAATGATGAGttt       >  1:497674/1‑76 (MQ=255)
                                                                aTCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGtttt      <  1:1026286/76‑1 (MQ=255)
                                                                 tCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGttttt     <  2:283959/76‑1 (MQ=255)
                                                                 tCCAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGttttt     <  2:283953/76‑1 (MQ=255)
                                                                  ccAGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATATTAATTTATTTAATGATGAGtttttt    >  1:927998/1‑76 (MQ=255)
                                                                    aGGGAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGTTTTTTAc  >  2:996197/1‑76 (MQ=255)
                                                                       |                                                                        
TAGTGCGGCTGTTATTTTCCATTTAATTTACCTTTCCTGATGTGAATGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGTTTTTTAC  >  REL606/702896‑703039

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