Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,051 A→G 9.3% Q103R (CAG→CGG) ‡ ygiB → DUF1190 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,0510AG 9.3% 159.0 / 5.1 53Q103R (CAG→CGG) ‡ygiBDUF1190 family protein
Reads supporting (aligned to +/- strand):  ref base A (19/29);  new base G (2/3);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCA  >  NC_000913/3179903‑3180199
                                                                                                                                                    |                                                                                                                                                    
tCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCa                                                                                                                                                      <  1:316514/149‑1 (MQ=255)
  aGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAgg                                                                                                                                                    >  4:49835/1‑149 (MQ=255)
      gCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtgg                                                                                                                                                >  1:136758/1‑149 (MQ=255)
        aaaCCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGATCGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCTGGGCtggtg                                                                                                                                              <  1:426892/149‑3 (MQ=255)
           cccAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCTGGGCtggtgcct                                                                                                                                           <  2:160516/149‑6 (MQ=255)
              aGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCAc                                                                                                                                        <  4:127028/149‑1 (MQ=255)
               ggCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCAcc                                                                                                                                       <  3:170425/149‑1 (MQ=255)
                gCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCTGGGCtggt                                                                                                                                               >  2:513389/1‑139 (MQ=255)
                gCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCTGGGCtggt                                                                                                                                               <  1:513389/140‑2 (MQ=255)
                    aaGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGaaa                                                                                                                                  >  2:530363/1‑149 (MQ=255)
                     aGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGaaaa                                                                                                                                 <  3:21219/149‑1 (MQ=255)
                       cgcCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAAcc                                                                                                                               <  2:271586/149‑1 (MQ=255)
                        gcCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCa                                                                                                                              <  1:278702/149‑1 (MQ=255)
                         ccGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAg                                                                                                                             <  3:400047/149‑1 (MQ=255)
                            aaTGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGgcg                                                                                                                          >  4:3456/1‑149 (MQ=255)
                               gTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAg                                                                                                                       <  3:249894/149‑1 (MQ=255)
                                   caccGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGccc                                                                                                                   >  4:288699/1‑149 (MQ=255)
                                   caccGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGccc                                                                                                                   >  2:74114/1‑149 (MQ=255)
                                     ccGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAg                                                                                                                 <  1:74114/149‑1 (MQ=255)
                                     ccGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAg                                                                                                                 <  3:553022/149‑1 (MQ=255)
                                           acaacaATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATcc                                                                                                           >  1:598320/1‑149 (MQ=255)
                                             aacaaTGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAg                                                                                                         >  2:254201/1‑149 (MQ=255)
                                               caaTGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCg                                                                                                       >  2:167670/1‑149 (MQ=255)
                                                 aTGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCggg                                                                                                     >  2:14139/1‑149 (MQ=255)
                                                   gcgcTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAg                                                                                                   >  3:559797/1‑149 (MQ=255)
                                                       tGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTc                                                                                               <  3:49835/149‑1 (MQ=255)
                                                         aaaGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTg                                                                                             >  1:191844/1‑149 (MQ=255)
                                                         aaaGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTg                                                                                             <  3:3456/149‑1 (MQ=255)
                                                                 cGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGc                                                                                     <  2:136758/149‑1 (MQ=255)
                                                                  gAACGTACTGCGCCGAAATACGCCCCCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCt                                                                                    <  1:240587/149‑1 (MQ=255)
                                                                        aCTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGc                                                                              <  2:65100/149‑1 (MQ=255)
                                                                          tGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCg                                                                            >  3:576087/1‑149 (MQ=255)
                                                                               cGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTAc                                                                       >  2:309921/1‑149 (MQ=255)
                                                                                     aCGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACatgatg                                                                 <  3:582739/149‑1 (MQ=255)
                                                                                            ccGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTc                                                          <  3:356178/149‑1 (MQ=255)
                                                                                                gAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGAt                                                      <  1:43586/149‑1 (MQ=255)
                                                                                                    aCTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGggc                                                  <  2:191844/149‑1 (MQ=255)
                                                                                                                aaTTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGa                                      <  1:287122/149‑1 (MQ=255)
                                                                                                                                gTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGc                      <  3:134438/149‑1 (MQ=255)
                                                                                                                                gTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGc                      <  1:348311/149‑1 (MQ=255)
                                                                                                                                       gccatgccAGCCTGGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTtctc               >  3:162345/7‑149 (MQ=255)
                                                                                                                                       cagGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTtctc               >  3:273926/1‑149 (MQ=255)
                                                                                                                                        agGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTctcc              >  2:511368/1‑149 (MQ=255)
                                                                                                                                        agGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTctcc              <  1:254201/149‑1 (MQ=255)
                                                                                                                                        acgCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTctcc              <  3:171312/147‑1 (MQ=255)
                                                                                                                                         ggCACCAGCCCAGGCTGGCATGGCACCAGACAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTctcct             <  4:299604/149‑1 (MQ=255)
                                                                                                                                         ggCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTctcct             <  1:14139/149‑1 (MQ=255)
                                                                                                                                             ccAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGa                                                                                           <  4:113873/67‑1 (MQ=255)
                                                                                                                                             ccAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGa                                                                                           >  3:113873/1‑67 (MQ=255)
                                                                                                                                              cAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGaaa        >  3:563700/1‑149 (MQ=255)
                                                                                                                                                gCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGaaaaa      <  2:598320/149‑1 (MQ=255)
                                                                                                                                                 cccAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGAtt                                    <  2:133315/118‑1 (MQ=255)
                                                                                                                                                 cccAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGAtt                                    >  1:133315/1‑118 (MQ=255)
                                                                                                                                                   cAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAAccc   <  4:273926/149‑1 (MQ=255)
                                                                                                                                                    aGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACcca  <  1:503337/149‑1 (MQ=255)
                                                                                                                                                    aGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACcca  <  4:372315/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCCA  >  NC_000913/3179903‑3180199

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