Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,876,374 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,3740GT98.2% 158.3 / ‑4.5 56intergenic (‑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/23);  total (33/23)
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.43e-01

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

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: