Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,716,360 C→A 9.2% W209C (TGG→TGT)  kipI ← possible allophanate hydrolase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,716,3600CA 9.2% 161.6 / 5.4 51W209C (TGG→TGT) kipIpossible allophanate hydrolase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (18/27);  new base A (3/2);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACCGTTGAAAAGAGACCTGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGATTATTCGCAATACCTACAG  >  CP000730/1716220‑1716496
                                                                                                                                            |                                                                                                                                        
gTACCGTTGAAAAGAGACCTGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATc                                                                                                                                          >  1:270449/1‑141 (MQ=255)
        gAAAAGAGACCTGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCTATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCTTTAAATTTTATTCATTCTcc                                                                                                                                  >  2:644747/1‑141 (MQ=255)
            agagACCTGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTACAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCt                                                                                                                              <  2:330053/141‑1 (MQ=255)
              agaCCTGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTc                                                                                                                            >  2:655774/1‑141 (MQ=255)
                 ccTGGTTGTAAAATCTTAATTGGCATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATACATTCTCCTGCTTCGTa                                                                                                                         >  2:693534/1‑141 (MQ=255)
                 ccTGGTTGTAAAATCTTAATTGCCATTTTCAATCACCACCCAGTCATCAACATTAAGGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATACATTCTCCTGCTTCGTa                                                                                                                         >  2:694271/1‑141 (MQ=255)
                 ccTGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATACATTCTCCTGCTTCGTa                                                                                                                         <  1:694271/141‑1 (MQ=255)
                 ccTGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATACATTCTCCTGCTTCGTa                                                                                                                         <  1:693534/141‑1 (MQ=255)
                  cTGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTac                                                                                                                        >  1:549817/1‑141 (MQ=255)
                   tGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTaca                                                                                                                       >  2:524744/1‑141 (MQ=255)
                       tGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTacattga                                                                                                                   <  1:263964/141‑1 (MQ=255)
                         tAAAATCTTAATTGACATGTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTacattgaca                                                                                                                 <  1:682910/141‑1 (MQ=255)
                                   aTTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCAcg                                                                                                       <  2:220565/141‑1 (MQ=255)
                                    ttGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCAcgc                                                                                                      >  2:306446/1‑141 (MQ=255)
                                     tGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGGTCAcgct                                                                                                     <  1:549087/141‑1 (MQ=255)
                                           tttCAATCACTACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCgctgc                                                                                               <  1:10692/141‑1 (MQ=255)
                                           tttCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCgctgc                                                                                               <  1:707640/141‑1 (MQ=255)
                                           tttCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCgctgc                                                                                               >  2:480272/1‑141 (MQ=255)
                                                aTCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCAATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAAt                                                                                          <  2:233919/141‑1 (MQ=255)
                                                  caccacCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATAc                                                                                        >  1:532607/1‑141 (MQ=255)
                                                   accaccCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACt                                                                                       <  2:399460/141‑1 (MQ=255)
                                                   accaccCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACt                                                                                       >  1:276317/1‑141 (MQ=255)
                                                          aatcatcaACATTAAAGTGGCCAGCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAACTTCTATCCATGCTCCTTCTGCCTACATTGACATTAGTTCACGCTCGCTGCTAAATACTTGTAACg                                                                                <  2:496403/139‑1 (MQ=255)
                                                            tcatcaACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGgt                                                                              >  1:374387/1‑141 (MQ=255)
                                                                caACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGc                                                                          >  1:156176/1‑141 (MQ=255)
                                                                             gCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTg                                                             >  1:451045/1‑141 (MQ=255)
                                                                               cATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCGCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGcc                                                           <  2:280606/141‑1 (MQ=255)
                                                                                 tCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCAt                                                         <  2:571919/141‑1 (MQ=255)
                                                                                   tGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATcc                                                       <  1:405888/141‑1 (MQ=255)
                                                                                     atatatCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGc                                                     <  1:440726/141‑1 (MQ=255)
                                                                                        tatCTCTTTCGATTTGTATAAATTTCTGTTTATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAg                                                  <  1:671062/141‑1 (MQ=255)
                                                                                            tctTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGgaga                                              >  2:335277/1‑141 (MQ=255)
                                                                                               gcaagagttattataTTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATc                                           <  2:255025/127‑1 (MQ=255)
                                                                                                  gATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGa                                        <  1:549876/141‑1 (MQ=255)
                                                                                                    tttGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATg                                      >  2:488742/1‑141 (MQ=255)
                                                                                                     ttGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATgg                                     >  1:144582/1‑141 (MQ=255)
                                                                                                          taAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGAtata                                <  2:56516/141‑1 (MQ=255)
                                                                                                           aaaTTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAg                               <  1:528154/141‑1 (MQ=255)
                                                                                                               ttCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCt                           <  1:644747/141‑1 (MQ=255)
                                                                                                                 cTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGt                         <  2:668304/141‑1 (MQ=255)
                                                                                                                  tGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGtt                        <  1:433902/141‑1 (MQ=255)
                                                                                                                    ttCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGACTCTGATGGATATAGTCCTGTTTg                      <  2:200112/141‑1 (MQ=255)
                                                                                                                      cATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGat                    <  1:524744/141‑1 (MQ=255)
                                                                                                                       aTCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGatt                   <  1:426764/141‑1 (MQ=255)
                                                                                                                             tGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCCTACGCCAGGAGAATCTGATGGATATAGTCCTGTTTGATTATTCGc             >  1:540092/1‑141 (MQ=255)
                                                                                                                              gCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGATTATTCGCa            <  1:605427/141‑1 (MQ=255)
                                                                                                                              gCATAAAATAGTATACATTCTCCTGCTTCGTACATTGACATAGGTTCACGCTCGCTGCTAAATACTTTTAACAGTGTGCGTCCAAAAATTTGCCATCAGCCAGGAAAATCAGATGAAAATAGCCCTGTTTGattatctgca            >  2:610676/1‑136 (MQ=255)
                                                                                                                               cATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGATTATTCGCaa           >  2:673679/1‑141 (MQ=255)
                                                                                                                                 tAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGATTATTCGCAATa         >  2:578619/1‑141 (MQ=255)
                                                                                                                                        gTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGATTATTCGCAATACCTACAg  >  2:678370/1‑141 (MQ=255)
                                                                                                                                        gTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGATTATTCGCAATACCTACAg  >  1:317071/1‑141 (MQ=255)
                                                                                                                                        gTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGATTATTCGCAATACCTACAg  <  2:577731/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
GTACCGTTGAAAAGAGACCTGGTTGTAAAATCTTAATTGACATTTTCAATCACCACCCAGTCATCAACATTAAAGTTGCCATCTGATATATCTCTTTCGATTTGTATAAATTTCTGTTCATCTATTGCATAAAATTGTATCCATTCTCCTGCTTCGTACATTGACATTGGTTCACGCTCGCTGCTAAATACTTTTAACGGTGTGCGTCCAATAATTTGCCATCCGCCAGGAGAATCTGATGGATATAGTCCTGTTTGATTATTCGCAATACCTACAG  >  CP000730/1716220‑1716496

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