Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 882,722 G→A 100% R5C (CGC→TGC)  deoR ← deoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913882,7220GA98.6% 228.3 / ‑4.3 72R5C (CGC→TGC) deoRdeoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (37/34);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCG  >  NC_000913/882625‑882815
                                                                                                 |                                                                                             
gTCATCTCCGAAACGCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTACCGCTGCGTTTCAATACGTGCAGCAGCTGCCCGATACGCTCTTCGCAAc                                                                                             >  2:1262512/1‑100 (MQ=255)
 tCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACg                                                                                            >  2:199437/1‑100 (MQ=255)
   aTCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgtg                                                                                          >  1:1559535/1‑100 (MQ=255)
   aTCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgtg                                                                                          >  1:1852930/1‑100 (MQ=255)
    tctcCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgtgt                                                                                         >  1:1468303/1‑100 (MQ=255)
      tcCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGttt                                                                                       >  2:466809/1‑100 (MQ=255)
       ccGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                      <  2:2457489/100‑1 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                <  2:1226075/100‑1 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                <  2:2284186/99‑1 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATc                                                                               >  2:1123720/1‑100 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATc                                                                               >  2:2443125/1‑100 (MQ=255)
               ccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATcc                                                                              >  1:2125428/1‑100 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                         <  2:886114/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                         >  1:1896437/1‑100 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                         <  1:1321586/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                         <  2:520848/100‑1 (MQ=255)
                     aGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGa                                                                        <  1:248777/100‑1 (MQ=255)
                     aGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCctc                                                                           <  2:774577/97‑1 (MQ=255)
                       ggcggcggcGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAAt                                                                      >  2:615846/1‑100 (MQ=255)
                        gcggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAAt                                                                      >  1:1754778/1‑99 (MQ=255)
                         cggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAg                                                                    >  1:576192/1‑100 (MQ=255)
                           gcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGtt                                                                  >  2:1020414/1‑100 (MQ=255)
                                 tCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGaa                                                            >  1:2147572/1‑100 (MQ=255)
                                   tttAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGc                                                          >  1:912367/1‑100 (MQ=255)
                                    ttgaGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTTCATAATCCCTCTGAATAGTTCTTGa                                                             >  1:577180/4‑96 (MQ=255)
                                       aGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGc                                                      <  1:1900409/100‑1 (MQ=255)
                                       aGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGc                                                      <  1:2086931/100‑1 (MQ=255)
                                        gATGTCACTTATCGCTGCCTTTTAATTCGTTCCTCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGc                                                      <  2:996214/99‑1 (MQ=255)
                                        gATGTAACTTAGCGCCGCGTTATAATTCTTGCAGCATCTGCCCGATACGCTCTTCGCAACGTGTTTGCATAATCCCTCTAACTAGTTATTGAAGcca                                                        >  2:295738/1‑95 (MQ=255)
                                         aTGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCg                                                    >  1:393781/1‑100 (MQ=255)
                                            tAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTc                                                 <  1:527839/100‑1 (MQ=255)
                                            tAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTc                                                 <  1:1842585/100‑1 (MQ=255)
                                             aaCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCa                                                >  1:493528/1‑100 (MQ=255)
                                             aaCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCa                                                >  1:1105034/1‑100 (MQ=255)
                                             aaCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCa                                                <  1:1491825/100‑1 (MQ=255)
                                              aCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCaa                                               >  2:1115004/1‑100 (MQ=255)
                                               cTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGGAGCGAGCCGCTCaa                                               <  2:2056642/99‑1 (MQ=255)
                                                    cGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                         <  1:1287852/100‑1 (MQ=255)
                                                    cGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                         <  1:1842141/100‑1 (MQ=255)
                                                     gCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTaca                                        >  2:1475974/1‑100 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                       >  2:1111576/1‑100 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                       >  1:2013199/1‑100 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                       >  1:2449862/1‑100 (MQ=255)
                                                         cGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttt                                    <  2:1939073/100‑1 (MQ=255)
                                                          gTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCACAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACtttt                                   >  1:1716134/1‑100 (MQ=255)
                                                            tttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTa                                 >  1:1045967/1‑100 (MQ=255)
                                                            tttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTa                                 <  2:750336/100‑1 (MQ=255)
                                                            tttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCACAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTa                                 <  2:1716122/100‑1 (MQ=255)
                                                            ttcaaTTCTTGCAGCTGCTGCCCGATACGGTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCACTCAATACTACACttttt                                  <  2:1354418/96‑1 (MQ=255)
                                                             ttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAg                                >  1:767107/1‑100 (MQ=255)
                                                              tAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                               <  2:1541349/100‑1 (MQ=255)
                                                              tAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                               >  1:1400001/1‑100 (MQ=255)
                                                              tAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGGTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                               <  2:1213092/100‑1 (MQ=255)
                                                               aaTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCa                              >  2:1515012/1‑100 (MQ=255)
                                                                aTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                             <  1:700191/100‑1 (MQ=255)
                                                                aTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                             <  1:2510226/100‑1 (MQ=255)
                                                                  tCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCagag                           <  2:2341847/100‑1 (MQ=255)
                                                                   cTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAgaga                          >  1:714327/1‑100 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                        >  1:2198924/1‑100 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                        >  2:2177006/1‑100 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                        <  2:953823/100‑1 (MQ=255)
                                                                     tGCAGCAGCTGCCCCATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                        <  2:2283815/100‑1 (MQ=255)
                                                                       cagcagCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAg                      <  1:526212/100‑1 (MQ=255)
                                                                       cagcagCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAg                      <  2:2442956/100‑1 (MQ=255)
                                                                          cagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCa                   <  2:1166410/100‑1 (MQ=255)
                                                                               gCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAc              >  2:2347980/1‑100 (MQ=255)
                                                                                 ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAccc            <  1:1235397/100‑1 (MQ=255)
                                                                                 ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAccc            >  2:872190/1‑100 (MQ=255)
                                                                                  cGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCa           >  2:1718804/1‑100 (MQ=255)
                                                                                          tctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGc   <  1:1111584/100‑1 (MQ=255)
                                                                                          tctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGc   <  2:2387091/100‑1 (MQ=255)
                                                                                           ctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCg  <  1:82587/100‑1 (MQ=255)
                                                                                                 |                                                                                             
GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCG  >  NC_000913/882625‑882815

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