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)
               ttACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGttt                                                                                                                                    <  1:525851/150‑1 (MQ=255)
               ttACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTtg                                                                                                                                   >  1:449610/1‑151 (MQ=255)
                    aaCTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTgatgat                                                                                                                              >  1:325279/1‑151 (MQ=255)
                            gTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCa                                                                                                                      <  1:757777/151‑1 (MQ=255)
                               tCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTTCCACCGTAATTAAGTAGCAGCGGTGACAAGTATTTTAATCTACGTAACTATTCGGAATATTTAGAATATACTAGCGGTTTTTTTTTTTAGTTTGATGATTTTTTTCaaa                                                                                                                    >  2:700767/1‑150 (MQ=255)
                                     aataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATTTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGa                                                                                                              <  2:194456/150‑1 (MQ=255)
                                     aataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGa                                                                                                              <  2:539728/150‑1 (MQ=255)
                                      ataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACt                                                                                                            >  2:375693/1‑151 (MQ=255)
                                       taCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtt                                                                                                           >  2:404379/1‑151 (MQ=255)
                                           tGGGTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtttttt                                                                                                       >  1:84790/1‑151 (MQ=255)
                                           tGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtttttt                                                                                                       >  2:184573/1‑151 (MQ=255)
                                             gCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTcc                                                                                                     <  1:18844/151‑1 (MQ=255)
                                                gATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCt                                                                                                    <  2:838070/149‑1 (MQ=255)
                                                gATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTg                                                                                                   >  2:840473/1‑150 (MQ=255)
                                                   tATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaa                                                                                                <  2:884955/150‑1 (MQ=255)
                                                      gCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGt                                                                                             >  1:475244/1‑150 (MQ=255)
                                                       cGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTc                                                                                            >  2:298199/1‑150 (MQ=255)
                                                        ggAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTc                                                                                            >  1:218768/1‑149 (MQ=255)
                                                         gAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGt                                                                                         <  1:401589/151‑1 (MQ=255)
                                                              aaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGat                                                                                     <  1:185693/151‑1 (MQ=255)
                                                              aaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGat                                                                                     <  1:886423/150‑1 (MQ=255)
                                                              aaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑‑TTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGatat                                                                                   <  2:900184/151‑1 (MQ=255)
                                                              aaataaTTTCCCACCGGAATTAAGTAGCAGCCGTGGAAAGTATTTTCATCTACGTAACTAATCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGata                                                                                    <  2:155370/151‑1 (MQ=255)
                                                                         cACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCaa                                                                         <  2:16246/151‑1 (MQ=255)
                                                                         cACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCaa                                                                         <  1:473208/151‑1 (MQ=255)
                                                                            cGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGt                                                                       >  1:494253/1‑150 (MQ=255)
                                                                              gAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAg                                                                    <  1:22559/151‑1 (MQ=255)
                                                                                   aaGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGGGTTTGATGATTTTTATAAAATAATAAATTTTTTCCTGAAAAGTAGGTAGTTATTCAAAGGAAAGTAAGGttt                                                                >  1:417917/1‑150 (MQ=255)
                                                                                     gTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTtat                                                             <  2:449060/151‑1 (MQ=255)
                                                                                          agcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACttt                                                        <  1:758278/151‑1 (MQ=255)
                                                                                                       aTTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTTCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGAtt                                            >  2:314283/1‑150 (MQ=255)
                                                                                                       aTTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTc                                           <  1:720574/151‑1 (MQ=255)
                                                                                                        ttttTATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTTGTGTTTTATGGTTTTTTTCAAATCATGCCTTTTTTCCTTAAAAGTCAGTAGATATTCATAGGCAAGAAAGGCTTTATACTTTTCTGACAGGAtt                                            >  2:545623/1‑150 (MQ=255)
                                                                                                          ttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAg                                         <  1:559152/150‑1 (MQ=255)
                                                                                                              tCTACGTAACTATTCAGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTg                                    >  2:385852/1‑151 (MQ=255)
                                                                                                               cTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGt                                   >  2:882792/1‑151 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtct                                 >  1:461764/1‑151 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtc                                  >  1:270363/1‑150 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGt                                   <  2:494204/149‑1 (MQ=255)
                                                                                                                  cGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtctc                                <  1:679738/151‑1 (MQ=255)
                                                                                                                   gTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCa                               <  2:690021/151‑1 (MQ=255)
                                                                                                                    tAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATATGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCa                               <  2:482860/150‑1 (MQ=255)
                                                                                                                    tAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTTGAGCTTGATGATTTTTATCAAATAATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATCGGCAAGTAAGGGTTTATACTTTGCTGACAGGATTCGGGCATGtacc                                >  1:218996/1‑146 (MQ=255)
                                                                                                                                  atTCAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGt                <  1:370438/151‑1 (MQ=255)
                                                                                                                                     cAGAATATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAAt             >  1:823527/1‑151 (MQ=255)
                                                                                                                                        aaTATACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGa           >  1:218285/1‑150 (MQ=255)
                                                                                                                                         atatACTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAAc         <  2:420068/151‑1 (MQ=255)
                                                                                                                                            taCTAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACaa       >  2:775075/1‑150 (MQ=255)
                                                                                                                                               tAGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTTAAAAGTCATTCGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACACAAgg   >  2:741187/1‑151 (MQ=255)
                                                                                                                                                aGCGGTTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAggg  >  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: