Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,556,769 (T)9→10 intergenic (‑449/‑224) yjiC ← / → iraD uncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,556,7601.T98.0% 189.6 / ‑2.9 50intergenic (‑440/‑233)yjiC/iraDuncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (25/24);  total (25/25)
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.90e-01

ATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGG  >  NC_000913/4556612‑4556905
                                                                                                                                                     |                                                                                                                                                  
aTCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑‑t                                                                                                                                                  <  1:758139/150‑1 (MQ=255)
      ttATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑‑ttttttt                                                                                                                                            <  2:338346/150‑1 (MQ=255)
      ttATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑‑ttttttt                                                                                                                                            <  2:376397/150‑1 (MQ=255)
        atatTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑‑tttttttttt                                                                                                                                         <  2:704908/151‑2 (MQ=255)
             attTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑‑tttttttttttgagt                                                                                                                                    <  1:177573/151‑7 (MQ=255)
               ttACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGttt                                                                                                                                    <  1:525851/150‑1 (MQ=255)
               ttACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTtg                                                                                                                                   >  1:449610/1‑151 (MQ=255)
                    aaCTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTgatgat                                                                                                                              >  1:325279/1‑151 (MQ=255)
                            gTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCa                                                                                                                      <  1:757777/151‑1 (MQ=255)
                               tCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTTCCACCGTAATTAAGTAGCAGCGGTGACAAGTATTTTAATCTACGTAACTATTCGGAATATTTAGAATATACTAGCGGT‑TTTTTTTTTTAGTTTGATGATTTTTTTCaaa                                                                                                                    >  2:700767/1‑150 (MQ=255)
                                     aataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATTTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGa                                                                                                              <  2:194456/150‑1 (MQ=255)
                                     aataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGa                                                                                                              <  2:539728/150‑1 (MQ=255)
                                      ataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACt                                                                                                            >  2:375693/1‑151 (MQ=255)
                                       taCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtt                                                                                                           >  2:404379/1‑151 (MQ=255)
                                           tGGGTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtttttt                                                                                                       >  1:84790/1‑151 (MQ=255)
                                           tGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtttttt                                                                                                       >  2:184573/1‑151 (MQ=255)
                                             gCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTcc                                                                                                     <  1:18844/151‑1 (MQ=255)
                                                gATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCt                                                                                                    <  2:838070/149‑1 (MQ=255)
                                                gATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTg                                                                                                   >  2:840473/1‑150 (MQ=255)
                                                   tATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaa                                                                                                <  2:884955/150‑1 (MQ=255)
                                                      gCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGt                                                                                             >  1:475244/1‑150 (MQ=255)
                                                       cGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTc                                                                                            >  2:298199/1‑150 (MQ=255)
                                                        ggAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTc                                                                                            >  1:218768/1‑149 (MQ=255)
                                                         gAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGt                                                                                         <  1:401589/151‑1 (MQ=255)
                                                              aaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGat                                                                                     <  1:185693/151‑1 (MQ=255)
                                                              aaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGat                                                                                     <  1:886423/150‑1 (MQ=255)
                                                              aaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGatat                                                                                   <  2:900184/151‑1 (MQ=255)
                                                              aaataaTTTCCCACCGGAATTAAGTAGCAGCCGTGGAAAGTATTTTCATCTACGTAACTAATCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGata                                                                                    <  2:155370/151‑1 (MQ=255)
                                                                         cACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCaa                                                                         <  2:16246/151‑1 (MQ=255)
                                                                         cACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCaa                                                                         <  1:473208/151‑1 (MQ=255)
                                                                            cGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGt                                                                       >  1:494253/1‑150 (MQ=255)
                                                                              gAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAg                                                                    <  1:22559/151‑1 (MQ=255)
                                                                                   aaGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGGGTTTGATGATTTTTATAAAATAATAAATTTTTTCCTGAAAAGTAGGTAGTTATTCAAAGGAAAGTAAGGttt                                                                >  1:417917/1‑150 (MQ=255)
                                                                                     gTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTtat                                                             <  2:449060/151‑1 (MQ=255)
                                                                                          agcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACttt                                                        <  1:758278/151‑1 (MQ=255)
                                                                                                       aTTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTTCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGAtt                                            >  2:314283/1‑150 (MQ=255)
                                                                                                       aTTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTc                                           <  1:720574/151‑1 (MQ=255)
                                                                                                        ttttTATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTTGTGTTTTATGGTTTTTTTCAAATCATGCCTTTTTTCCTTAAAAGTCAGTAGATATTCATAGGCAAGAAAGGCTTTATACTTTTCTGACAGGAtt                                            >  2:545623/1‑150 (MQ=255)
                                                                                                          ttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAg                                         <  1:559152/150‑1 (MQ=255)
                                                                                                              tCTACGTAACTATTCAGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTg                                    >  2:385852/1‑151 (MQ=255)
                                                                                                               cTACGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGt                                   >  2:882792/1‑151 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtct                                 >  1:461764/1‑151 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtc                                  >  1:270363/1‑150 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGt                                   <  2:494204/149‑1 (MQ=255)
                                                                                                                  cGTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtctc                                <  1:679738/151‑1 (MQ=255)
                                                                                                                   gTAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCa                               <  2:690021/151‑1 (MQ=255)
                                                                                                                    tAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATATGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCa                               <  2:482860/150‑1 (MQ=255)
                                                                                                                    tAACTATTCGGAATATTCAGAATATACTAGCGGT‑TTTTTTTTTGAGCTTGATGATTTTTATCAAATAATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATCGGCAAGTAAGGGTTTATACTTTGCTGACAGGATTCGGGCATGtacc                                >  1:218996/1‑146 (MQ=255)
                                                                                                                                  atTCAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGt                <  1:370438/151‑1 (MQ=255)
                                                                                                                                     cAGAATATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAAt             >  1:823527/1‑151 (MQ=255)
                                                                                                                                        aaTATACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGa           >  1:218285/1‑150 (MQ=255)
                                                                                                                                         atatACTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAAc         <  2:420068/151‑1 (MQ=255)
                                                                                                                                            taCTAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACaa       >  2:775075/1‑150 (MQ=255)
                                                                                                                                               tAGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTTAAAAGTCATTCGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACACAAgg   >  2:741187/1‑151 (MQ=255)
                                                                                                                                                aGCGGT‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAggg  >  1:59980/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGG  >  NC_000913/4556612‑4556905

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