Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 837,642 A→T 8.4% intergenic (+206/+23) ybiC → / ← ybiJ putative dehydrogenase/DUF1471 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913837,6420AT 8.4% 206.8 / 4.0 59intergenic (+206/+23)ybiC/ybiJputative dehydrogenase/DUF1471 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base T (3/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAG  >  NC_000913/837494‑837784
                                                                                                                                                    |                                                                                                                                              
ttAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTa                                                                                                                                                >  4:601894/1‑149 (MQ=255)
        tatGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGc                                                                                                                                           <  4:1742924/146‑1 (MQ=255)
        tatGCGTCGGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGc                                                                                                                                           >  3:1742924/1‑146 (MQ=255)
                  tCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTa                                                                                                                              >  3:1765536/1‑149 (MQ=255)
                     gCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGa                                                                                                                           >  1:32014/1‑149 (MQ=255)
                           agacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTaaa                                                                                                                                              >  3:1042365/1‑124 (MQ=255)
                           agacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTaaa                                                                                                                                              <  4:1042365/124‑1 (MQ=255)
                           agacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGcc                                                                                                                                          <  4:1870828/128‑1 (MQ=255)
                           agacagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGcc                                                                                                                                          >  3:1870828/1‑128 (MQ=255)
                              cagGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACaaaatcgtc                                                                                                                    >  3:553413/1‑141 (MQ=255)
                                  gAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                         >  4:342490/1‑122 (MQ=255)
                                  gAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                         <  3:342490/122‑1 (MQ=255)
                                   aaaTACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAAtt                                                                                                             >  4:1558680/1‑149 (MQ=255)
                                         ttGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                       <  3:601894/149‑1 (MQ=255)
                                            tgCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAAGGCCTTATCAGTCCTACGAATTACTTATAAATTACCGCa                                                                                                       >  4:1137713/1‑146 (MQ=255)
                                            tgCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                       <  3:1137713/146‑1 (MQ=255)
                                                 gACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACtt                                                                                                                    >  1:448133/1‑128 (MQ=255)
                                                 gACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACtt                                                                                                                    <  2:448133/128‑1 (MQ=255)
                                                    gTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTaa                                                                                            <  1:870805/149‑1 (MQ=255)
                                                         gTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAt                                                                                       >  3:1008179/1‑149 (MQ=255)
                                                         gTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAt                                                                                       >  1:1698588/1‑149 (MQ=255)
                                                          tAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                       <  4:1557844/132‑1 (MQ=255)
                                                          tAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                       >  3:1557844/1‑132 (MQ=255)
                                                          tAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAtt                                                                                      <  3:1558680/149‑1 (MQ=255)
                                                            aaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAAt                                                                                                                         <  3:713180/112‑1 (MQ=255)
                                                            aaCGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAAt                                                                                                                         >  4:713180/1‑112 (MQ=255)
                                                               gCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgt                                                                                 <  4:1765536/149‑1 (MQ=255)
                                                               gCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAtt                                                                                      <  4:1713630/144‑1 (MQ=255)
                                                               gCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAtt                                                                                      >  3:1713630/1‑144 (MQ=255)
                                                                 cTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                <  1:1136955/148‑1 (MQ=255)
                                                                 cTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTACTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                >  2:1136955/1‑148 (MQ=255)
                                                                    aTCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggtg                                                                            >  1:611797/1‑149 (MQ=255)
                                                                      ccGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGc                                                                          <  2:611797/149‑1 (MQ=255)
                                                                           cAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAg                                                                     <  2:1698588/149‑1 (MQ=255)
                                                                             aCGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAg                                                                     >  4:69533/1‑147 (MQ=255)
                                                                             aCGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAg                                                                     <  3:69533/147‑1 (MQ=255)
                                                                              cGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGa                                                                  <  4:1008179/149‑1 (MQ=255)
                                                                                     gATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGt                                                            <  3:317770/148‑1 (MQ=255)
                                                                                     gATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGt                                                            >  4:317770/1‑148 (MQ=255)
                                                                                                   gATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCgg                                             <  2:32014/149‑1 (MQ=255)
                                                                                                     tAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGaa                                                                                                                          <  3:557298/70‑1 (MQ=255)
                                                                                                     tAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGaa                                                                                                                          <  3:829829/70‑1 (MQ=255)
                                                                                                     tAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGaa                                                                                                                          >  4:557298/1‑70 (MQ=255)
                                                                                                     tAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGaa                                                                                                                          >  4:829829/1‑70 (MQ=255)
                                                                                                        gACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTCCGAATTa                                                                                                                       <  2:1018883/70‑1 (MQ=255)
                                                                                                        gACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTa                                                                                                                       >  1:1018883/1‑70 (MQ=255)
                                                                                                                            ggCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGttt                    >  1:219188/1‑149 (MQ=255)
                                                                                                                              cATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggt                                                                             >  4:346649/1‑90 (MQ=255)
                                                                                                                              cATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggt                                                                             <  3:346649/90‑1 (MQ=255)
                                                                                                                                      gacGGATAAGGCGTTTACGCCGCATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTtggtg                                                                            >  4:1808219/3‑83 (MQ=25)
                                                                                                                                      gCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAgg          >  1:692079/1‑149 (MQ=255)
                                                                                                                                            taaggcgTTTACGCCGCATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCt                                           >  3:1361306/4‑110 (MQ=255)
                                                                                                                                            taaggcgTTTACGCCGCATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCt                                           <  2:1456583/107‑1 (MQ=255)
                                                                                                                                            taaggcgTTTACGCCGCATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCt                                           >  1:1456583/4‑110 (MQ=255)
                                                                                                                                            taaggcgTTTACGCCGCATCCGTCCTACGAAGTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCt                                           <  4:1361306/107‑1 (MQ=255)
                                                                                                                                              cggcgTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAg  >  2:1393721/1‑149 (MQ=255)
                                                                                                                                                   tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                >  1:223455/1‑66 (MQ=255)
                                                                                                                                                   tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                >  1:1705882/1‑66 (MQ=255)
                                                                                                                                                   tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                >  1:78033/1‑66 (MQ=255)
                                                                                                                                                   tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                <  2:1705882/66‑1 (MQ=255)
                                                                                                                                                   tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                <  2:223455/66‑1 (MQ=255)
                                                                                                                                                   tAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTAttgttgtt                                                                                <  2:78033/66‑1 (MQ=255)
                                                                                                                                                    aaaCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                       <  4:550535/42‑1 (MQ=255)
                                                                                                                                                    aaaCGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCa                                                                                                       >  3:550535/1‑42 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TTAGAGCATATGCGTCTGTCGGCAAACAGACAGGGAAATACTTGTGCTGGACGTAGCGTAAACGCCTGATCCGTCCAACGGTTCGGATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGCATCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGAATTACTTATAAATTACCGCAGTACCGCTTAATTTATTGTTGTTGGTGGCGGAAGTGATGCTGTATCCACTGGCACCGGCTGCCGCGGCTTTCTCAGCCAGTTTCGCTTCCAGGGCATCGAG  >  NC_000913/837494‑837784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑