Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,251,804 C→T 7.7% intergenic (‑21/‑290) hmrA ← / → luxS M20D subfamily peptidase/S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,251,8040CT 7.7% 139.1 / 2.6 55intergenic (‑21/‑290)hmrA/luxSM20D subfamily peptidase/S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

GACGCAAATATTTCTTCATTACCAAGTTCAGGACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAATTTTCTAATAACTAAAAGCCTTTCATAAATT  >  CP000730/2251668‑2251944
                                                                                                                                        |                                                                                                                                            
gaCGCAAATATTTCTTCATTAACCAGTTGAGGACGTTCATGAACTCTATGACTGATTTCGATATAACAATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTa                                                                                                                                          >  2:110555/1‑141 (MQ=255)
gaCGCAAATATTTCTTCATAGCCAAGTTCAGGACGTTCATGAATTCTATGACTGGTTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTAGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTa                                                                                                                                          >  1:364395/1‑141 (MQ=255)
   gCAAATATTTCTTCATTACCAAGTTCAGGGCGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTATCCTacac                                                                                                                                       >  1:537449/1‑141 (MQ=255)
   gCAAATATTTCTTCATTACCAAGTTCAGGACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTATCCTaccc                                                                                                                                       >  1:6981/1‑139 (MQ=255)
    cAAATATTTCTTCATTACCAAGTTCAGGACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACt                                                                                                                                      >  1:460157/1‑141 (MQ=255)
        ttgtgcttcATTACCAAGTTCAGGACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCAtt                                                                                                                                  <  2:417945/136‑1 (MQ=255)
        tctttttttcttcccAAGTTCAGGACGTTCATGAATTCTATTACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCAtt                                                                                                                                  <  2:548298/128‑1 (MQ=255)
              ttcATTACCAAGTTCAGGACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAatc                                                                                                                            >  1:321467/1‑141 (MQ=255)
                    aCCAAGTTCAGGACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCcac                                                                                                                      >  1:350393/1‑141 (MQ=255)
                    aCCAAGTTCAGGACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCcac                                                                                                                      <  2:3904/141‑1 (MQ=255)
                    aCCAAGTTCAGGACGTTAATTAATTCTCTGCCTGATTTCGATCTAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCcac                                                                                                                      <  1:502718/141‑1 (MQ=255)
                         gTTCAGGACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAatat                                                                                                                 <  2:538284/141‑1 (MQ=255)
                         gTTCAGGACGGTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAatat                                                                                                                 <  1:221237/141‑1 (MQ=255)
                              ggACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGt                                                                                                            >  1:124590/1‑141 (MQ=255)
                               gACGTTCATGAATTGTATGACAGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGCGAATTTGTTGTTTTGCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGtt                                                                                                           >  2:384594/1‑141 (MQ=255)
                               gACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCTCTTTTTTTTTACTATCACACTCCTACCCTACACTCATCTTATTCATTGACAATtgtttgtt                                                                                                           >  2:284231/1‑141 (MQ=255)
                               gACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGtt                                                                                                           <  2:4921/141‑1 (MQ=255)
                               gACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGtt                                                                                                           >  2:182839/1‑141 (MQ=255)
                               gACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTACTTTTTACCCTGCACTCATTTTAATCATCCACAATATTTTGtt                                                                                                           >  1:240583/1‑141 (MQ=255)
                                      aTGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCaa                                                                                                    >  1:117370/1‑141 (MQ=255)
                                        taaTTCTATGACTGATTACGAGATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCCTTTCACCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCaaaa                                                                                                  <  2:422251/140‑1 (MQ=255)
                                        gAATTCTATGACTGATATCGGTATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGGTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTAATTTTGATCATCCACAGTATTTGTTTCTTTCgaaa                                                                                                  >  2:282311/1‑141 (MQ=255)
                                         aaTTCTATGACTTATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAAt                                                                                                 <  2:224607/141‑1 (MQ=255)
                                         aaTTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTCTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAAt                                                                                                 <  1:216523/141‑1 (MQ=255)
                                           ttCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGa                                                                                               >  1:273593/1‑141 (MQ=255)
                                               ttGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTACTTTCAAAATGAatta                                                                                           <  1:277436/140‑1 (MQ=255)
                                                     gATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATAGTTTGTTTTTTCACAATGAATTATTACtt                                                                                     >  1:346700/1‑141 (MQ=255)
                                                          cGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCATAGAATTTTGTTAGTCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCt                                                                                >  1:16347/1‑141 (MQ=255)
                                                             tataACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATc                                                                             >  1:269664/1‑141 (MQ=255)
                                                              ataACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTATATCg                                                                            >  1:236850/1‑141 (MQ=255)
                                                                           tttGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTATCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCatg                                                               <  2:6981/141‑1 (MQ=255)
                                                                           tttGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTATCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCatg                                                               <  2:537449/141‑1 (MQ=255)
                                                                           tttGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCatg                                                               >  2:283652/1‑141 (MQ=255)
                                                                            ttGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCatgt                                                              >  1:500402/1‑140 (MQ=255)
                                                                              gTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCatgatt                                                            >  1:156748/1‑140 (MQ=255)
                                                                               tttCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGt                                                           <  2:269664/141‑1 (MQ=255)
                                                                                  cTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCAt                                                        <  1:309392/141‑1 (MQ=255)
                                                                                   tatataATCGAGAATTTGTTGTGTTTCACTCACTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGATTTTATCTCATGATGTCATc                                                       <  2:202352/141‑1 (MQ=255)
                                                                                        aaTCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGtt                                                  >  1:558883/1‑141 (MQ=255)
                                                                                        aaTCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGtt                                                  <  1:143191/141‑1 (MQ=255)
                                                                                              gaATTTGTTGTTTGTCTCTCATTTTTTACTATCTCCCTTTTCCCCTACACTCATTTTAATCATCCAAAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCtt                                            <  2:238195/141‑1 (MQ=255)
                                                                                                     ttgttTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTATTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTAc                                     <  1:214945/141‑1 (MQ=255)
                                                                                                      tgttTTTCACTCCTTTTTTACTATCTCCCTTTTACCCTCCACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACa                                    <  2:452944/141‑1 (MQ=255)
                                                                                                       gttTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACCCTCATTTTAATCATCCACAATCTTTTGTTCTTTCAAAATGAATCATTACTTATTCTTTCGGTTTTCTCTCCTGATGTCATCTTGTTTTTCTTTATTTACac                                   >  1:327652/1‑140 (MQ=255)
                                                                                                        tttttCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCCATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACaaa                                  <  2:117370/141‑1 (MQ=255)
                                                                                                         ttttCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACaaaa                                 <  1:12188/141‑1 (MQ=255)
                                                                                                            tCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAAttt                              <  2:460157/141‑1 (MQ=255)
                                                                                                             cACTCATTTTTTACTATCTCGCTTTTACGCTACACTCATTTTAATCATCCACAACATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAAtttt                             >  1:363719/1‑141 (MQ=255)
                                                                                                              aCTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAATTTTc                            >  2:176789/1‑141 (MQ=255)
                                                                                                                      tttACTATCTCCCCTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAATTTTCTAATAACt                    <  2:101577/141‑1 (MQ=255)
                                                                                                                       ttatcaTCTCCCTTTTACCCTACCCTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTACTCTATCGGTTTTATCTCTTGATGTCATCTAGTTTTTCTTTATTTACAAAATTTTCTAATAACTa                   <  1:303276/136‑1 (MQ=255)
                                                                                                                          ccaTCTCCCTTTCACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTGTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAATTTTCTAATAACTaaaa                <  2:101520/139‑1 (MQ=255)
                                                                                                                              ttcCCTTGTAACCTGCACTAATTTTACAGGTGCACGACATTTTGTTCTTTCAAAATGAATTATTACGTATTCTATCGGTCTTATCTCATGATGTCATCTAGTTTTTCGTTATTGACAAAATTTTCTAATAACTAAAAGCCt            <  2:253131/140‑1 (MQ=255)
                                                                                                                               tcCTTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAATTTTCTAATAACTAAAAGCCtt           <  2:46851/141‑1 (MQ=255)
                                                                                                                                     ttACCCTACACTGATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTAATTATTCTATCGGTTTTATCTCATGATGTGCTCTAGTTTTTCTTTATAGACAAAATTTTCTAATAACTAAAAGCGTTTCATaa     >  2:212703/1‑141 (MQ=255)
                                                                                                                                       tcccTCCTCTTATTTTAAACATGCGCAAGATTGTGTCGTTTCAACATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAATTTTCTAATAACTAAAAGCCTTTCATAAAt   <  1:115874/140‑1 (MQ=255)
                                                                                                                                        cccTACACTCATTTTAATCATCCACAATATTCTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAATTTTCTAATAACTAAAAGCCTTTCATAAAtt  >  1:318974/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                            
GACGCAAATATTTCTTCATTACCAAGTTCAGGACGTTCATGAATTCTATGACTGATTTCGATATAACTATATTTATTTGTTTCTATATAATCGAGAATTTGTTGTTTTTCACTCATTTTTTACTATCTCCCTTTTACCCTACACTCATTTTAATCATCCACAATATTTTGTTCTTTCAAAATGAATTATTACTTATTCTATCGGTTTTATCTCATGATGTCATCTAGTTTTTCTTTATTTACAAAATTTTCTAATAACTAAAAGCCTTTCATAAATT  >  CP000730/2251668‑2251944

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑