Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 669,362 A→G 7.3% D182D (GAT→GAC)  cobC ← putative alpha‑ribazole‑5'‑P phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913669,3620AG 7.3% 168.1 / 3.7 55D182D (GAT→GAC) cobCputative alpha‑ribazole‑5'‑P phosphatase
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base G (2/2);  total (26/29)
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.55e-01

AAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTA  >  NC_000913/669216‑669503
                                                                                                                                                  |                                                                                                                                             
aaaaGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTCCCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGatat                                                                                                                                              >  8:110395/1‑148 (MQ=255)
  aaGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCa                                                                                                                                            <  5:264290/148‑1 (MQ=255)
  aaGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCa                                                                                                                                            <  7:27177/148‑1 (MQ=255)
  aaGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCa                                                                                                                                            <  6:58067/148‑1 (MQ=255)
             cGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCa                                                                                                                                 <  5:43115/148‑1 (MQ=255)
                  aaTATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACAcc                                                                                                                            <  6:90020/148‑1 (MQ=255)
                  aaTATCAGCCTGCCCGGGTCTGGCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACAcc                                                                                                                            <  3:161428/148‑1 (MQ=255)
                                   gTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGaaaa                                                                                                           >  3:105912/1‑148 (MQ=255)
                                            tGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCcaca                                                                                                    >  3:1451/1‑146 (MQ=255)
                                            tGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCcaca                                                                                                    >  3:74373/1‑146 (MQ=255)
                                                cGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCtt                                                                                              >  3:288264/1‑148 (MQ=255)
                                                    gCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCgg                                                                                           >  3:183191/1‑147 (MQ=255)
                                                     cGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCggc                                                                                          <  7:224029/147‑1 (MQ=255)
                                                        cTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCggcaggc                                                                                      <  6:225825/148‑1 (MQ=255)
                                                              cGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCggcagg                                                                                       <  5:316333/141‑1 (MQ=255)
                                                              cGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCggcagg                                                                                       >  6:316333/1‑141 (MQ=255)
                                                              cGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGcc                                                                                 >  6:227564/1‑147 (MQ=255)
                                                              cGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGcc                                                                                 >  7:247377/1‑147 (MQ=255)
                                                                tAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAAt                                                                              >  3:66157/1‑148 (MQ=255)
                                                                tAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAAt                                                                              >  4:272037/1‑148 (MQ=255)
                                                                  aaaaaaCAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCGACATAGCTTCGGCAGGCATGCCAAtt                                                                             <  4:66157/147‑1 (MQ=255)
                                                                  aaaaaaCAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCggcaggc                                                                                      >  8:38142/1‑138 (MQ=255)
                                                                  aaaaaaCAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCggcaggc                                                                                      <  7:38142/138‑1 (MQ=255)
                                                                           aaaaGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGa                                                                   <  4:288264/148‑1 (MQ=255)
                                                                               gTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTaa                                                               <  4:105912/148‑1 (MQ=255)
                                                                               gTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTaa                                                               <  7:110395/148‑1 (MQ=255)
                                                                                tCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCaa                                                                                                                   <  1:135119/95‑1 (MQ=255)
                                                                                tCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCaa                                                                                                                   >  2:135119/1‑95 (MQ=255)
                                                                                        tctcGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTcc                                                                                                                                  >  3:171634/1‑72 (MQ=255)
                                                                                        tctcGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTcc                                                                                                                                  <  4:171634/72‑1 (MQ=255)
                                                                                         ctcGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAg                                                      <  4:183191/147‑1 (MQ=255)
                                                                                         ctcGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAg                                                      <  3:272037/147‑1 (MQ=255)
                                                                                             aCCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTacac                                                 >  5:138441/1‑148 (MQ=255)
                                                                                                 cGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCggcaggc                                                                                      >  1:286607/1‑107 (MQ=255)
                                                                                                 cGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCggcaggc                                                                                      <  2:286607/107‑1 (MQ=255)
                                                                                                  gATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTg                                             <  6:162172/147‑1 (MQ=255)
                                                                                                  gATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTg                                             >  5:162172/1‑147 (MQ=255)
                                                                                                          ggCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGGTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTgac                                    >  4:160540/1‑148 (MQ=255)
                                                                                                          ggCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTgac                                    <  8:247377/148‑1 (MQ=255)
                                                                                                           gCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCa                                                       <  8:291728/128‑1 (MQ=255)
                                                                                                           gCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCa                                                       >  7:291728/1‑128 (MQ=255)
                                                                                                                 gAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTaa                              >  4:96014/1‑147 (MQ=255)
                                                                                                                    gACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAata                          >  8:32941/1‑148 (MQ=255)
                                                                                                                    gACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAata                          >  2:266648/1‑148 (MQ=255)
                                                                                                                      cgcgTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACAGCCTGATGGCTGACGACTAAAATAtt                        <  5:227564/148‑1 (MQ=255)
                                                                                                                      cgcgTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATAtt                        <  7:32941/148‑1 (MQ=255)
                                                                                                                        cgTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGCTCATCTCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGTTTACCAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCt                      <  8:149499/148‑1 (MQ=255)
                                                                                                                         gTAGCGTCGCGAATTTTTGGTTGATGTCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCt                      >  5:106534/1‑147 (MQ=255)
                                                                                                                         gTAGCGTCGCGAATTTTTGGTTGATGTCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCt                      >  5:111328/1‑147 (MQ=255)
                                                                                                                          tAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGa                    >  8:217205/1‑148 (MQ=255)
                                                                                                                           aGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGa                    >  2:94497/1‑147 (MQ=255)
                                                                                                                             cGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCAGa                    >  3:269060/1‑145 (MQ=255)
                                                                                                                               tCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTg               <  1:94497/148‑1 (MQ=255)
                                                                                                                                   gAATTTTTGGTTGATGTCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGa           <  6:111328/148‑1 (MQ=255)
                                                                                                                                   gAATTTTTGGTTGATGTCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGa           <  6:106534/148‑1 (MQ=255)
                                                                                                                                            gTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAATGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTa  <  4:1451/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
AAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGTCGCGAATTTTTGGTTGATATCAATGGCGCTCCAACACCCCTGGTCAACGCGAAAATGCCACATAGCTTCGGCAGGCATGCCAATTAAACGGGCGATTAACAGACTCAGTACACCCTGATGGCTGACGACTAAAATATTCTGATAGTGCTGAAATTCACTA  >  NC_000913/669216‑669503

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑