Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,050 C→T 9.8% Q103* (CAG→TAG) ‡ ygiB → DUF1190 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,0500CT 9.8% 148.2 / 5.7 51Q103* (CAG→TAG) ‡ygiBDUF1190 family protein
Reads supporting (aligned to +/- strand):  ref base C (18/28);  new base T (2/3);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCC  >  NC_000913/3179903‑3180198
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                                              <  1:513389/140‑2 (MQ=255)
                gCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCTGGGCtggt                                                                                                                                              >  2:513389/1‑139 (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                                                                                                                  >  2:74114/1‑149 (MQ=255)
                                   caccGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGccc                                                                                                                  >  4:288699/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                                                                                            <  3:3456/149‑1 (MQ=255)
                                                         aaaGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTg                                                                                            >  1:191844/1‑149 (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)
                                                                                                                                                   |                                                                                                                                                    
TCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTCCTCGAAAAACCC  >  NC_000913/3179903‑3180198

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: ‑