Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,169,384 G→A 7.4% intergenic (+54/+28) 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,3840GA 7.4% 138.7 / 3.2 54intergenic (+54/+28)bhsA/ldtCDUF1471 domain‑containing multiple stress resistance outer membrane protein BhsA/L,D‑transpeptidase YcfS
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base A (2/2);  total (26/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.35e-01

AGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGA  >  NC_000913/1169243‑1169525
                                                                                                                                             |                                                                                                                                               
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:1056212/143‑1 (MQ=255)
                                                             cATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGa                                                                                 <  1:233631/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)
                                                                                                                                            cGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGa  <  1:1005355/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
AGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGG‑‑GCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGA  >  NC_000913/1169243‑1169525

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: