Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 745,091 T→A 15.8% C283* (TGT→TGA ybgK → putative allophanate hydrolase, subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913745,0910TA15.8% 186.3 / 20.5 74C283* (TGT→TGAybgKputative allophanate hydrolase, subunit 2
Reads supporting (aligned to +/- strand):  ref base T (40/22);  new base A (7/5);  total (47/27)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATC  >  NC_000913/744992‑745189
                                                                                                   |                                                                                                  
gACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAgtgt                                                                                                    >  1:739061/1‑100 (MQ=255)
gACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAgtgt                                                                                                    <  1:747432/100‑1 (MQ=255)
gACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAgtgt                                                                                                    >  1:957729/1‑100 (MQ=255)
  ccaccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTc                                                                                                  >  1:2573581/1‑100 (MQ=255)
    accGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCAc                                                                                                >  1:1399961/1‑100 (MQ=255)
     ccGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACt                                                                                               >  2:1831/1‑100 (MQ=255)
      cGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTg                                                                                              <  1:2238835/100‑1 (MQ=255)
      cGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTg                                                                                              <  2:2263026/100‑1 (MQ=255)
       ggtggtTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTgg                                                                                             <  2:2423878/100‑1 (MQ=255)
       ggtggtTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTgg                                                                                             <  2:2564184/100‑1 (MQ=255)
        gtggtTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGa                                                                                            <  1:1721678/100‑1 (MQ=255)
           gtTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGAACACTGGacaa                                                                                         <  2:1930893/100‑4 (MQ=255)
            ttACCCGCGTATTGCTTGTATCATTGAGGCTGATATATACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATATTGTCCAGTGAACACTGGacaaa                                                                                        <  1:536226/100‑5 (MQ=255)
              aCCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGc                                                                                      <  2:1176801/100‑1 (MQ=255)
               cccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTTCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCa                                                                                     >  1:2153780/1‑100 (MQ=255)
               cccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGc                                                                                      <  2:2348056/99‑1 (MQ=255)
               cccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCa                                                                                     >  1:647951/1‑100 (MQ=255)
               cccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCa                                                                                     >  2:460781/1‑100 (MQ=255)
               cccGCGGAGTGCCCGTATCATTGAGGCTGATATCTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGAACACTGGacaaaatg                                                                                     <  1:229852/100‑8 (MQ=255)
                ccGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGAACACTGGacaaaatga                                                                                    <  2:1379863/100‑9 (MQ=255)
                 cgcgTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACt                                                                                   <  2:1350686/100‑1 (MQ=255)
                  gcgTATTGCCTTTCTCATTGATGCTGATATGTTCCATCTTGCGTAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGAACTCTGGataaaatga                                                                                    <  1:1539126/98‑9 (MQ=255)
                      aTTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAg                                                                              <  2:1074374/100‑1 (MQ=255)
                       ttGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGc                                                                             <  1:2426617/100‑1 (MQ=255)
                         gCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCAc                                                                           <  2:1604496/100‑1 (MQ=255)
                         gCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCAc                                                                           >  2:972454/1‑100 (MQ=255)
                             gTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCa                                                                       >  1:1139653/1‑100 (MQ=255)
                              tATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCaa                                                                      >  2:709967/1‑100 (MQ=255)
                               aTCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCaa                                                                      <  2:2573600/99‑1 (MQ=255)
                                tCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGa                                                                    <  2:1002644/100‑1 (MQ=255)
                                 cATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGAt                                                                   <  2:2150660/100‑1 (MQ=255)
                                       ggCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCaa                                                             >  1:1040805/1‑100 (MQ=255)
                                         cTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACg                                                           <  2:555696/100‑1 (MQ=255)
                                         cTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACg                                                           <  1:2130959/100‑1 (MQ=255)
                                          tGATCTGTACCACCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGATCACTGGACGCGGCACTAAAAGCACGGCAATATCAGCAAc                                                            >  1:992198/1‑98 (MQ=255)
                                              atGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtt                                                      >  1:73518/1‑100 (MQ=255)
                                                gTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtttt                                                    >  1:791124/1‑100 (MQ=255)
                                                gTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtttt                                                    >  1:340320/1‑100 (MQ=255)
                                                gTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtttt                                                    >  1:1655329/1‑100 (MQ=255)
                                                gTACCATCTGGCGCAAATTCCGCTCGGTCAACCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTAtttt                                                    >  2:1489962/1‑100 (MQ=255)
                                                  aCCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGa                                                  <  2:1682314/100‑1 (MQ=255)
                                                         ngcgcAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTa                                           >  1:427995/2‑100 (MQ=255)
                                                         ggCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTa                                           >  2:2213341/1‑100 (MQ=255)
                                                         ggCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTa                                           >  2:2290675/1‑100 (MQ=255)
                                                          gcgcAAATTCCGCTCGGTAAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCCCTAAAAGCACGGCAAGATCAGCAACGTCATTTTGAACAAtt                                            <  1:820745/98‑1 (MQ=255)
                                                           cgcAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGc                                         >  1:434209/1‑100 (MQ=255)
                                                           cgcAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGc                                         >  2:2211027/1‑100 (MQ=255)
                                                           cgcAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGc                                         >  2:813603/1‑100 (MQ=255)
                                                                aTTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTggc                                    >  2:260797/1‑100 (MQ=255)
                                                                aTTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTggc                                    <  1:799244/100‑1 (MQ=255)
                                                                  tCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTggcgg                                  <  2:1130928/100‑1 (MQ=255)
                                                                   ccGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTggcggc                                 >  2:797601/1‑100 (MQ=255)
                                                                    cGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCt                                >  2:677822/1‑100 (MQ=255)
                                                                    cGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCt                                >  2:2120594/1‑100 (MQ=255)
                                                                    cGCTCCTTCACCCGATTCAATTTGACCAGGGTTCACTAGCAGAGGAACTAAAAGCACGGCAAGATCAACTACGTTATATTGacaaatcacc                                         >  1:2139850/1‑82 (MQ=255)
                                                                     gCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTg                               >  2:152877/1‑100 (MQ=255)
                                                                      cTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGc                              <  2:953159/100‑1 (MQ=255)
                                                                        cGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGcac                            >  2:2378426/1‑100 (MQ=255)
                                                                        cGGTCAGCCGATTCATTTTGTCCAGTGATCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGcac                            >  1:300109/1‑100 (MQ=255)
                                                                         ggTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGcaca                           >  1:472559/1‑100 (MQ=255)
                                                                               ccGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGaaa                     >  1:182930/1‑100 (MQ=255)
                                                                                 gATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAAt                   >  1:930914/1‑100 (MQ=255)
                                                                                 gATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAAt                   >  2:91118/1‑100 (MQ=255)
                                                                                 gATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAAt                   >  2:91214/1‑100 (MQ=255)
                                                                                 gATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAAt                   >  2:2238723/1‑100 (MQ=255)
                                                                                   ttCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTg                 <  1:1151452/100‑1 (MQ=255)
                                                                                      gcctcttCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc              >  2:2225283/7‑100 (MQ=255)
                                                                                      gcctcttCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc              >  2:1008529/7‑100 (MQ=255)
                                                                                      aTTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGAcc              >  2:1297670/1‑100 (MQ=255)
                                                                                        tttGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGTAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTg            >  1:707828/1‑100 (MQ=255)
                                                                                         tcttCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGa           >  2:1583952/4‑100 (MQ=255)
                                                                                          cttCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGaa          >  2:313711/3‑100 (MQ=255)
                                                                                           gTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAAc         >  1:672248/1‑100 (MQ=255)
                                                                                            tCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACg        >  1:272132/1‑100 (MQ=255)
                                                                                            tCCAGTGAACACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACg        >  2:2014156/1‑100 (MQ=255)
                                                                                             ccAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGc       >  1:2249827/1‑100 (MQ=255)
                                                                                             ccAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGc       >  2:155146/1‑100 (MQ=255)
                                                                                                gtgtTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGa    >  1:1844903/1‑100 (MQ=255)
                                                                                                 tgtTCACTGGAAGAGTCACTAACAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGAt   <  1:1807922/100‑1 (MQ=255)
                                                                                                  gtTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATc  >  2:1083112/1‑100 (MQ=255)
                                                                                                  gtTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATc  >  2:64577/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GACCACCGGTGGTTACCCGCGTATTGCCTGTATCATTGAGGCTGATATGTACCATCTGGCGCAAATTCCGCTCGGTCAGCCGATTCATTTTGTCCAGTGTTCACTGGAAGAGGCACTAAAAGCACGGCAAGATCAGCAACGTTATTTTGAACAATTAGCGTGGCGGCTGCACAATGAAAATTGACCTGAACGCCGATC  >  NC_000913/744992‑745189

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