Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,380 T→C 8.6% intergenic (+50/+32) bhsA → / ← ldtC DUF1471 domain‑containing multiple stress resistance outer membrane protein BhsA/L,D‑transpeptidase YcfS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,169,3800TC 8.6% 134.7 / 3.5 53intergenic (+50/+32)bhsA/ldtCDUF1471 domain‑containing multiple stress resistance outer membrane protein BhsA/L,D‑transpeptidase YcfS
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (2/3);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

AGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAA  >  NC_000913/1169243‑1169522
                                                                                                                                         |                                                                                                                                                
aGCGGATTAGATGGGCGCAAAATCTTTCCGTATTACTTCTCGAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATGCGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCgt                                                                                                                                             <  1:772042/143‑1 (MQ=255)
  cGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACTTGTCGGAATCGCCCCCATGGGGTAGCCCGTATCAAGAAACAAGCATTAACCCTCATTCATGCCTCCTACTGCTGATTTTTTCCCCGCGACCTGCCgtgt                                                                                                                                           >  1:524302/1‑143 (MQ=255)
       gggaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                      <  2:814759/141‑1 (MQ=255)
       gagaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTTCATGGAACATCAGTAATTTATAAATAAGCATTAACCCGCATCAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg                                                                                                                                      <  1:540071/143‑1 (MQ=255)
       gagaTGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCCGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg                                                                                                                                      <  1:865957/143‑1 (MQ=255)
              gcgcAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAtttt                                                                                                                               <  1:91349/143‑1 (MQ=255)
               cgcAAAATCTTTCCGTCTTATTTCTGTAACCGGTCCGAATTCCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                              <  2:813426/143‑1 (MQ=255)
               cgcAAAATCTTTCCGTATTACTTCCGTAACCGGTCCGAATACCCCCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCAGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAttttt                                                                                                                              <  2:815717/143‑1 (MQ=255)
                gcAAAATCTTGCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCACGGGGAtttttt                                                                                                                             <  2:509058/143‑1 (MQ=255)
                    aaTCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGCTTTTTTtctc                                                                                                                         >  1:697412/1‑143 (MQ=255)
                        tttCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCggg                                                                                                                     <  1:698780/143‑1 (MQ=255)
                           ccGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCtt                                                                                                                  <  1:873342/143‑1 (MQ=255)
                             gTATTACTTCTGTAACCGGTCAGAATACCCTCCATGGAATAGCAGTAATTTATAAATAAGCATTAACACTCATTAATGCCTGCTACTGCTGATAATTCCCGCGCGTCAGGCCGTGTCGCCCGGGTTTTATTCGCCGGGcatca                                                                                                                >  2:798433/1‑139 (MQ=255)
                              tATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTAc                                                                                                               >  1:735338/1‑143 (MQ=255)
                                ttACTGCAGTAACCGGGCCGAATACCCTACATGGAACAGAAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCTCGACATGCCGTGTCGCGGCGATTTTTTTATCCGGGCTTTACAg                                                                                                             <  1:653640/143‑1 (MQ=255)
                                    ttCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAGATAAGCATTAACGCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCTGGGCTTTattacgtt                                                                                                         >  1:1076093/1‑136 (MQ=255)
                                         tAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGAATTAACCCTCATTAATGCCTGCTCCTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTGTTTATCCGGGCTTTACAGCGTTTGTgg                                                                                                    >  2:375860/1‑143 (MQ=255)
                                                tCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACt                                                                                           >  2:1005355/1‑143 (MQ=255)
                                                 ccGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGGGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGGTTGTGG‑‑GCTCACtt                                                                                          >  1:1086634/1‑143 (MQ=255)
                                                         ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGAcgg                                                                                  >  1:538810/1‑143 (MQ=255)
                                                         ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGAcgg                                                                                  >  1:597633/1‑143 (MQ=255)
                                                          cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGAcggc                                                                                 >  2:653127/1‑143 (MQ=255)
                                                             cATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGa                                                                              <  1:233631/143‑1 (MQ=255)
                                                             cATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGa                                                                              <  1:1056212/143‑1 (MQ=255)
                                                              ctGGAACAGCAGTAATTTATAAATAAGCAGTAGCCCTCATTCATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGAc                                                                             <  2:394134/142‑1 (MQ=255)
                                                                      gcagTAATTTCTAGTTACGCATTAGCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCTCTTGATGACGGCGGACATCCACCg                                                                     <  2:526881/143‑1 (MQ=255)
                                                                         gTAATTTATAAATAGGCATTACCCCTCTTTAATGCCTGCTACTGCTGTTTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCa                                                                  <  2:289989/143‑1 (MQ=255)
                                                                          tAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACGGCGTTTGTGG‑‑GCTCTCTTGCTGAGGGCCGGCGTCCATCggtct                                                                 >  2:695023/1‑140 (MQ=255)
                                                                                  aaataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGTTGACGGCGGACATCCACCGGCATCCCGGAAc                                                         >  1:881609/1‑143 (MQ=255)
                                                                                      aaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGAATTTTTTATCCGGGCTTTACAGCGTTTGAGGAAGAACACTTCTGCGGCGGACAAGCCCCGGCATCCCGGAACggcc                                                     >  1:818673/1‑141 (MQ=255)
                                                                                         cactaccccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGGCGCGGGGAGTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatc                                                  <  1:447240/137‑1 (MQ=255)
                                                                                         cATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGGCGGCGGGCCTCCACCGGCATCCCGGAACGGCcatc                                                  >  1:716298/1‑143 (MQ=255)
                                                                                           ttAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTGTTTATCCTGGCTTTACAGCGTTTGTGG‑‑GGTGACGTGATGACGGCGGCCATCGACCGGCACCCAGGAACGGAcgacca                                                >  2:392271/1‑138 (MQ=255)
                                                                                             agcccTCATTAATGGCTGCTCCTGCTGTTTTTTTCCCCGCGACATACCGTGTCGCGGGTATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccatc                                              <  2:101382/141‑1 (MQ=255)
                                                                                                ccTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATcaca                                           <  1:1057772/143‑1 (MQ=255)
                                                                                                    gttAATGCCTGCTACTGCTGATTATTTCCCGGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTGCAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGtt                                       <  2:1130333/142‑1 (MQ=255)
                                                                                                       aaTGCCTGCTACTGCTGATGTTTGCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCa                                    <  1:538195/143‑1 (MQ=255)
                                                                                                        aTGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCAt                                   <  1:707249/143‑1 (MQ=255)
                                                                                                            cTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACGTGTTGCGTCCCt                               >  2:157120/1‑143 (MQ=255)
                                                                                                              gCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGTCGGACATCCACCGGCATCCCGGAACGGAAACACATCACATGTTGCATTACTTc                             >  1:641678/1‑143 (MQ=255)
                                                                                                              gCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTc                             >  1:631110/1‑143 (MQ=255)
                                                                                                               cTACTGCTGTGTTTTTTCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTGACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCa                            <  2:717576/143‑1 (MQ=255)
                                                                                                                     ctgATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTATTATCCCGGATTTGGGTCGATTGTGG‑‑GCTCATACGCTGACGGCGGACATCCACAGGCATCCCGGTCCCGCCATCCATCACATGTTGCATCACTTCAGCGTCa                      >  1:523781/1‑143 (MQ=255)
                                                                                                                       gATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCCCCGGCATCCCGGCACGGACATCCATCACCTGTTGCCTCACTTCAGCGTCAGt                    >  2:236052/1‑143 (MQ=255)
                                                                                                                        aTTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCGCTTGCTGACGGCGGACATCCACCGGCATCGCGGAACGGACATCCATCACATGTTGCATCACTGCACCGTccgct                   >  2:319704/1‑139 (MQ=255)
                                                                                                                        aTTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGGGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCGTCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGtt                   >  1:363485/1‑143 (MQ=255)
                                                                                                                         ttttttCCCCGGGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGttt                  <  2:1075025/143‑1 (MQ=255)
                                                                                                                         tttgttCCCCGCGACACGCCGGGGCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGttt                  <  2:962457/143‑1 (MQ=255)
                                                                                                                                 ccGCGACACGGCATGTCGCGGGGATTTTTTTCTCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGACATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCCGTTTGTGCTGCa          >  1:167212/1‑143 (MQ=255)
                                                                                                                                 ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑ACTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGAGTAAATTAGTGCTGCa          >  1:821824/1‑143 (MQ=255)
                                                                                                                                  cgcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCAt         <  1:26938/143‑1 (MQ=255)
                                                                                                                                  cgcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCAt         >  1:117591/1‑143 (MQ=255)
                                                                                                                                  cgcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACAGCCATCACAGCTTGCATCGCTTCAGCGTCAGTTTGTCCTCCac         >  2:833692/1‑142 (MQ=255)
                                                                                                                                         tGCCGTGTCGCGGGGATTTTTTTAGCCGGGCGTTAGAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa  <  1:327509/143‑1 (MQ=255)
                                                                                                                                         tGCCGGGTCGGGGGGAGTTTTTTATCCGGGTTTTACAGCGGTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa  <  2:315031/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                
AGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAA  >  NC_000913/1169243‑1169522

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