Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,138 (A)6→7 coding (15/594 nt) bdcR → transcriptional repressor for divergent bdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,1321.A98.5% 214.2 / ‑3.9 67coding (9/594 nt)bdcRtranscriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (34/32);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

ACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGT  >  NC_000913/4474006‑4474264
                                                                                                                               |                                                                                                                                     
aCGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAGAAAACAAtc                                                                                                                           >  1:1108538/1‑139 (MQ=255)
  gACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAAtctc                                                                                                                         <  2:1435001/139‑1 (MQ=255)
  gACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAAtctc                                                                                                                         <  1:1170744/139‑1 (MQ=255)
       ccACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGtt                                                                                                                    >  2:1598957/1‑139 (MQ=255)
       ccACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGtt                                                                                                                    >  1:713387/1‑139 (MQ=255)
           cGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAg                                                                                                                <  1:1338579/139‑1 (MQ=255)
             aGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              <  2:1108538/139‑1 (MQ=255)
                           cTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                <  2:713387/139‑1 (MQ=255)
                              aCCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGc                                                                                             >  1:603050/1‑139 (MQ=255)
                                   tAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTg                                                                                        <  1:1598957/139‑1 (MQ=255)
                                     aaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:237725/139‑1 (MQ=255)
                                             cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  2:974361/1‑131 (MQ=255)
                                             cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:974361/131‑1 (MQ=255)
                                                      tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATTGTCACTAAA‑AAAACAATCTcg                                                                                                                        <  2:1041125/88‑1 (MQ=255)
                                                      tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTcg                                                                                                                        >  1:1041125/1‑88 (MQ=255)
                                                      tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGc                                                                     >  2:1468367/1‑139 (MQ=255)
                                                      tCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGc                                                                     >  2:624663/1‑139 (MQ=255)
                                                        cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAgg                                                                                                               >  1:757313/1‑95 (MQ=255)
                                                        cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAgg                                                                                                               <  2:757313/95‑1 (MQ=255)
                                                        cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTg                                                                                        <  2:1668712/118‑1 (MQ=255)
                                                        cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTg                                                                                        >  1:1668712/1‑118 (MQ=255)
                                                              ttattaTGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCg                                                                       >  1:65006/1‑129 (MQ=255)
                                                              ttattaTGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCg                                                                       <  2:65006/129‑1 (MQ=255)
                                                                attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCtt                                                           >  1:1345456/1‑139 (MQ=255)
                                                                   atGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttttt                                                        >  2:326055/1‑139 (MQ=255)
                                                                          aTGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 <  2:171073/139‑1 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              <  2:785940/77‑1 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              >  1:785940/1‑77 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAGAGTGCTTTTTCACCaaaa                                                >  1:1167482/1‑139 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttt                                                          >  2:1562865/1‑129 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCttt                                                          <  1:1562865/129‑1 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                >  1:840910/1‑139 (MQ=255)
                                                                           tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                >  1:1427860/1‑139 (MQ=255)
                                                                              aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTcg                                                                                                                        <  1:943048/64‑1 (MQ=255)
                                                                              aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTcg                                                                                                                        >  2:943048/1‑64 (MQ=255)
                                                                                tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGtt                                           <  2:1345456/139‑1 (MQ=255)
                                                                                  cAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGc                                                                        >  2:1712759/1‑108 (MQ=255)
                                                                                  cAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGc                                                                        <  1:1712759/108‑1 (MQ=255)
                                                                                     aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      >  2:54704/1‑91 (MQ=255)
                                                                                     aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                      <  1:54704/91‑1 (MQ=255)
                                                                                       cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                           <  1:1276957/84‑1 (MQ=255)
                                                                                       cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAAGCCCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                           >  2:1276957/1‑84 (MQ=255)
                                                                                          ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                 >  1:1636107/1‑139 (MQ=255)
                                                                                          ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                 <  2:1427860/139‑1 (MQ=255)
                                                                                               aCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              <  1:173657/57‑1 (MQ=255)
                                                                                               aCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGgt                                                                                                              >  2:173657/1‑57 (MQ=255)
                                                                                               aCTTTCTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAgtg                            >  2:309948/1‑139 (MQ=255)
                                                                                                    cTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 <  2:147303/113‑1 (MQ=255)
                                                                                                    cTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                 >  1:147303/1‑113 (MQ=255)
                                                                                                    cTGTCAACAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCt                       >  1:129544/1‑139 (MQ=255)
                                                                                                         aaCAGGTTTTGTGATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                     <  1:759154/56‑1 (MQ=255)
                                                                                                         aaCAGGTTTTGTGATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                     >  2:759154/1‑56 (MQ=255)
                                                                                                         aaCAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCa                                                                                >  1:1039639/1‑77 (MQ=255)
                                                                                                         aaCAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCa                                                                                <  2:1039639/77‑1 (MQ=255)
                                                                                                           cAGGTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa                <  2:129544/139‑1 (MQ=255)
                                                                                                              gTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGc                        <  1:1020848/128‑1 (MQ=255)
                                                                                                              gTTTTGTAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTCTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGc                        >  2:1020848/1‑128 (MQ=255)
                                                                                                                    tAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa            <  2:563824/134‑1 (MQ=255)
                                                                                                                    tAATGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa            >  1:563824/1‑134 (MQ=255)
                                                                                                                     aaTGGTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt      <  2:603050/139‑1 (MQ=255)
                                                                                                                        ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                <  2:226609/46‑1 (MQ=255)
                                                                                                                        ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                >  1:226609/1‑46 (MQ=255)
                                                                                                                        ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAg                              >  2:1690782/1‑112 (MQ=255)
                                                                                                                        ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAg                              <  1:1690782/112‑1 (MQ=255)
                                                                                                                        ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa            <  2:170476/130‑1 (MQ=255)
                                                                                                                        ggTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa            >  1:170476/1‑130 (MQ=255)
                                                                                                                         gTCACTAAA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGt  >  2:824595/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                                     
ACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AA‑AAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGT  >  NC_000913/4474006‑4474264

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