Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,231,778 G→T intergenic (+64/‑66) yohK → / → cdd LrgB family inner membrane protein/cytidine/deoxycytidine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,231,7780GT96.4% 196.7 / 0.7 55intergenic (+64/‑66)yohK/cddLrgB family inner membrane protein/cytidine/deoxycytidine deaminase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (29/24);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

CGCTAGCGCTGGTGTTATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGC  >  NC_000913/2231632‑2231917
                                                                                                                                                  |                                                                                                                                           
cGCTAGCGCTGGTGTTATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATca                                                                                                                                         <  1:8037/151‑1 (MQ=255)
  cTAGCGCTGGTGTTATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATcac                                                                                                                                        >  1:829623/1‑150 (MQ=255)
    aGCGCTGGTGTTATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTACTGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACat                                                                                                                                      >  2:911505/1‑150 (MQ=255)
         tGGTGTTATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATaa                                                                                                                                  >  2:527524/1‑149 (MQ=255)
          ggTGTTATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGc                                                                                                                               <  1:318148/151‑1 (MQ=255)
             gtTATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCca                                                                                                                             <  1:204618/150‑1 (MQ=255)
              ttATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCcac                                                                                                                            >  1:268389/1‑150 (MQ=255)
               tATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTGCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCcaca                                                                                                                           <  2:966807/150‑1 (MQ=255)
                aTGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAAc                                                                                                                         <  2:27021/151‑1 (MQ=255)
                     ggATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGtt                                                                                                                    <  1:975204/151‑1 (MQ=255)
                      gATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTc                                                                                                                   <  1:571754/151‑1 (MQ=255)
                      gATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTc                                                                                                                   <  2:1129096/151‑1 (MQ=255)
                        tAATTACTCAGCTGATCGCACCGTTCCTTTTCCAGATTATTCTCGCAGTAATGGGCTAAAATCTGCGACGCGTTGCGCATTTTTGATGTATGTTTCACGCGTAGCATAATTAATGAGATTCATATGACATATATGGTCACAACGGGTTCg                                                                                                                  >  1:1072530/1‑150 (MQ=255)
                         aaTTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGt                                                                                                                 >  1:364965/1‑150 (MQ=255)
                             aCTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACt                                                                                                            >  2:160549/1‑151 (MQ=255)
                             aCTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACt                                                                                                            >  2:572601/1‑151 (MQ=255)
                                        cGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCAtt                                                                                                 >  2:177285/1‑151 (MQ=255)
                                                  cTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTAt                                                                                       >  1:187731/1‑151 (MQ=255)
                                                   ttttACCGATTATTCTTGAAGTAATGGGTTAAAATTTGCGATGCGTCGCGCATTTTTGATGTCTGTTTGACGGGTTGCATAATTTATGAGATTCATATCACATATTAAGCCACACCGGGTTCGTAAATTGTTATCCCATTACATGAtt                                                                                         >  1:971487/1‑148 (MQ=255)
                                                       cccGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGc                                                                                  <  2:969933/151‑1 (MQ=255)
                                                       cccGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGc                                                                                  >  1:935461/1‑151 (MQ=255)
                                                         cGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCa                                                                                 >  1:1205240/1‑150 (MQ=255)
                                                          gATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCaa                                                                                >  1:297337/1‑150 (MQ=255)
                                                          gATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAAc                                                                               >  2:1159309/1‑151 (MQ=255)
                                                                   ggCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCa                                                                       <  1:795741/150‑1 (MQ=255)
                                                                   ggCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCAt                                                                      <  1:289502/151‑1 (MQ=255)
                                                                      aGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCa                                                                   <  2:182192/151‑1 (MQ=255)
                                                                       gTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCAc                                                                  >  2:595356/1‑151 (MQ=255)
                                                                       gTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCAc                                                                  <  2:784133/151‑1 (MQ=255)
                                                                                tAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCaa                                                          >  1:392116/1‑150 (MQ=255)
                                                                                 aaaaTTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAAc                                                        >  1:283729/1‑151 (MQ=255)
                                                                                  aaaTTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAATCCACACCGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAAcc                                                       >  2:1078239/1‑151 (MQ=255)
                                                                                   aaTTTGCGATGCGTCGGGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCg                                                      <  2:634333/151‑1 (MQ=255)
                                                                                    aTTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGc                                                     >  2:802562/1‑151 (MQ=255)
                                                                                    aTTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGc                                                     <  2:268394/151‑1 (MQ=255)
                                                                                    aTTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGc                                                     >  2:123027/1‑151 (MQ=255)
                                                                                      ttGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtt                                                   >  1:389175/1‑151 (MQ=255)
                                                                                      ttGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCt                                                    <  2:940590/150‑1 (MQ=255)
                                                                                       tGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtt                                                   <  1:853368/150‑1 (MQ=255)
                                                                                        gCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtttt                                                 >  1:848945/1‑151 (MQ=255)
                                                                                             gCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCa                                            <  1:1242031/151‑1 (MQ=255)
                                                                                                  gcgcATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTg                                       <  2:391112/151‑1 (MQ=255)
                                                                                                        ttttGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGAt                                  >  2:863418/1‑150 (MQ=255)
                                                                                                           tGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACt                              >  2:1258046/1‑151 (MQ=255)
                                                                                                               gtatgtTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGc                           <  1:326218/150‑1 (MQ=255)
                                                                                                                tatgtTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCa                          >  2:91992/1‑150 (MQ=255)
                                                                                                                   gtTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATc                       <  2:870024/150‑1 (MQ=255)
                                                                                                                   gtTTCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCt                      <  2:1150910/151‑1 (MQ=255)
                                                                                                                    tttCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTg                     >  1:593083/1‑151 (MQ=255)
                                                                                                                     ttCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTg                     <  2:1109954/150‑1 (MQ=255)
                                                                                                                     ttCACGCGTTGCATAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGc                    >  2:370730/1‑151 (MQ=255)
                                                                                                                                 aTAATTAATGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTAt        >  2:399969/1‑151 (MQ=255)
                                                                                                                                   aattaatGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTc      >  2:79251/1‑151 (MQ=255)
                                                                                                                                      taatGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTg    <  1:138715/150‑1 (MQ=255)
                                                                                                                                       aatGAGATTCATATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGc  <  2:610409/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
CGCTAGCGCTGGTGTTATGCGGGATAATTACTTCGCTGATCGCACCGTTCCTTTTCCCGATTATTCTGGCAGTAATGGGCTAAAATTTGCGATGCGTCGCGCATTTTTGATGTATGTTTCACGCGTTGCATAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGC  >  NC_000913/2231632‑2231917

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