Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,183,656 T→A intergenic (‑5/+60) yegX ← / ← thiD putative family 25 glycosyl hydrolase/hydroxy‑methylpyrimidine kinase and hydroxy‑phosphomethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,183,6560TA98.4% 215.6 / ‑3.7 62intergenic (‑5/+60)yegX/thiDputative family 25 glycosyl hydrolase/hydroxy‑methylpyrimidine kinase and hydroxy‑phosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (26/35);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

CGTCTGACGCGGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATA  >  NC_000913/2183559‑2183752
                                                                                                 |                                                                                                
cGTCTGACGCGGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATc                                                                                                <  2:669995/100‑1 (MQ=255)
   cTGACGCGGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTatcatc                                                                                             <  1:575403/100‑1 (MQ=255)
    tGACGCGGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATcc                                                                                            >  2:690334/1‑100 (MQ=255)
     gACGCGGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATccc                                                                                           <  1:101987/100‑1 (MQ=255)
     gACGCGGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATccc                                                                                           >  1:1244520/1‑100 (MQ=255)
      aCGCGGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCt                                                                                          >  1:35861/1‑100 (MQ=255)
       cgcgGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTg                                                                                         <  2:934008/100‑1 (MQ=255)
       cgcgGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTGTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTg                                                                                         >  1:385955/1‑100 (MQ=255)
        gcgGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGa                                                                                        >  2:195428/1‑100 (MQ=255)
           gATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATg                                                                                     <  2:325185/100‑1 (MQ=255)
           gATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATg                                                                                     <  2:749212/100‑1 (MQ=255)
               aaTTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAg                                                                                 <  1:228318/100‑1 (MQ=255)
                   gCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGcc                                                                             >  1:328771/1‑100 (MQ=255)
                    cAACGATGGAAATGAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCa                                                                            >  1:875917/1‑100 (MQ=255)
                       cGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAg                                                                         <  1:839764/100‑1 (MQ=255)
                       cGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAg                                                                         <  1:1227321/100‑1 (MQ=255)
                       cGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAg                                                                         <  2:905677/100‑1 (MQ=255)
                          tGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCt                                                                      <  1:724394/100‑1 (MQ=255)
                           ggAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTg                                                                     <  2:911568/100‑1 (MQ=255)
                           ggAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTg                                                                     <  2:707804/100‑1 (MQ=255)
                            gAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGa                                                                    >  1:295027/1‑100 (MQ=255)
                             aaaTTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGa                                                                    >  1:1261990/1‑99 (MQ=255)
                             aaaTTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGAc                                                                   <  2:411649/100‑1 (MQ=255)
                             aaaTTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGAc                                                                   <  1:164371/100‑1 (MQ=255)
                               aTTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACaa                                                                 >  2:486977/1‑100 (MQ=255)
                                 tAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATa                                                               <  1:1000960/100‑1 (MQ=255)
                                    cccAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAAcc                                                            >  1:1283936/1‑100 (MQ=255)
                                    cccAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAAcc                                                            <  1:852325/100‑1 (MQ=255)
                                    cccAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAAcc                                                            <  2:1264304/100‑1 (MQ=255)
                                     ccAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCa                                                           >  2:569558/1‑100 (MQ=255)
                                     ccAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCa                                                           >  1:876099/1‑100 (MQ=255)
                                      cAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCaa                                                          <  2:348960/100‑1 (MQ=255)
                                      cAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCaa                                                          <  1:1277433/100‑1 (MQ=255)
                                      cAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCaa                                                          <  1:707386/100‑1 (MQ=255)
                                          agCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAagca                                                      >  1:378165/1‑100 (MQ=255)
                                          agCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAagca                                                      >  2:557166/1‑100 (MQ=255)
                                              cATAATAGAGAGGTGAGTTTTTTTTTTTTTTTTATTCTTAGATGCATTGGTATCATTCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCta                                                  >  2:523960/1‑100 (MQ=255)
                                              cATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCta                                                  >  1:733236/1‑100 (MQ=255)
                                               ataataGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCtat                                                 <  2:353475/100‑1 (MQ=255)
                                                taataGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCtata                                                >  1:871635/1‑100 (MQ=255)
                                                    agagagGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTAc                                            <  1:1180713/100‑1 (MQ=255)
                                                       gagGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACgg                                          >  2:382236/1‑99 (MQ=255)
                                                         ggTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGc                                       <  2:1230367/100‑1 (MQ=255)
                                                           tGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTc                                     <  1:557170/100‑1 (MQ=255)
                                                                 tttttCGACTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCaccacc                               <  2:1185958/100‑1 (MQ=255)
                                                                         cTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTgg                       <  2:875911/100‑1 (MQ=255)
                                                                         cTGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTgg                       <  2:884739/100‑1 (MQ=255)
                                                                          tGGTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGa                      <  1:705955/100‑1 (MQ=255)
                                                                           ggTAATTCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGaa                     >  1:862324/1‑100 (MQ=255)
                                                                                ttCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAgtggt                <  1:195430/100‑1 (MQ=255)
                                                                                ttCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAgtggt                <  1:1151989/100‑1 (MQ=255)
                                                                                 tCTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAgtggtg               <  1:1160375/100‑1 (MQ=255)
                                                                                  cTTAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGa              >  2:270674/1‑100 (MQ=255)
                                                                                    tAGCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAAc            >  2:522570/1‑100 (MQ=255)
                                                                                      gCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCg          >  2:1293450/1‑100 (MQ=255)
                                                                                      gCTGCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCg          >  1:298574/1‑100 (MQ=255)
                                                                                         gCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGAc       >  2:465117/1‑100 (MQ=255)
                                                                                         gCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGAc       <  2:1283925/100‑1 (MQ=255)
                                                                                         gCATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGAc       <  1:210437/100‑1 (MQ=255)
                                                                                          cATTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGAcc      <  1:184524/100‑1 (MQ=255)
                                                                                          cATTCGTATCATCCCTGAATGTCACTGCCAAAGGCTGACAATAACCAAAGCAGCTATCGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGAcc      <  2:777960/100‑1 (MQ=255)
                                                                                           aTTGGTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCa     <  1:721241/100‑1 (MQ=255)
                                                                                              ggTATCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATa  >  2:429418/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CGTCTGACGCGGATAAATTGCAACGATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATA  >  NC_000913/2183559‑2183752

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