Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,876,376 G→T intergenic (‑47/‑205) cysD ← / → iap sulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,876,3760GT98.2% 156.8 / ‑4.5 55intergenic (‑47/‑205)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (31/22);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

CGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGT  >  NC_000913/2876240‑2876515
                                                                                                                                         |                                                                                                                                           
cGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACtata                                                                                                                                       <  2:858847/143‑1 (MQ=255)
cGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACtata                                                                                                                                       <  1:831958/143‑1 (MQ=255)
   aTTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTGGTATTTGATCCATAACCGTTCCTTTGCAATCCCGCTATTTTCTTGCACTCAGATCTTTCTTCTATAgtg                                                                                                                                    >  1:871351/1‑141 (MQ=255)
     ttGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGCTTCTACTATAGGGAg                                                                                                                                  >  2:866089/1‑143 (MQ=255)
        agaATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGt                                                                                                                               >  1:412525/1‑143 (MQ=255)
        agaATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGt                                                                                                                               >  1:218938/1‑143 (MQ=255)
         gaATTCTGCCGCCACCTCGCGAATGATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTa                                                                                                                              <  1:829078/143‑1 (MQ=255)
          aaTTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCCGCTGTTTCTACTATAAGGAGCGTaa                                                                                                                             >  1:845780/1‑143 (MQ=255)
          aaTTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGATATTTTCTTGCCTTGAAATGTTTCTACTATAGGGAGCGTaa                                                                                                                             <  2:834322/143‑1 (MQ=255)
            ttCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAaga                                                                                                                           <  2:371608/143‑1 (MQ=255)
                gccgccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCTGCTAGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAAc                                                                                                                       <  1:649738/143‑1 (MQ=255)
                gccgccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAAc                                                                                                                       <  1:519216/143‑1 (MQ=255)
                 ccgccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACg                                                                                                                      >  1:757573/1‑143 (MQ=255)
                 ccgccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTCTGCCATACCGCTATTTTCTTGCCATCAGATGTTTCTACTACAGGGAGCGTAAGAGAACg                                                                                                                      >  2:468902/1‑143 (MQ=255)
                  cgccACCTCGCGAAGAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGCACGa                                                                                                                     >  2:429931/1‑143 (MQ=255)
                         tCGCGAATAATGTGGATGCTTTCCGCCTGCAGTTGCCGCAGGTGAGTAAGTGGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTATGTTGCCATCACAAGCTTCAAATATAAGGACCGCACCAGAACaaacacacg                                                                                                              >  2:450244/1‑137 (MQ=255)
                           gcgAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTGGATCCATAACCGTTCCTTTGCAATACCGCTAGTTTCTTGCCATCAGCTGTTTCTACTCTAGTGCGCGTCACAGACCGAAAGaccttg                                                                                                            >  1:288391/1‑138 (MQ=255)
                                     tgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAGGGGAACGAATGAAATTACCATTTAGaa                                                                                                  >  1:701804/1‑143 (MQ=255)
                                           cTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTCACAATTAGACTGagta                                                                                            >  2:303072/1‑143 (MQ=255)
                                                   cTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTa                                                                                    <  1:871962/143‑1 (MQ=255)
                                                    tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTaa                                                                                   <  1:110323/143‑1 (MQ=255)
                                                     ccAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAAc                                                                                  >  1:772479/1‑143 (MQ=255)
                                                      cAGTGGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTACTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACg                                                                                 <  1:905806/143‑1 (MQ=255)
                                                        gTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGa                                                                               <  1:813623/143‑1 (MQ=255)
                                                              gCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCCATACAACCAAAACGCATGGAGTAGTTGTAAAGACGAAGCGGAGCGTAAGAGAACGAATGAAATTGCCGATTAGAATGAGTAGTTCCTTAACGGAATAACg                                                                         >  1:442832/1‑143 (MQ=255)
                                                                 ggTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGAtt                                                                      >  2:36762/1‑143 (MQ=255)
                                                                  gtgtggaaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATGACGAttt                                                                     >  2:792300/7‑143 (MQ=255)
                                                                    gAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTgg                                                                   <  1:622035/143‑1 (MQ=255)
                                                                             gTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCCACTAGCATGCGTACTTCCTTAACGGAATAACGATTTGGCAAAGCTaa                                                          >  1:483965/1‑143 (MQ=255)
                                                                                 ttGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATc                                                      >  2:487122/1‑143 (MQ=255)
                                                                                 ttGATCCATAACCGTTCCTTTGCAATACCGCTATTTCCTTGCCATCAGATGTTTCTACAGTAGGGAGCGGACGCGAGCGAATGAAATCACCACTTAGAATGCGGAGCTCCTTACCGGGAGAACGAGTTGCGAAAGCTAATATc                                                      >  2:490483/1‑143 (MQ=255)
                                                                                   gATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTCAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTTGCAAAGCTAATATCaa                                                    >  2:348628/1‑143 (MQ=255)
                                                                                    aTCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaa                                                   >  2:859931/1‑143 (MQ=255)
                                                                                     tCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGACAAAGCTAATATCaaaa                                                  >  2:642062/1‑143 (MQ=255)
                                                                                           aCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGCTTTGGCAAAGCTAATAGCAAAAAGTGCt                                            >  1:112138/1‑143 (MQ=255)
                                                                                           aCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGCAGCGTCAGAGAACGAATGAGATTACCAACTAGAATGAGTAGTTCCTTAACGGCATCACGCATTGGCAACCCTAACGTCaaattgtggt                                            >  1:431949/1‑136 (MQ=255)
                                                                                              gTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCACTCCACGCCCACATTTGTTTAGGGGAATACCGATTTGGCAAAGCTAATATCAAAACGCGCTGaa                                         >  1:945809/1‑143 (MQ=255)
                                                                                               ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAGTGAAATTACCAATTAGCATGAGTAGGTCCTTAGCGGCATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTAAAg                                        >  1:72057/1‑143 (MQ=255)
                                                                                                  cTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATACCGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCa                                     >  1:300309/1‑143 (MQ=255)
                                                                                                  cTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCa                                     >  2:596414/1‑143 (MQ=255)
                                                                                                   tttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAc                                    >  1:373900/1‑143 (MQ=255)
                                                                                                   tttGCAATACCGCTATTTTCTTACCATCAGATGTTTCTACTAGAGAGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTGACTTAACGGAATAACGATGTGGGACAGCTAATATCAAGAAGTGCTTAAGGCAc                                    >  2:248627/1‑143 (MQ=255)
                                                                                                      gCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCgg                                 <  1:362232/143‑1 (MQ=255)
                                                                                                            cctgtATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCTATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                           <  2:63228/139‑1 (MQ=255)
                                                                                                            ccGCTATTTTGTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                           <  1:424110/143‑1 (MQ=255)
                                                                                                             cGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCgg                          <  1:21400/143‑1 (MQ=255)
                                                                                                              gCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCggg                         <  1:331704/143‑1 (MQ=255)
                                                                                                                   tttCTTGCCATCAGATGGTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGttt                    <  1:305588/143‑1 (MQ=255)
                                                                                                                     gtCTTGCCATCAGATGGTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTCACGCACCGGATTTCGGGCGTTTa                   <  2:835897/142‑1 (MQ=255)
                                                                                                                           ccATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGa             >  2:840437/1‑143 (MQ=255)
                                                                                                                             aTCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGCGTAGTTCCTTAACGGAATAACGATCTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAtc           >  2:474019/1‑142 (MQ=255)
                                                                                                                                ataTGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATGAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCGTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGa        <  1:51885/141‑1 (MQ=255)
                                                                                                                                  agttttCTAGTAGAGGGACCGTAAGAGAACGAATGAAATTACCAACTAGAATTAGTAGTTTGTTAACGGAATAACGCTTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGaaa      <  2:58525/140‑1 (MQ=255)
                                                                                                                                  aTGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGCTTAGGTAGATTTGaaa      >  1:589102/1‑143 (MQ=255)
                                                                                                                                      ttCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGt  <  1:468902/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
CGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGT  >  NC_000913/2876240‑2876515

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