Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 882,722 G→A 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,7220GA97.7% 140.4 / ‑4.0 43R5C (CGC→TGC) deoRdeoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (19/23);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGC  >  NC_000913/882625‑882814
                                                                                                 |                                                                                            
gTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAAc                                                                                            <  1:774988/100‑1 (MQ=255)
 tCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACg                                                                                           <  1:1181671/100‑1 (MQ=255)
 tCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACg                                                                                           <  1:98665/100‑1 (MQ=255)
   aTCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgtg                                                                                         >  2:1284544/1‑100 (MQ=255)
      tcCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGttt                                                                                      >  2:992382/1‑100 (MQ=255)
       ccGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                     >  1:1244776/1‑100 (MQ=255)
       ccGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                     <  2:1238466/100‑1 (MQ=255)
        cGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTAATCGCTGCGTTTTAATTCTTTCATCAGATGCCCGATACGCTCTTCGCAACGTGTTTcc                                                                                    <  2:743632/100‑1 (MQ=255)
          aaaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCAt                                                                                  <  2:515832/100‑1 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCTCGATACGCTCTTCGCAAGGTGTTTCCATagt                                                                               <  1:1063630/100‑3 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                        <  2:666597/100‑1 (MQ=255)
                     ngggCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGa                                                                       >  1:439262/2‑100 (MQ=255)
                        gcggcggcgTCTTTAAGATGTAACGTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAAt                                                                     <  2:73791/99‑1 (MQ=255)
                         cggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAg                                                                   <  1:1323813/100‑1 (MQ=255)
                           gcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGtt                                                                 >  1:1208100/1‑100 (MQ=255)
                            cggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTa                                                                <  2:78804/100‑1 (MQ=255)
                              gcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAtt                                                              <  1:205446/100‑1 (MQ=255)
                               cgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTg                                                             >  1:109756/1‑100 (MQ=255)
                                  cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                          <  2:258587/100‑1 (MQ=255)
                                           gTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCt                                                 <  2:1353337/100‑1 (MQ=255)
                                           gTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCt                                                 <  2:552076/100‑1 (MQ=255)
                                            tAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTc                                                >  1:943058/1‑100 (MQ=255)
                                                 tATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtac                                           >  1:1561395/1‑100 (MQ=255)
                                                   tCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtacta                                         >  1:276623/1‑100 (MQ=255)
                                                    cGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                        <  2:429252/100‑1 (MQ=255)
                                                    cGCTGCGTTTTAATTCTTCCAGCAGCTGCCCGACACGCTCTTCGCAACGTGTTTCCATAATCCCCCAGAATAGCTATTAAAGCGAGCCGCTCAAcactac                                        >  2:66803/1‑100 (MQ=255)
                                                     gCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTaca                                       >  1:1524544/1‑100 (MQ=255)
                                                           ttttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATCCTACACtttt                                  >  1:1432545/1‑99 (MQ=255)
                                                            tttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCACAATACTACACTTTTTa                                >  2:17911/1‑100 (MQ=255)
                                                               aaTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCa                             >  1:710714/1‑100 (MQ=255)
                                                                aTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                            <  2:1100464/100‑1 (MQ=255)
                                                                  tCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCagag                          <  2:1561385/100‑1 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                       <  1:940731/100‑1 (MQ=255)
                                                                      gcagcaGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGTagagtt                        >  2:1030747/1‑96 (MQ=255)
                                                                          cagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCa                  <  2:378066/100‑1 (MQ=255)
                                                                           agcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCAc                 >  1:216350/1‑100 (MQ=255)
                                                                               gCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAc             >  2:1531242/1‑100 (MQ=255)
                                                                                cccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAcc            >  2:1357487/1‑100 (MQ=255)
                                                                                cccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAcc            <  2:1244771/100‑1 (MQ=255)
                                                                                 ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAccc           >  2:1449920/1‑100 (MQ=255)
                                                                                 ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAccc           >  1:168690/1‑100 (MQ=255)
                                                                                      aCGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCt      <  1:385015/100‑1 (MQ=255)
                                                                                          tctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGc  <  1:1396105/100‑1 (MQ=255)
                                                                                                 |                                                                                            
GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGC  >  NC_000913/882625‑882814

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: